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Comparative  Analyses  of  Aryl  Hydrocarbon  Receptor  Structure  and  Function  in 

Marine  Mammals 


by 

Joy  M.  Lapseritis 


ABSTRACT 

Marine  mammals  possess  high  body  burdens  of  persistent  organic  pollutants, 
including  PCBs  and  dioxin-like  compounds  (DLC).  Chronic  environmental  or 
dietary  exposure  to  these  chemicals  can  disrupt  the  function  of  reproductive  and 
immune  systems,  as  well  as  cause  developmental  defects  in  laboratory  animals. 
The  aryl  hydrocarbon  receptor  (AHR)  is  a  ligand-activated  transcription  factor, 
mediating  the  expression  of  a  suite  of  genes  in  response  to  exposure  to  DLC  and 
structurally  related  chemicals.  Species-specific  differences  in  AHR  structure  can 
affect  an  organism’s  susceptibility  to  the  effects  of  DLC.  The  structures  and 
functions  of  several  cetacean  AHRs  were  investigated  using  in  vitro  molecular 
cloning  and  biochemical  techniques.  Using  a  novel  combination  of  remote 
biopsy  and  molecular  cloning  methods,  RNA  was  extracted  from  small 
integument  samples  from  living  North  Atlantic  right  whales  to  identify  the  cDNA 
sequence  for  AHR  and  other  genes  of  physiological  importance.  Biopsy-derived 
RNA  was  found  to  be  of  higher  quality  than  RNA  extracted  from  stranded 
cetaceans,  and  proved  a  good  source  for  identifying  cDNA  sequences  for 
expressed  genes.  The  molecular  sequences,  binding  constants,  and 
transcriptional  activities  for  North  Atlantic  right  whale  and  humpback  whale  AHRs 
cDNAs  were  determined  using  in  vitro  and  cell  culture  methods.  Whale  AHRs 
are  capable  of  specifically  binding  dioxin  and  initiating  transcription  of  reporter 
genes.  The  properties  of  these  AHRs  were  compared  with  those  from  other 
mammalian  species,  including  human,  mouse,  hamster,  and  guinea  pig,  and 
other  novel  marine  mammal  AHRs,  using  biochemical,  phylogenetic,  and 
homology  modeling  analyses.  The  relative  binding  affinities  for  some  marine 
mammal  AHRs  fall  between  those  for  the  high-affinity  mouse  AHR6'1  and  the 
lower  affinity  human  AHR.  Species-specific  variability  in  two  regions  of  the  AHR 
ligand  binding  domain  were  identified  as  having  the  greatest  potential  impact  on 
AHR  tertiary  structure,  yet  does  not  sufficiently  explain  differences  observed  in 
ligand  binding  assays.  Additional  studies  are  necessary  to  link  exposure  to 
environmental  contaminants  with  potential  reproductive  effects  in  marine 
mammals,  especially  via  interactions  with  steroid  hormone  receptor  pathways. 
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1.  ENVIRONMENTAL  CONTAMINANTS:  PERSISTENT  ORGANIC 
POLLUTANTS 


1. 1  Xenobiotic  Compounds 

The  post-industrial  era  environment  is  burdened  with  a  variety  of 
persistent  organic  pollutants  (POPs),  including  some  of  the  most  toxic  xenobiotic 
substances  known.  Innumerable  xenobiotic  compounds  exist  in  the  environment, 
and  are  defined  simply  as  'chemicals  foreign  to  a  living  organism’.  In  most  cases 
however,  ‘xenobiotics’  refers  specifically  to  anthropogenic  contaminants 
(toxicants),  excluding  toxins  that  may  be  exogenously  produced  by  other 
organisms  or  by  geological  means.  Many  of  the  most  problematic  xenobiotic 
compounds  tend  to  be  persistent  both  in  fatty  tissues  of  organisms  and  in  marine 
sediments  owing  to  their  lipophilic  properties  and  resistance  to  environmental 
and  biologic  degradation.  Of  those  POPs  that  are  subject  to  metabolic 
degradation,  some  can  become  more  potent  (and  dangerous)  than  the  parent 
compound.  2,3,7,8-tetrachlorodibenzo-p-dioxin  (TCDD  or  dioxin)  is  the  most 
toxic  xenobiotic  chemical,  and  is  the  benchmark  compound  for  many 
experimental  studies  in  order  to  classify  the  relative  toxicity  of  structurally  related 
compounds  such  as  PAHs  (polycyclic  aromatic  hydrocarbons),  halogenated 
aromatic  hydrocarbons  (HAHs),  and  planar  HAHs  (PHAHs)  such  as 
polychlorinated  biphenyls  (PCBs)  and  dioxins.  Trophic  transfer  contributes  to 
bioaccumulation  and  biomagnification,  sometimes  leading  to  large  contaminant 
body  burdens  in  wildlife  and  humans. 

1.2.  Sources,  Levels,  &  Fates 

From  post-World  War  II  until  approximately  the  1970s,  the  production  and 
use  of  POPs  was  widespread  in  the  United  States.  Some  POPs  were 
intentionally  manufactured  as  pesticides,  fungicides,  and  herbicides,  whereas 
others  are  by-products  of  industrial  activities,  such  as  combustion,  pulp  mills  and 
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paper  production,  and  electrical  power  plants.  Concerns  about  effects  on  human 
and  wildlife  health  helped  to  restrict  the  use  of  these  compounds  in  the  USA  and 
Europe.  Nevertheless,  the  environmental  persistence  of  even  small  amounts  of 
dioxins,  dibenzofurans,  PCBs,  and  other  organochlorines  remains  problematic, 
revealing  more  of  the  sublethal  effects  of  POPs  in  the  absence  of  extreme 
exposures  (Giesy  et  al.  1994).  Unfortunately,  POP  production  and  use  in  less 
developed  countries  is  less  controlled  and  leads  to  continual  inputs  of  xenobiotic 
contaminants  in  the  global  system.  Global  POP  contamination  is  understood  to 
result  in  mass  transfer  of  airborne  pollutants  to  the  open  oceans,  especially  to 
arctic  and  tropical  regions,  leading  to  the  conclusion  that  the  marine  environment 
is  a  ‘sink1  for  environmental  contaminants  (Ali  et  at.',  MacKay  1982;  Tanabe  et  al. 
1994).  PCBs  exist  in  the  tissues  of  most  marine  animals,  including  prey  species 
such  as  plankton  and  fish  (Ali  et  al Hoekstra  et  al.  2003).  Marine  mammals 
accumulate  high  concentrations  of  POPs  through  trophic  transfer  and 
biomagnification  (Hoekstra  et  al.  2003)  as  well  as  via  lactational  transfer  across 
generations  (Tanabe  et  al.  1994), 

The  majority  of  toxicological  studies  of  marine  mammals  simply  measure 
POP  residue  levels;  however,  direct  comparisons  of  their  impacts  are 
complicated  by  geographic  area  and  species-specific  factors  including  trophic 
position,  metabolic  rates,  and  seasonal  changes  in  blubber  thickness.  High 
tissue  burdens  of  PCBs  and  DDTs  have  been  reported  in  many  odontocete  (up 
to  700  pg/g  sum  blubber  PCBs;  striped  dolphins  (Kannan  et  al.  1993),  pilot 
whales  (Weisbrod  et  al.  2000b),  beluga  whales  (Muir  et  al.  1996),  killer  whales 
(Ross  et  al.  2000))  and  pinniped  species  (up  to  5000  pg/g  sum  blubber  PCBs; 
California  sea  lions  (DeLong  etal.  1973),  grey  seals  (Blomkvist  et  al.  1992)). 
Mysticete  POP  levels  range  from  trace  levels  up  to  27.45  pg/g  PCBs  and  47.70 
pg/g  DDT  (Weisbrod  et  al.  2000a;  Woodley  et  al.  1991).  Analyses  of  blubber 
biopsies  have  revealed  the  presence  of  PCBs,  DDT  and  other  organochlorines  in 
right  whales  (Weisbrod  ef  al.  2000a;  Woodley  et  al.  1991).  Concentrations  of 
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POPs  in  polar  bear  blubber  are  measurable  in  the  mg/kg  range  (Norstrom  and 
Muir  1994),  and  can  reach  incredibly  high  levels  in  the  far  north  and  east  Arctic 
limits  of  their  range  (Norstrom  et  al.  1998). 


1.3.  Potential  Impacts 

Typical  exposures  of  marine  mammals  to  POPs  seldom  result  in  direct 
lethality.  However,  among  the  most  obvious  effects  of  POPs  on  wildlife  is 
reproductive  impairment,  ranging  from  conception  to  early  neonatal  development 
and  to  gonadal  maturation.  Although  these  effects  are  often  variable  and  difficult 
to  compare  due  to  differences  in  exposures  and  physiological  status,  nearly 
every  organism  studied  has  shown  some  form  of  reproductive  dysfunction 
associated  with  elevated  POP  concentration.  Eggshell  thinning,  suppressed 
vitellogenesis,  malformed  reproductive  organs,  reduced  numbers  of  offspring, 
and  complete  failure  to  reproduce  are  documented  in  fish,  birds,  and  mammals 
(Colburn  and  Clement  1992).  Other  potential  impacts  of  POPs  include 
suppression  of  immune  function,  increased  risk  of  tumors,  and  porphyria. 

Exposure  to  POPs  has  been  linked  to  numerous  sublethal  effects  on 
reproduction  and  immune  function  in  marine  mammals  (Reijnders  2003; 

Reijnders  1986a;  Tanabe  et  al.  1994).  Reproductive  impairment  due  to  PCB  and 
DDT  exposure  has  been  suggested  in  beluga  whales  in  the  St.  Lawrence  River 
population  (Martineau  et  al.  1987;  Masse  etal.  1986).  North  Atlantic  right  whales 
(NARW)  exhibit  symptoms  of  reproductive  dysfunction  (low  calving  rate,  high 
neonatal  mortality  (Knowlton  etal.  1994)),  however  efforts  to  connect  POP 
burdens  with  actual  reproductive  impairment  have  been  limited  by  sampling 
techniques  and  other  physiological  factors  (Woodley  et  al.  1991).  Reproductive 
impairment  by  PCBs  and  DDTs  has  been  hypothesized  and  experimentally 
tested  in  pinnipeds  (Helle  et  al.  1976a;  1976b;  Le  Boeuf  and  Bonnell  1971; 
Reijnders  1986b).  Field  observations  of  Baltic  seals  have  correlated  high  blubber 
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PCB  burdens  with  pregnancy  ratios  and  uterine  pathologies  (Helle  et  al.  1976a; 
1976b),  and  high  PCB  levels  in  California  sea  lions  are  associated  with 
premature  births  (DeLong  ef  al.  1973). 

Experimental  studies  of  terrestrial  mammals  provide  important  links 
between  field  observations  in  marine  mammals  and  quantifiable  data.  Mink 
suffer  reproductive  impairment  and  increased  kit  mortality  from  even  low  PCB 
exposures  (Heaton  et  al.  1995;  Resturn  et  at .  1998).  Experimental  work  with 
rhesus  monkeys  exposed  to  PCBs  has  shown  resorption,  abortion,  and  stillbirth 
rates  consistent  with  mink  data  (Barsotti  et  al.  1976).  Fetal  death  and  stillbirth 
can  be  induced  in  sows  treated  with  low  levels  of  PCBs  (Hansen  ef  al.  1975). 

The  low  concentrations  of  PCBs  used  in  these  studies  are  consistent  with 
environmental  levels. 

2.  A  MECHANISM  OF  ACTION 
2. 1.  Aryl  Hydrocarbon  Receptor 

Although  the  exact  mechanisms  of  action  for  these  POPs  are  not  entirely 
clear,  dioxin-like  compounds  proximally  bind  to  the  aryl  hydrocarbon  receptor 
(AHR),  a  soluble,  ligand-activated  transcription  factor,  and  a  protein  in  the  bHLH- 
PAS  (basic  helix-loop-helix  Period-ARNT-Singleminded)  superfamily  causing 
alterations  in  gene  expression  and  toxicity  (Burbach  et  al.  1992;  Poland  et  al. 
1976).  Upon  ligand  binding,  AHR  is  transported  into  the  nucleus  to  dimerize  with 
ARNT  (Ah  Receptor  Nuclear  Translocator).  In  the  nucleus.  TCDD-AHR-ARNT 
activates  transcription  of  genes  with  promoters  containing  xenobiotic  response 
elements  (XREs;  Figure  1 ).  Best  known  for  its  interaction  with  dioxin,  the  AHR 
participates  in  a  number  of  cellular  activities  including  induction  of  cytochrome 
P450  1A  (CYP1A),  cell  cycle  control,  and  teratogenesis,  but  an  endogenous 
ligand  has  not  been  identified.  Because  the  AHR  has  been  predominantly 
defined  by  its  role  in  xenobiotic  metabolism,  yet  is  also  apparently  involved  with 
basic  cell  physiology,  this  mechanism  of  action  provides  an  excellent  framework 
for  integrating  toxicology  with  basic  biology. 
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Studies  of  TCDD  exposure  in  animals  as  diverse  as  chicken,  monkey, 
cow,  hamster,  and  other  vertebrates  have  shown  a  common  pattern  of  toxic 
effects  and  biochemical  induction  implicating  an  AHR-dependent  mechanism  of 
action  (Hahn  1998;  Poland  and  Knutson  1982).  Species-specificity,  as  well  as 
dose-dependent  effects  and  duration  of  exposure,  complicate  direct  comparison 
of  toxic  endpoints,  yet  many  of  the  lesions  observed  in  vertebrates  exposed  to 
TCDD  and  HAH  exhibit  tissue-specific  pathologies  including  hyperplasias  and 
altered  differentiation,  or  hypoplasias,  atrophy,  and  necrosis  of  epithelial  tissues 
(Poland  and  Knutson  1982).  Short  of  the  lethality  observed  from  acute  doses, 
laboratory  animals  display  a  variety  of  short-term  toxic  responses  ranging  from 
wasting  syndrome  to  immune,  hepatic,  and  endocrine  pathologies  (Pohjanvirta 
and  Tuomisto  1994).  Experiments  with  AHR-null  mice  have  shown  that  the  AHR 
is  required  for  TCDD-  and  benzo[a]pyrene-induced  toxicity  (Fernandez-Salguero 
et  al.  1996;  Mimuraefa/,  1997).  These  “knock-out”  experiments  substantiate 
over  thirty  years  of  work  investigating  the  mechanism  action  for  TCDD  and 
related  HAH. 

Studies  of  mouse  and  human  AHRs  have  helped  to  define  the  basic 
sequence  structure  and  specific  functional  domains  for  AHRs  in  general 
(Burbach  et  al.  1992;  Dolwick  et  al.  1993b;  Fukunaga  and  Hankinson  1996; 
Fukunaga  et  al.  1995;  Mimura  ef  al.  1994).  Ten  to  eleven  exons  compose  AHR 
cDNAs  of  approximately  3.0  kb  in  length,  encoding  proteins  of  850  -  950  amino 
acids  and  95-110  kDaltons.  DNA  binding  and  heterodimerization  with  cofactors 
predominantly  occurs  in  the  amino-terminal  portion  of  the  protein,  whereas  the 
more  centrally  located  PAS-B  domain  controls  ligand  binding,  The  carboxyl- 
terminus  tends  to  include  numerous  glutamine  residues  and  is  critical  for 
transcriptional  activity  (Figure  2).  The  bHLH,  PAS-A  and  PAS-B  domains  are 
highly  conserved  across  vertebrate  species,  and  are  easily  recognizable  in  amino 
acid  alignments. 
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3.  SPECIES  OF  INTEREST 


3. 1  Marine  Mammals 

Marine  mammal  toxicology  has  been  limited  by  numerous  logistical  and 
ethical  challenges.  Marine  mammals  are  protected  in  the  United  States  by  the 
Endangered  Species  Act  and  the  Marine  Mammal  Protection  Act,  and  with  a  few 
notable  exceptions,  other  nations  also  provide  some  level  of  protection  for  these 
species.  Protection  is  mandated  due  to  exploitative  commercial  hunting  of  some 
species  to  near  (or  total)  extinction,  in  addition  to  more  contemporary  threats 
such  as  habitat  degradation,  high  mortalities  from  ship  strikes,  epizootic 
episodes,  and  a  plethora  of  health  issues.  Marine  mammals  are  sensitive  targets 
for  POPs  because  of  the  large  proportion  of  lipids  per  body  mass,  trophic 
position,  and  extensive  life  spans  that  characterize  this  group.  Most  of  the  data 
used  for  studying  marine  mammal  biology  and  toxicology  originates  from 
stranded  animals,  skewing  the  perspective  towards  sick  or  unhealthy  individuals. 
Healthy  specimens  can  be  collected  from  subsistence  hunts  (North  America)  and 
scientific  whaling  (Japan),  but  are  a  limited  and  undesirable  means  for  acquiring 
information  on  biological  or  physiological  condition  of  marine  mammals. 
Notwithstanding  these  challenges,  the  dramatic  physiological  adaptations  in 
marine  mammals  to  their  environment  provide  a  unique  perspective  for 
mammalian  toxicology  and  tempt  scientific  investigation.  Recent  advances  in 
molecular  biology  and  acquisition  of  small  biopsies  promise  to  help  negotiate 
these  issues. 

These  logistical,  ethical,  and  legal  challenges  limit  the  methods  available 
to  test  toxicological  questions  in  marine  mammals.  In  vivo  studies  can  provide 
the  most  definitive  results,  but  such  work  is  essentially  impossible  due  to  the 
reasons  stated  above  (with  several  exceptions,  for  example  see  Reijnders 
(1986)).  Less  invasive  environmental  studies  using  stranded  (sick  or  dead) 
marine  mammals  are  common,  but  research  progress  is  constrained  by  the 
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random  numbers,  species,  and  physical  condition  of  stranded  animals.  The 
conclusions  of  these  analyses  must  be  carefully  derived  to  account  for  sampling 
bias.  Live  biopsy  of  healthy,  wild  marine  mammals  (Lambertsen  1987)  has 
provided  a  means  to  reduce  bias  and  improve  sample  quality;  however,  small 
skin  and  blubber  samples  may  not  be  representative  of  total  toxicant  burden  or 
other  factors  being  measured  (Woodley  et  al.  1991). 

Because  cetaceans  have  been  shown  to  retain  more  PCBs  than  other 
animals,  including  lower  chlorinated  congeners  that  are  typically  rapidly 
metabolized,  it  has  been  suggested  that  cetaceans  have  a  reduced  ability  to 
metabolize  certain  PCB  isomers  (Tanabe  et  al.  1988).  Comparisons  of  “marine 
mammal"  biology  are  complicated  in  that  taxonomic  relationships  within  this 
artificial  (habitat-based)  group  are  not  monophyletic,  and  trophic  position  varies 
widely.  Such  variability  among  mammals  can  cause  problems  in  making  valid 
comparisons. 

The  AHR  has  been  characterized  in  an  odontocete  (beluga)  (Jensen  and 
Hahn  2001)  and  two  pinnipeds  (harbor  and  Baikal  seals)  (Kim  and  Hahn  2002; 
Kim  et  al.  2002).  These  studies  have  illustrated  high  conservation  of  amino  acid 
sequence,  and  beluga  and  harbor  seal  AHR  have  shown  specificity  and  relatively 
high  affinity  of  AHR  for  TCDD  in  vitro.  The  extent  of  sequence  and  functional 
conservation  of  AHR  among  marine  mammal  clades,  and  in  comparison  with 
terrestrial  mammals,  has  yet  to  be  examined. 

3,2  Surrogate  Species 

The  use  of  surrogate  species  presents  an  attractive  alternative  to  the 
methodological  challenges  faced  in  addressing  marine  mammal  toxicology. 
Human  studies  -  limited  by  similar  ethical  and  legal  issues  -  have  benefited 
tremendously  from  experimental  work  on  rodents,  for  example.  However  in 
order  to  be  truly  useful,  surrogate  species  must  be  carefully  chosen  and  the 
results  validated  for  applicability  to  the  target  species.  Furthermore,  species- 
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specific  characteristics  (either  of  the  surrogate  or  target  species)  must  be 
examined  and  experimentally  controlled.  Because  mechanisms  of  PCB-related 
reproductive  toxicity  are  not  well  understood  and  species-specific  differences 
may  exist,  it  is  important  to  choose  a  surrogate  species  with  high  sensitivity  to 
conservatively  estimate  toxicant  effects. 

The  mink  has  been  heralded  as  an  ideal  surrogate  species  for  mammals 
because  of  its  extreme  sensitivity  to  POPs  and  ease  of  handling.  Because  mink 
enzymatic  activities  are  more  similar  to  marine  than  terrestrial  mammals,  mink  is 
considered  a  particularly  good  surrogate  for  small  cetacean  species  (Schwacke 
et  al.).  However,  its  applicability  to  these  or  other  marine  mammal  species  in 
regards  to  AHR  sensitivity  has  not  yet  been  validated. 

Another  alternative  is  the  use  of  in  vitro  studies,  which  permit  a 
reductionist  approach  and  perhaps  yield  the  most  discrete  results.  Focusing  on 
the  molecular  and  cellular  level,  in  vitro  studies  include  assays  using  amplified 
genes,  native  extracted  or  recombinant  proteins,  and  experimentation  in  cultured 
cell  lines.  For  example,  in  vitro  protein  interactions  with  potential  ligands  or  other 
proteins  can  identify  functional  characteristics  that  may  be  important  at  the 
organismal  level,  or  in  the  least  can  help  to  formulate  hypotheses  for  future 
study.  Combined  with  careful  inclusion  of  surrogate  species  and  comparison 
with  model  species  (such  as  rodents),  the  results  of  in  vitro  assays  such  as 
ligand  binding  studies  can  be  used  to  infer  in  vivo  conditions. 

3.3  Potential  Risks  Specific  to  North  Atlantic  Right  Whale 

Environmental  contaminants  have  been  measured  in  NARW  blubber 
biopsies,  feces,  and  prey.  Some  authors  caution  that  seasonal  changes  in  lipid 
deposition  and  mobilization  influence  the  concentrations  of  contaminants  by 
effectively  diluting  measurements  during  periods  of  blubber  deposition  (Woodley 
et  al.  1991).  Others  have  shown  fluctuations  based  on  sampling  efforts  over 
seasons  and  years  (Weisbrod  et  al.  2000a).  Persistent  organic  pollutants 
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measured  in  the  most  superficial  blubber  are  a  record  of  recalcitrant  chemicals  or 
metabolic  products,  and  may  not  accurately  represent  the  recent  exposures  of 
NARW  to  environmental  contaminants.  Recent  work  on  bottlenose  dolphins 
shows  that  toxicant  mobilization  is  dynamic  and  more  closely  associated  with  the 
deepest  blubber  strata  (Montie  2006).  Deep  tissue  biopsies  are  more  likely  to 
sample  adipose  tissues  involved  in  the  storage  and  mobilization  of 
organochlorines  in  cetaceans  (Aguilar  1985). 

But  what  is  the  actual  daily  or  annual  exposure  of  an  adult  NARW  to  AHR 
agonists?  Because  profiles  of  blubber  contaminants  mirror  the  organochlorine 
signatures  of  specific  prey  and  regions,  it  is  assumed  that  exposure  to 
environmental  contaminants  is  primarily  dietary  (Weisbrod  et  al.  2000a).  We 
estimated  the  daily  toxicant  exposure  of  NARW  based  on  an  adult  body  mass  of 
40,000  kg  (Kenney  et  al.  1986),  NARW  ingestion  rates  and  feeding  ecology 
(Baumgartner  and  Mate  2003;  Durbin  et  al.  2002),  and  measurements  of  prey 
contaminant  loads  (Moore  et  al.  in  preparation;  Weisbrod  et  al.  2000a).  Durbin  et 
al.  (2002)  estimated  two  ingestion  rates  for  NARW.  Assuming  a  weight-specific 
feeding  rate,  they  calculated  a  high  ingestion  rate  of  9.41  x  10e  copepods  per 
day.  By  tagging  NARW  to  track  feeding  dive  depths  and  durations,  and  sampling 
zooplankton  to  identify  life  cycle,  concentration,  and  distribution  (MOCNESS  and 
OPC),  they  determined  a  lower  ingestion  rate  of  4.61  x  108  copepods  per  day, 
and  determined  that  NARW  specifically  select  stage  V  Calanus  finmarchicus 
copepodites  (CV)  (Baumgartner  and  Mate  2003;  Durbin  et  al.  2002),  Estimating 
a  wet  weight  of  1 .7  mg  per  CV  (Durbin  et  al.  2002),  a  single  NARW  consumes 
between  7.8  x  108  and  1.6  x  109  mg  of  CV  per  day.  Using  a  mean  estimate  of 
daily  consumption  (1,20  x  109  mg  CV  per  day)  and  mean  values  of  total  PAH  and 
total  PCBs  measured  in  zooplankton  from  the  Bay  of  Fundy,  Gulf  of  Maine,  Cape 
Cod  Bay,  and  Georges  Bank  (Moore  ef  al.  in  preparation;  Weisbrod  ef  al.  2000a), 
we  estimate  that  daily  dosage  of  PAH  and  PCBs  for  a  feeding,  adult  NARW  is 
1.53  x  10'3  and  1,44  x  1C4  mg/kg,  respectively.  The  calculated  daily  dosage  of  a 
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single  PAH,  BaP,  is  9.33  x  10'6  mg/kg.  The  concentrations  of  the  most  potent 
PCBs  (non-ortho-substituted  PCBs  77,  126,  156  and  169)  were  lower  than  the 
limits  of  detection  of  the  method  (Weisbrod  et  al.  2000a),  therefore  total  PCB 
exposure  is  a  slight  underestimate,  but  even  very  small  concentrations  of  these 
PCBs  may  have  significant  effects. 

Three  mono-ortho-substituted  PCB  congeners,  for  which  toxic  equivalency 
factors  (TEFs)  exist,  were  found  in  measurable  amounts  in  NARW  biopsies:  PCB 
105,  1 14,  and  118.  The  TEFs  for  these  PCBs  have  recently  been  adjusted  to 
0.00003,  representing  the  relative  potency  of  these  individual  compounds  to 
produce  biological  (or  toxic)  effects  relative  to  2,3,7,8-TCDD  (van  den  Berg  et  at. 
2006).  Following  the  assumption  of  additivity,  the  toxic  equivalency  (TEQ)  is 
calculated  as  the  sum  of  the  product  of  each  congener’s  concentration  and  its 
TEF.  For  PCBs  105,  114,  and  118,  the  mean  TEQ  for  NARW  biopsies  sampled 
during  summer  and  winter  is  0.01 1,  and  0.032  ng/g  sample  lipid  content, 
respectively  (mean  summer  lipid  content  =  14.8%,  winter  =  3.6%).  Because  the 
TEFs  for  these  PCBs  are  equal,  the  higher  concentration  of  PCB  118  accounts 
for  the  majority  of  the  dioxin-like  activity  in  each  calculation,  as  well  as  the 
seasonal  difference.  Substituting  WHO  2005  TEFs  with  relative  potencies  (REP) 
empirically  derived  from  beluga  whale  in  vitro  AHR  binding  for  PCB  105  and  1 18 
(Jensen  2000)  increases  the  TEQs  above  by  -2. 5-fold.  Although  van  den  Berg 
et  al.  (2006)  highlight  the  wide  variation  in  reported  REP  values  for  mono-ortho- 
substituted  PCBs,  and  attribute  it  to  low-level  contamination  by  more  potent 
dioxin-like  compounds,  Jensen’s  estimate  may  be  a  more  accurate  assessment 
of  cetacean  sensitivities  to  dioxin-like  PCBs. 

Complicating  the  calculation  of  NARW  exposures  to  environmental 
contaminants  is  the  combination  of  this  chronic  low  exposure  and  the  potential 
for  more  acute  (re)exposures  during  fasting  and  migration  during  which  lipid 
mobilization  reduces  the  percent  lipid  approximately  4-fold,  and  concentrations  of 
total  PCBs  increases  --6-fold  in  the  superficial  blubber  (Weisbrod  et  al.  2000a). 
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While  these  are  indeed  very  small  dosages,  chronic  low  exposure  may  still 
produce  biological  effects.  The  predicted  timing  of  ovulation  for  NARW  occurs 
during  the  time  when  persistent  organic  pollutants  are  at  their  highest 
concentration  in  blubber,  increasing  the  potential  for  reproductive  disruption. 

For  comparison,  dietary  exposure  to  BaP  at  low  chronic  dosages  can 
cause  adverse  reproductive  and  developmental  effects  in  mice.  Administration  of 
as  little  as  10  mg/kg  BaP  during  gestation  decreased  gonad  weights  and  reduced 
fertility  and  reproductive  capacity  in  offspring  of  both  sexes  by  impairing  germ  cell 
development  (Kristensen  etal.  1995;  MacKenzie  and  Angevine  1981).  Another 
study  found  that  high  doses  of  BaP  {120  mg/kg/d)  during  gestation  also  caused 
stillbirths,  fetal  resorption,  and  malformations  regardless  of  AHR  allelomorph 
(high  and  low  affinity  AH  receptors)  (Legraverend  et  al.  1984).  Although  in  each 
of  these  studies  exposure  to  BaP  was  much  higher  than  our  estimate  of  the 
exposures  to  NARW,  these  experimental  exposures  were  limited  to  less  than  ten 
days;  longer  chronic  exposures  in  NARW  may  have  similar  effects. 

Because  NARW  population  growth  differs  so  dramatically  from  the 
southern  right  whale  and  other  cetaceans,  it  is  important  to  compare  NARW  with 
other  Cetaceans.  In  comparison  to  NARW,  sympatric  HW  have  more  diverse 
prey  sources,  and  depend  more  heavily  on  schooling  fish  such  as  sand  lance, 
herring,  and  mackerel,  and  occasionally  euphausiids,  in  the  western  North 
Atlantic  (Main  et  al.  1982;  Payne  ef  al.  1990).  These  planktivorous  fishes  have 
different  life  histories,  potentially  exposing  HW  to  different  sources  of  persistent 
environmental  contaminants.  Furthermore,  while  it  is  well  established  that  bony 
fishes  possess  mechanisms  of  detoxication  (Andersson  and  Fortin  1992;  Hahn  et 
al.  1998)  that  could  effectively  reduce  their  PAH  accumulation,  copepods 
(especially  in  diapause)  may  not.  Thus,  HW  may  have  both  more  diversity  of 
DLC  and  less  exposure  to  PAH  in  their  diet  than  NARW.  In  contrast,  BW  in  the 
Gulf  of  St.  Lawrence  suffer  from  heavily  contaminated  food  sources,  and 
themselves  have  high  DLC  burdens  (Martineau  et  al.  1987).  This  BW  population 
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has  a  high  rate  of  disease,  including  a  high  incidence  of  tumors,  reduced  immune 
responses,  and  a  low  reproductive  rate  (Bailey  and  Singer  1995;  DeGuise  etal. 
1995;  Martineau  et  at.  1988;  Martineau  et  at.  1995).  Studies  of  BW  may 
represent  effects  of  extremely  high  and  chronic  exposures,  whereas  NARW  have 
low  daily  exposures  punctuated  by  periods  of  acute  (re)exposure  during  fasting 
or  reproductively  sensitive  times. 

4.  SPECIFIC  OBJECTIVES 

Given  the  obstacles  to  studying  marine  mammal  toxicology,  how  can  we 
investigate  the  specific  factors  and  mechanisms  involved  in  biological  responses 
to  POPs?  Can  molecular  techniques  typically  employed  for  laboratory  systems 
be  applied  in  a  minimally  invasive  manner  to  infer  the  responses  of  marine 
mammals?  Have  mysticetes  evolved  different  mechanisms  in  response  to  dioxin 
like  compounds  that  make  them  more  or  less  susceptible  to  the  effects  of  these 
chemicals?  How  do  marine  mammals  compare  to  other  mammals  in  their 
sensitivity  to  POPs? 

We  began  this  project  with  a  specific  interest  in  the  reproductive  health  of 
NARW.  We  hypothesized  that  NARW  may  be  particularly  sensitive  to  the  effects 
of  POPs.  However,  the  highly  endangered  NARW  is  difficult  to  study,  and  fresh 
tissue  sources  are  rarely  available;  the  exception  is  that  integument  biopsies  are 
routinely  collected  for  population  genetic  studies.  We  sought  to  determine 
whether  RNA  extracted  from  NARW  integument  biopsies  could  be  used  to  clone 
genes  of  toxicological  and  physiological  interest.  We  also  compare  biopsy- 
derived  RNA  with  RNA  from  tissues  collected  from  stranded  cetaceans,  the  most 
frequently  used  and  most  widely  available  cetacean  tissue.  These  results  are 
described  in  Chapter  2. 

The  successful  amplification  of  a  fragment  of  AHR  cDNA  lead  to  further 
characterization  of  the  NARW  in  comparison  with  another  new  mysticete  AHR 
from  humpback  whale,  and  the  previously  described  beluga  whale  (Jensen  and 
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Hahn  2001),  mouse  {Burbach  et  al.  1992),  and  human  (Dolwick  et  al.  1993a) 
AHRs  (Chapter  3).  In  vitro  biochemical  and  cell-culture  assays  were  employed  to 
describe  the  binding  and  transcriptional  activities  of  these  cetacean  receptors. 

We  also  sought  to  place  these  and  other  marine  mammal  AHRs  in  context 
with  mammalian  species  more  frequently  used  in  laboratory  experiments.  While 
the  application  of  the  comparative  approach  in  an  in  vivo  system  is  cumbersome, 
the  extension  of  in  vitro  binding  assays,  combined  with  sequence  analysis  and 
three-dimensional  homology  modeling,  comprise  a  more  tractable  series  of 
analyses  for  this  purpose  (Chapter  4). 
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Figure  1.  The  intracellular  mechanism  of  action  for  the  aryl  hydrocarbon 
receptor. 

Once  bound  with  ligand,  the  cytosolic  AHR  protein  enters  the  nucleus  and 
exchanges  the  Ara9  and  Hsp90  cofactors  for  ARNT.  The  ligand-AHR-ARNT 
complex  locates  xenobiotic  recognition  elements  (XRE)  in  the  promoters  of 
responsive  genes,  and  recruits  the  transcription  initiation  complex  (TIC). 
Transcription  and  translation  of  target  genes  including  cypla  and  cyp2b  follows. 
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Figure  2,  Structural  and  functional  domains  of  the  AHR. 

Aryl  hydrocarbon  receptor  cDNA  (top;  intron/exon  boundaries  =  black  triangles) 
and  protein  structure  characteristics  with  PAS-family  and  AHR-specific  structural 
motifs  indicated,  with  approximate  locations  of  functional  domains  (bottom) 
based  on  the  murine  Ahr  (Burbach  et  al.  1992;  Dolwick  et  ai.  1993;  Fukunaga  et 
al.  1995;  Mimura  et  al.  1994).  Boxes  represent  protein  coding  regions; 
untranslated  regions  of  the  cDNA  are  shown  as  lines.  LBD,  ligand  binding 
domain;  HSP90,  heat-shock  protein  90  dimerization  region;  ARNT,  ARNT 
dimerization  binding  domain;  DBD,  DNA  binding  domain;  TAD,  transcriptional 
activation  domain. 
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Chapter  2. 

Genes  of  toxicological  response  cloned  from  biopsies  of  North  Atlantic 
right  whales  (Eubalaena  glacialis ):  comparison  with  stranded  cetacean 
samples 


35 


ABSTRACT 


The  North  Atlantic  right  whale  (NARW)  is  the  most  endangered  large  whale 
species  with  a  total  population  size  estimated  at  around  350  individuals.  Vessel 
collisions  and  fishing  gear  entanglements  are  the  leading  causes  of  mortality,  but 
variable  reproductive  success  contributes  to  the  species’  failure  to  reach  pre¬ 
whaling  numbers.  Physiological  or  toxicological  stressors  may  play  causative 
roles  in  the  observed  reproductive  variability,  but  are  difficult  to  measure  in  large, 
wild  cetaceans.  Polycyclic  aromatic  hydrocarbons  induce  metabolic  enzymes 
(cytochrome  P450s;  CYP1A1)  through  an  aryl  hydrocarbon  receptor  (AHR) 
mediated  pathway.  Minimally  invasive  integument  biopsies  from  live  cetaceans 
have  given  researchers  access  to  tissue  samples  from  many  individuals.  Here, 
we  isolated  RNA  from  these  samples  to  identify  and  clone  biomarkers  of  toxicant 
response  and  nutritional  status.  The  quality  of  biopsy-derived  RNA  was  better 
than  that  isolated  from  necropsied  skin  and  blubber  samples  of  stranded 
Delphinids  (white-beaked  dolphin,  Atlantic  white-sided  dolphin,  pilot  whale, 
common  dolphin),  and  a  humpback  whale  (HW).  A  fragment  of  an  AHR  cDNA 
was  amplified  by  RT-PCR,  cloned,  and  sequenced  from  NARW,  HW,  and  the 
Delphinids.  A  CYP1 A1  fragment  was  similarly  identified  from  NARW  RNA. 
Estrogen  receptor-alpha  and  aromatase  (CYP19)  were  also  investigated  using 
the  same  techniques.  RNA  from  stranded  cetaceans  showed  a  pattern  of 
degradation  related  to  time  of  sampling  post-mortem;  however,  AHR  and  leptin 
cDNAs  were  amplified  and  cloned  from  some  of  these  tissues.  Although  biopsies 
provide  the  best  source  of  RNA,  necropsies  conducted  on  cetaceans  up  to  24  hr 
post-mortem  are  still  of  value.  Here  we  show  that  RNA  from  biopsies  can  be 
used  to  identify  biomarkers  of  physiological  condition  and  toxicological  response. 
cDNAs  derived  from  this  source  can  be  cloned  for  use  in  gene  expression  and  in 
vitro  functional  studies  to  better  understand  the  molecular  structure,  function,  and 
regulation  of  proteins  involved  in  reproductive  success. 
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1.  INTRODUCTION 


Right  whales,  Eubalaena  glacialis,  were  once  so  abundant  in  the  North 
Atlantic  Ocean  that  they  were  considered  the  preferred  commercial  whaling 
target  Because  they  are  slow  moving  and  float  when  dead,  right  whales  were 
dubbed  the  “right  whale”  to  hunt.  Centuries  of  hunting,  coupled  with  life  history 
characteristics  such  as  slow  reproductive  rates,  contributed  to  the  rapid  decline 
and  near  extinction  of  the  North  Atlantic  right  whale  (NARW).  Unlike  other 
species  (e.g.  grey  whale),  the  NARW  population  has  not  recovered  since  the 
cessation  of  whaling.  Currently  the  NARW  population  is  estimated  around  350 
individuals,  and  while  some  molecular  evidence  suggests  the  population  is 
slightly  larger  (Frasier  2005;  Frasier  ef  at.  2007),  modeling  suggests  that 
extinction  could  be  imminent  within  the  next  150  to  200  years  (Fujiwara  and 
Caswell  2001).  Ship  strikes  and  entanglement  in  fishing  gear  are  now  the 
leading  causes  of  right  whale  mortality;  however,  a  recent  decline  in  calving  rate 
suggests  the  disruption  of  other  biological  factors.  Curiously,  twelve  percent  of 
adult  female  NARW  have  never  reproduced  (Brown  et  al.  1994;  Kraus  ef  al. 

2001 ;  Kraus  et  al.  2007).  The  observed  reproductive  variability  might  be 
explained  by  physiological  (e.g.  nutritional  status)  or  toxicological  (e.g.  responses 
to  contaminants)  factors;  however,  measuring  such  factors  in  large  whales  using 
traditional  in  vivo  methods  is  exceptionally  difficult. 

Numerous  logistical,  ethical,  and  legal  issues  preclude  the  use  of  standard 
experimental  techniques  in  the  study  of  marine  mammal  biology.  Whaling 
collections  have  historically  provided  much  of  the  information  about  cetacean 
physiology,  anatomy,  and  population  status.  Necropsies  of  stranded  specimens 
are  also  useful  in  studying  anatomy,  physiology,  and  life  history  of  cetaceans; 
however,  interpretation  of  these  samples  is  sometimes  confounded  by  poor 
health  condition  of  stranded  individuals  and  decomposition  of  carcasses.  The 
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use  of  minimally  invasive  techniques  for  studying  physiological  parameters  in 
living  whales  has  increased  in  recent  years. 

The  development  of  minimally  invasive  remote  biopsy  has  permitted  the 
sampling  of  integument  from  live,  healthy  cetaceans  (Aguilar  and  Nadal  1984; 
Lambertson  1987).  Molecular  techniques  were  initially  applied  to  karyotype, 
identify,  sex,  and  assess  paternity  of  individual  whales  {Amos  and  Hoelzel  1990; 
Arnason  etal.  1985;  Gowans  et  al,  2000;  Hoelzel  et  al.  1991;  Palsboll  et  al. 

1992).  More  recently  microsatellite  loci  have  been  employed  to  examine  patterns 
of  reproductive  success  and  to  challenge  estimates  of  population  size  in  NARW 
(Frasier  2005;  Frasier  et  al.  2007).  Biochemical  analyses  of  lipid  content, 
contaminant  load,  enzyme  activity  and  localization,  and  stable  isotope  ratios 
have  been  conducted  on  many  cetacean  species  through  the  use  of  biopsy 
collection  (Aguilar  and  Borrell  1991;  Angell  et  al.  2004;  Fossi  et  al.  1992;  Fossi  et 
al.  2000;  Gauthier  et  al.  1996;  Hooker  et  al.  2001;  Ross  etal.  2000;  Todd  et  al. 
1997;  Weisbrod  et  al.  2000a).  Integumentary  biopsies  have  also  been  used  for 
tissue  culture  (Godard  et  al.  2004)  and  cell  culture  of  odontocete  fibroblasts 
(Godard  etal.  2004;  Marsili  etal.  2000). 

Nucleic  acids  in  general,  and  ribonucleic  acid  (RNA)  in  particular,  are 
sensitive  to  degradation  by  endogenous  and  environmental  nucleases. 

Therefore,  sample  preservation  is  critical  to  obtaining  high  quality  material  for 
genetic  analyses,  and  collection  of  tissue  for  this  purpose  during  cetacean 
necropsies  is  often  neglected  except  under  the  most  optimal  conditions  (<6  hours 
post-mortem,  cool  temperatures,  clean  and  convenient  sampling  location).  When 
we  began  this  work,  cetacean  biopsies  had  been  used  for  isolation  and  analysis 
of  DNA,  but  there  were  no  reports  of  RNA  isolated  from  biopsies.  The 
application  of  existing  gene  cloning  techniques  to  RNA  collected  from  biopsies  is 
emerging  as  a  new  approach  to  understand  cetacean  biology. 

Previous  studies  have  addressed  toxicological  and  reproductive  issues  in 
marine  mammals  at  the  level  of  protein  expression  and  activity.  Variations  in  aryl 
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hydrocarbon  receptor  (AHR)  structure  and  function  contribute  to  species-specific 
susceptibility  to  polycyclic  aromatic  hydrocarbons  (PAH)  and  planar  halogenated 
aromatic  hydrocarbons  (PHAH).  AHR  molecular  structure,  binding  capacity,  and 
tissue-specific  expression  has  been  characterized  for  beluga  whale  (Jensen 
2000;  Jensen  and  Hahn  2000;  Jensen  and  Hahn  2001),  harbor  and  Baikal  seals 
(Kim  and  Hahn  2002;  Kim  ef  al.  2002),  and  in  several  other  marine  mammal 
species  (Lapseritis  and  Hahn  2001;  Lapseritis  and  Hahn  in  preparation).  Some 
AHR  ligands  induce  cytochrome  P450  1A  (CYP1A)  enzymes,  a  biomarker  of 
PAH/PHAH  exposure,  and  CYP1A  protein  expression  and  activity  have  been 
shown  in  many  cetacean  species  (Angell  etal.  2004;  Fossi  ef  al.  2003;  Goksoyr 
et  al.  1988;  Goksoyr  et  al.  1986;  White  etal.  1994;  White  etal.  2000;  Wilson  et 
al.  2005a).  Another  cytochrome  P450  (CYP19/aromatase)  is  among  the 
enzymes  responsible  for  conversion  of  androgens  to  estrogens,  a  critical  aspect 
of  reproductive  function,  and  has  been  studied  in  white-sided  dolphin  (Wilson  et 
al.  2005b).  For  cetaceans  and  other  mammalian  species,  successful 
reproduction  depends  on  a  number  of  biological  factors,  including  sufficient 
energy  (lipid)  storage  (Lockyer  1986).  Mammalian  fat  deposition  and  energy 
utilization  are  strongly  influenced  by  leptin,  a  protein  that  also  regulates 
numerous  reproductive  functions  at  least  in  part  through  CYP19  and  estrogen 
receptor  alpha  (ER«)  expression  (Catalano  ef  al.  2003;  Catalano  ef  al.  2004; 
Chehab  2000). 

In  order  to  address  questions  concerning  the  possible  role  of  persistent 
organic  pollutants  in  the  decline  of  the  NARW,  we  sought  to  test  the  feasibility  of 
using  integument  biopsies  as  a  source  of  RNA  to  clone  genes  for  contaminant 
susceptibility  and  other  biomarkers  related  to  physiological  condition.  We 
extracted  mRNA  (polyadenylated  messenger  RNA)  from  integument  biopsies  of 
NARW,  and  amplified  complementary  DNA  (cDNA)  sequences  for  several 
expressed  genes:  AHR,  CYP1A1,  CYP19,  ERa,  leptin,  and  actin.  We  also 
compared  the  quality  of  RNA  from  biopsies  to  that  obtained  from  samples 
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collected  during  necropsies  of  stranded  delphinid  cetaceans.  We  show  that  gene 
expression  can  be  detected  at  the  level  of  mRNA  in  cetaceans  using  minimally 
invasive  techniques,  and  demonstrate  a  technique  to  better  understand 
physiological  processes  active  in  whales,  including  toxicological  responses. 

2.  METHODS 

2. 1  Tissue  collection 

Integumentary  biopsies  were  collected  using  retrievable  biopsy  darts  fired 
from  a  crossbow.  Two  biopsies  consisting  primarily  of  epidermal  tissue  were 
collected  from  live  NARW  in  the  Bay  of  Fundy  (field  identification  "K"  &  “R”,  both 
collected  4/6/1999)  using  standard  techniques  (Brown  ef  al.  1991;  Lambertson 
1987)  in  conjunction  with  other  sampling  efforts.  The  samples  were  stored  in 
DMSO  and  immediately  placed  in  liquid  nitrogen  dry  shippers.  Samples  were 
stored  in  liquid  nitrogen  for  34  months  until  processing. 

Post-mortem  surgical  biopsies  were  obtained  from  Delphinid  cetaceans 
found  dead,  or  stranded  live  and  euthanized,  on  the  shores  of  Cape  Cod,  MA, 
USA,  with  the  cooperation  of  the  Cape  Cod  Stranding  Network  (CCSN), 
Integument  including  epidermal,  dermal,  and  subdermal  (skin  and  blubber) 
tissues  were  obtained  from  individual  Lagenorhychus  albirostris  (white  beaked 
dolphin;  WBD),  L.  acutus  (Atlantic  white-sided  dolphin;  AWSD),  Globicephalus 
melas  (pilot  whale;  PW),  and  Delphinus  delphis  (common  dolphin;  CD)  during 
standard  necropsies,  replicating  the  type  of  sample  obtained  from  biopsy  darts. 
Integument,  mammary  gland,  and  ovary  samples  were  collected  shortly  after  the 
death  of  a  juvenile  humpback  whale  (HW;  Megaptera  novaeangliae)  that 
stranded  live  on  a  Cape  Cod  beach.  CCSN  identification  numbers  link  collected 
tissues  with  individual  records  including  necropsy  reports  (Table  1).  Tissues 
were  flash  frozen  and  stored  in  liquid  nitrogen  for  one  to  eleven  months  until 
processed. 


40 


2.2  RNA  extraction 

Total  RNA  was  extracted  using  RNA  STAT60™  (Tel-Test,  Inc.),  with  care 
taken  to  eliminate  ribonucleases  and  to  prevent  cross-species  contamination  by 
diligent  workspace  and  equipment  cleaning,  and  extraction  on  separate  days. 
Briefly,  skin  and  blubber  biopsies  (0.3  - 1 .0  g)  were  pulverized  by  hand  in  a 
mortar  and  pestle  with  continuous  addition  of  liquid  nitrogen.  Tissue  powders 
were  then  homogenized  in  giass  tubes  with  STAT60™,  extracted  with 
chloroform,  and  precipitated  with  isopropanpol/ethanol  according  to  the 
manufacturer’s  protocol.  Total  RNA  quality  was  assessed  by  denaturing  gel 
electrophoresis,  and  concentrations  were  determined  by  spectrophotometry. 
mRNA  was  purified  from  North  Atlantic  right  whale  total  RNA  samples  using  the 
Nucleotrap  mRNA  Purification  Kit  (Clontech)  or  Micro  Poly  Purist  (Ambion). 

2.3  RT-PCR 

Complementary  DNA  was  generated  with  OmniScript  reverse 
transcriptase  (OmniScript  RT  kit,  Qiagen)  using  random  hexamer  primers  and  1 
iig  RNA  (mRNA  for  NARW,  total  RNA  for  all  other  species).  Amplification  was  by 
PCR  with  degenerate,  gene-specific  primers  (Table  2)  and  AmpliTaq  Gold  DNA 
polymerase  (50  pL  reactions:  IX  PCR  Gold  Buffer,  2  mM  MgCh,  200  tiM  dNTPs, 
0. 2-2.0  uM  each  gene  specific  primer,  1.25  U  AmpliTaq  Gold;  Applied 
Biosystems).  Thermocycling  conditions  for  the  PCR  step  were  generally:  95°C 
10  min;  35  cycles:  95°C  15  sec,  50°C  30  sec,  72°C  60  sec;  72°C  7  min;  rapid 
cooling  to  4°C.  PCR  products  were  visualized  by  electrophoresis  on  a  1% 
agarose  gel  with  ethidium  bromide,  then  gel  extracted  for  cloning  (Gene  Clean  II, 
Biol 01 ,  Inc.). 

Primers  and  PCR  conditions  used  to  successfully  amplify  AWSD  CYP19 
in  a  previous  study  were  kindly  supplied  by  J.Y.  Wilson  (Wilson  et  al.  2005b). 

Two  pairs  of  primers  were  used  in  separate  PCR  reactions  (CYP19  2-23f/CYP19 


41 


810-830r  and  CYP19  672-691f/CYP1 91 353-1 373r;  0.4  uM  each)  in  an  initial  PCR 
attempt  utilizing  AmpliTaq  Gold  polymerase,  and  reaction  conditions  as  follows: 
95°C  7  min;  35  cycles:  95°C  30  sec,  50°C  30  sec,  72°C  90  sec;  72°C  10  min; 
rapid  cooling  to  4°C.  Following  unsuccessful  amplification,  the  reaction  was 
repeated  using  Advantage  2  polymerase  (50  pL  reaction:  IX  PCR  Buffer,  IX 
dNTP  Mix,  IX  Advantage  2  polymerase  mix;  Clontech  Laboratories,  Inc.),  1.0  pM 
each  primer,  and  reaction  conditions  as  follows:  95°C  1  min;  30  cycles:  95°C  5 
sec,  60°C  1  min;  72°C  7  min;  rapid  cooling  to  4°C,  Delphinid  cDNA  was  used  in 
the  first  PCR;  NARW  total  RNA  was  used  in  both  sets  of  reactions,  and  human 
placental  cDNA  (Clontech)  was  used  as  a  positive  control  in  the  Advantage  2 
reactions. 

2.4  Cloning  and  Sequencing 

Individual  gene  products  resulting  from  RT-PCR  in  each  species  were 
handled  separately  and  with  the  same  general  methods.  Purified  PCR  products 
were  cloned  into  the  pGEM-T  Easy©  plasmid  vector  (Promega),  and  propagated 
in  E.  coli  JM1 09  cells.  Plasmid  DNA  was  isolated  from  6-12  colonies  per  gene 
using  the  boiling-lysis  method  (Sambrook  et  al.  1989).  Inserts  were 
independently  verified  by  restriction  endonuclease  digestion  (EcoRI;  Promega), 
and  with  internal  PCR  using  nested  primers  for  AHR  (D143F/D328R;  data  not 
shown).  Two  to  four  positive  clones  were  sequenced  bi-directionally  to  obtain 
consensus  sequences.  Delphinid  AHR  and  actin  were  sequenced  using  Ml  3 
primers  and  Big  Dye  Terminator  Ready  Mix©  with  standard  parameters,  purified 
using  G-50  Sephadex  with  Centri-Sep  columns  or  isopropanol  precipitation,  and 
visualized  on  an  ABI  377  DNA  sequencer  (Applied  Biosystems);  other 
sequencing  reactions  were  performed  by  the  University  of  Maine  DNA 
Sequencing  Facility  (Orono,  ME).  Sequences  were  manually  edited  using 
Sequencher  4.2.2  (Gene  Codes  Corp.)  and  MacVector  8.1.1  (Accelrys,  Inc.) 
software,  and  examined  by  NCBI/GenBank  tBLASTx  for  gene  identification 
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{Altschul  et  al.  1997).  Predicted  amino  acid  sequences  were  determined  for 
each  sequence,  and  aligned  with  known  cDNA  sequences  from  other  species 
using  Clustalw  in  MacVector  8.1 .1 . 


3.  RESULTS 

3. 1  RNA  Extraction 

We  first  sought  to  establish  that  RNA  could  be  obtained  from  NARW  integument 
biopsies,  and  to  compare  its  quality  with  that  collected  from  necropsy  samples. 
Using  standard  RNA  extraction  methods,  we  obtained  good  total  and  mRNA 
yields  from  small  tissue  samples  (0.3  -  1 .0  g)  (Table  1).  RNA  quality  was 
determined  by  electrophoresing  equal  volumes  of  extracted  total  RNA  on  a 
formaldehyde  gel  (Figure  1 ).  The  total  RNA  extracted  from  biopsied  NARW 
tissues  was  of  higher  quality  than  that  extracted  from  necropsy  samples,  as 
shown  by  the  increasing  clarity  and  brightness  of  the  18S  and  28S  ribosomal 
RNA  bands.  A  trend  in  RNA  degradation  appears  related  to  time  of  sampling 
post-mortem,  but  may  also  reflect  individual  or  species-specific  differences. 

Spectrophotometry  was  used  to  evaluate  purity  and  determine  quantities 
extracted  from  each  sample.  The  ratio  of  light  absorbance  at  260  versus  280  nm 
(260/280)  provides  a  measure  of  RNA  purity,  with  an  optimum  ratio  of  2.0;  in 
practice,  260/280  ratios  of  1 .8  -  2.0  are  common,  indicating  some  protein 
contamination  (Sambrook  and  Russell  2001).  NARW  biopsy-derived  RNA 
samples  produced  an  average  ratio  of  1 .79,  compared  with  the  average 
necropsy-derived  RNA  260/280  ratio  of  1 .66.  The  greatest  yield  of  total  RNA 
was  extracted  from  one  NARW  biopsy  (1087  ug/g  tissue).  This  is  consistent  with 
expected  yields  from  fresh  tissues  (1000-1500  ug/g  skeletal  muscle  or  brain, 
1000-4000  pg/g  placenta;  RNA  Stat-60™,  Tel-Test,  Inc.).  Interestingly,  the  yield 
of  total  RNA  from  stranded  cetaceans  that  were  sampled  less  than  approximately 
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15  hours  post-mortem  was  similar  to  that  acquired  from  live  biopsy  samples 
(Table  1).  mRNA  was  only  extracted  from  HW  and  NARW  integument  samples, 
providing  yields  consistent  with  those  regularly  obtained  in  this  laboratory 
(-2.5%).  NARW  biopsies  yielded  better  RNA  than  that  from  necropsied 
Delphintds,  both  in  terms  of  quality  and  yield. 


3.2  Gene  Specific  PCR  Products  &  Sequence  Analysis 

Our  next  objective  was  to  determine  whether  RNA  from  NARW 
integument  biopsies  could  be  used  to  clone  cDNAs  encoding  proteins  related  to 
contaminant  effects  and  physiological  condition.  Actin  was  chosen  as  a  positive 
control  given  its  ubiquitous  high  expression  in  mammalian  tissues,  and  was 
successfully  amplified  from  all  samples  using  general  PCR  conditions  and 
degenerate  primers  (Figure  2).  Actin  was  sequenced  from  AWSD,  HW,  and 
NARW  to  verify  specific  PCR  amplification  (GenBank  Accession  #EF091836  - 
EF091838,  respectively). 

3.2. 1  Aryl  Hydrocarbon  Receptor  (AHR) 

The  AHR  is  a  ligand-activated  transcription  factor  that  regulates  the 
expression  of  a  number  of  genes  in  response  to  binding  PAH  and  PHAH,  such  as 
benzo[a]pyrene,  non-ortho-substituted  polychlorinated  biphenyls  (PCBs),  and 
dioxin-like  compounds  (Chen  and  Bunce  2004).  With  the  assistance  of 
numerous  cofactors,  ligand-bound  AHR  binds  response  elements  in  specific 
gene  promoters,  activating  transcription  of  CYP1A1  and  other  genes  (Hankinson 
2005;  Schmidt  and  Bradfield  1996;  Whitlock  1999). 

An  initial  AHR  cDNA  fragment  of  approximately  690  bp  was  isolated  and 
sequenced  from  the  necropsied  samples  as  described  in  Methods  using  general 
PCR  conditions  and  degenerate  AHR  primers.  PCR  products  from  white-beaked 
dolphin,  Atlantic  white-sided  dolphin,  and  pilot  whale  were  cloned  and 
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sequenced,  and  alignments  and  BLAST  searches  confirmed  AHR-specific 
amplification  (GenBank  Accession  #EF091 831  -  EF091833,  respectively).  Once 
we  were  convinced  that  our  primers  could  amplify  cetacean  AHR,  a  larger 
fragment  was  amplified  from  all  RNA  samples  (Figure  2).  The  NARW  PCR 
product  was  cloned  and  sequenced  (Figure  3),  and  subsequently  the  entire 
NARW  AHR  cDNA  was  fully  cloned  and  sequenced  (GenBank  Accession 
#EF091829;  Chapter  3).  These  data  show  that  AHR  is  expressed  in  cetacean 
integument,  including  NARW  superficial  biopsies. 

3.2.2  Cytochrome  P450  1A1  (CYP1A1) 

CYP1A1  expression  is  induced  in  an  AHR-dependent  manner,  and  is  a 
biomarker  of  exposure  to  PAH  and  PHAH  (Nelson  et  al.  1996).  The  CYP1A1 
enzyme  acts  to  hydroxylate  the  compounds  responsible  for  its  expression, 
facilitating  excretion  and  detoxification;  however,  this  mechanism  of  detoxication 
can  proceed  via  reactive  intermediates,  and  occasionally  the  metabolites  exhibit 
increased  toxicity  (Whitlock  1999).  Two  fragments  were  amplified  from  NARW 
mRNA  using  primers  designed  for  dolphin  CYP1A1  (Montie  2006),  and  the 
smaller  (-370  bp)  was  confirmed  by  sequencing  to  be  CYP1A1  (Figure  4; 
GenBank  Accession  #EF091834).  Expression  of  CYP1A1  in  NARW  superficial 
integument  biopsies  suggests  that  AHR  protein  is  transcriptionally  active  in  this 
tissue,  and  that  AHR  ligands  were  present. 

3.2.3  Estrogen  Receptor  alpha  (ERcx) 

Estrogen  receptors  are  also  ligand-dependent  transcription  factors, 
regulating  the  expression  of  a  suite  of  downstream  genes  of  reproductive 
function.  Two  forms  exist  in  mammals,  ERa  and  ERp,  with  overlapping  functions 
(Hall  et  al.  2001 ).  Degenerate  primers  for  a  large  portion  of  the  ERa  cDNA  were 
designed  from  alignment  of  published  sequences  (human,  mouse,  rat,  cow,  and 
pig).  Although  ERa  can  be  expressed  in  human  skin  (Kwon  et  al.  2004),  it  is 
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more  abundantly  expressed  in  reproductive  tissues  (Kuiper  et  al.  1997). 
Therefore,  RT-PCR  was  conducted  using  NARW  RNA  as  well  as  RNA  isolated 
from  humpback  whale  integument,  mammary,  and  ovary  tissue.  Amplification  of 
appropriately  sized  product  was  only  detected  from  HW  mammary-derived  RNA 
(Figure  5A),  but  its  sequence  remains  unconfirmed.  Although  a  lack  of  ER« 
expression  in  cetacean  integument  cannot  be  excluded,  more  degenerate 
primers  designed  to  a  more  conserved  region  of  ERs,  or  less  stringent  PCR 
conditions  might  be  more  successful. 

3.2.4  Aromatase/Cytochrome  P450  19  (CYP19) 

Binding  of  estradiol  with  estrogen  receptors  induces  expression  of  CYP19, 
an  enzyme  that  converts  androstenedione  and  testosterone  to  estrogens 
(Kinoshita  and  Chen  2003).  This  pathway  is  not  only  critical  for  hormone 
synthesis  and  regulation,  but  chemical  cross-talk  between  the  ER  and  AHR 
pathways  provides  potential  for  endocrine  disruption  along  these  two  axes  (Baba 
eta!.  2005;  Matthews  and  Gustafsson  2006), 

Two  attempts  were  made  to  amplify  CYP19,  each  using  two  separate  sets 
of  primers  previously  used  to  amplify  aromatase  in  cetaceans.  Human  placental 
cDNA  produced  bright  bands  of  expected  size  for  both  primer  pairs,  whereas 
multiple  bands  were  visible  for  the  NARW  reactions  (Figure  5B).  Two  faint  bands 
(-800  and  -900  bp)  roughly  equivalent  to  the  band  seen  for  the  human  positive 
control  were  recovered  from  one  NARW  Advantage  2  reaction  (CYP19  2- 
23f/CYP19  810-830r);  however,  sequencing  revealed  non-specific  products. 

3.2.5  Leptin 

Lastly,  we  attempted  to  amplify  a  cDNA  for  NARW  leptin.  Leptin  is  the 
product  of  the  ob/LEP  gene  and  a  regulatory  hormone  secreted  by  many  cell 
types  including  adipocytes  and  placental  tissue  (Friedman  and  Halaas  1998; 
Magarinos  ef  al.  2006).  Leptin  plays  important  roles  in  nutrition  and  reproductive 
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physiology,  controlling  fat  deposition,  energy  expenditure,  pubertal  onset,  and 
fertility  (Chehab  2000).  Oligonucleotide  primers  were  designed  in  the  5'  and  3’ 
UTR  of  published  leptin  { ob )  cDNAs  using  an  alignment  of  sequences  from 
human,  pig,  mouse,  rat,  cat,  dog,  harbor  seal,  and  grey  seal  (Ob-FI  &  Ob-RI , 
Table  1 ).  Primers  used  for  amplification  of  seal  leptin  were  also  synthesized  and 
used  under  the  standard  PCR  conditions  above  (S2  &  AS2,  Table  1 ;  Hammond 
ef  al,  2005).  Only  the  latter  primers  yielded  visible  products;  bands  from  WBD 
and  other  delphinids  were  of  expected  size,  while  those  for  NARW  were  larger 
than  expected  (Figure  6A).  Although  cloning  and  sequencing  of  the  NARW 
products  yielded  non-specific  sequence,  the  WBD  band  was  confirmed  as  leptin 
(Figure  6B;  GenBank  Accession  #EF091835). 


4.  DISCUSSION 

4. 1  Application  of  Molecular  Biology  Methods  to  Cetacean  Biopsies 

Establishing  non-destructive  methods  for  understanding  cetacean 
physiology  and  toxicology  is  difficult  given  the  obstacles  of  studying  these  large 
wild  animals,  but  is  particularly  critical  for  endangered  species.  Applying  the 
tools  of  molecular  biology  toward  biopsy  samples  is  an  approach  that  can 
circumvent  some  of  these  obstacles  by  enabling  the  use  of  in  vitro  experiments 
ranging  from  gene  expression  to  functional  biochemical  studies.  It  is  hoped  that 
biomarkers  in  integument  samples  could  provide  information  on  individual  health 
and  reproductive  parameters  or  activity  of  other  internal  organs;  this  goal  may  be 
achieved  through  a  combination  of  complementary  approaches.  Here  we  show 
that  RNA  from  integument  biopsies  can  be  used  to  identify  biomarkers  of 
physiological  condition  and  toxicological  response.  cDNAs  derived  from  this 
source  can  be  cloned  for  use  in  gene  expression  and  in  vitro  assays  (see 
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Chapter  3)  to  better  understand  the  molecular  structure,  function,  and  regulation 
of  proteins  involved  in  reproductive  success. 

4.2  Cetacean  Integument  Characteristics 

Although  many  general  characteristics  of  mammalian  integument  are 
observed  in  cetaceans,  several  fundamental  differences  are  worth  noting  (Berta 
and  Sumich  1999).  Sweat  and  sebaceous  glands  are  not  present  in  cetacean 
integument,  and  hair  is  virtually  absent;  sparse  vibrissae  are  localized  in  the 
anterior  sensory  regions  and  mouth  in  whales,  and  are  minimally  present  in 
dolphins.  The  functional  replacement  for  thermoinsulatory  hair  in  the  aquatic 
environment  is  an  expanded,  specialized  hypodermis  (blubber).  Cetacean 
epidermis  is  also  relatively  thick  compared  with  that  of  other  mammals,  yet 
characteristically  high  sloughing  rates  reduce  the  depth  of  the  topmost  layer  of 
dead  cells  (stratum  corneum).  Epidermal  replacement  has  been  calculated  at  9 
times  that  of  humans  (Bergstresser  and  Taylor  1977),  with  surface  renewal 
occurring  about  every  2  hours  (Geraci  ef  al.  1986).  The  interface  between 
cetacean  epidermis  and  dermis  is  also  unusual  in  the  elaborate  development  of 
rete  ridges  and  dermal  papillae,  which  may  have  evolved  in  response  to  the  large 
hydrodynamic  and  mechanical  stresses  placed  on  cetacean  skin,  but  also  serve 
to  reduce  the  distance  between  the  mitotically  active  stratum  germinativum  (the 
deepest  epidermal  layer)  and  the  skin’s  surface.  Thus,  rapidly  dividing  and 
differentiating  cells  are  positioned  close  to  the  skin’s  surface  and  are  easily 
sampled  by  remote  biopsy;  these  cells  may  be  performing  important 
physiological  processes  that  can  be  identified  using  the  techniques  we  describe 
here. 

4.3  Higher  quality  RNA  is  found  in  biopsies  than  from  standings  and  can  be 
used  to  clone  genes  expressed  in  skin  and  blubber. 
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Our  results  show  that  integument  from  biopsies  of  live  NARW  and 
stranded  delphinids  can  be  used  to  collect  RNA  of  sufficient  quality  and  quantity 
to  clone  expressed  genes.  Genes  specific  to  chemical  susceptibility  are 
accessible  using  these  techniques,  and  can  be  developed  as  reagents  for  in  vitro 
toxicological  investigation  (Lapseritis  and  Hahn  in  preparation).  We  have  shown 
that  combining  careful  tissue  sampling  already  underway  with  standard  molecular 
biology  techniques  is  a  useful  way  to  identify  species-specific  cDNA  sequences 
in  cetaceans. 

RNA  extracted  from  biopsies  from  living  whales  was  of  higher  quality  than 
RNA  isolated  from  post-mortem  tissues.  RNA  is  particularly  sensitive  to  tissue 
decomposition,  and  stranded  cetaceans  are  often  exposed  to  extreme 
temperatures  that  destroy  cellular  material,  particularly  labile  nucleic  acids. 
Attempts  were  not  made  to  isolate  RNA  from  post-mortem  right  whales  because 
carcasses  are  typically  located  in  remote  areas  that  preclude  sampling,  or  are 
identified  at  advanced  stages  of  decomposition  (and  often  lacking  skin 
altogether).  Because  biopsies  are  routinely  collected  from  NARW  during  annual 
surveys  and  for  other  study  purposes,  and  very  small  subsamples  can  be  used  to 
support  multiple  uses,  biopsies  represent  a  reliable  source  of  fresh  tissue  for 
nucleic  acid  analysis.  Until  recently,  the  use  of  biopsies  for  RNA  had  not  been 
reported  (Lapseritis  and  Hahn  2004;  Spinsanti  ef  a/.  2006). 

Biopsy  samples  provided  the  best  RNA  for  RT-PCR  and  permitted 
detection  of  several  cDNA  sequences  from  the  NARW.  Ideally,  biopsy  samples 
would  uniformly  sample  integument  perpendicular  to  the  surface  of  the  animal, 
but  in  reality  many  biopsies  consist  of  superficial  tissue  as  the  darts  obliquely 
strike  low-profile  animals  upon  surfacing  or  diving.  Our  biopsy  samples  included 
an  oblique  sample  with  undetectable  hypodermis  by  gross  visual  examination. 
Nonetheless,  we  were  able  to  detect  expressed  genes  in  epidermal  and  dermal 
tissues  indicating  that  living  cells  could  be  recovered  from  the  superficial  layers. 
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It  should  be  noted  that  despite  being  partially  degraded,  RNA  recovered 
from  necropsy  samples  was  still  of  sufficient  quality  to  amplify  partial  and  full 
sequences  of  expressed  genes.  Knowing  the  progression  of  RNA  degradation  is 
important  for  researchers  dependent  on  this  tissue  source,  and  for  individuals 
sampling  cetaceans  at  necropsies.  The  quality  of  RNA  appears  to  decrease  with 
increasing  time  sampled  post-mortem.  Other  factors  may  also  contribute  to  this 
difference:  species-specific  differences;  individual  factors  including  health 
condition  and  life  history;  trauma  associated  with  stranding;  and  environmental 
conditions  at  stranding.  For  species  rarely  observed  alive  or  for  stranding  events 
of  unusual  circumstances,  tissue  samples  collected  from  stranded  specimens 
may  still  be  useful  for  identifying  genes  related  to  physiological  condition,  even 
when  sampling  occurs  as  late  as  18  hours  post-mortem.  Such  samples  are  most 
likely  to  be  useful  for  qualitative  rather  than  quantitative  measurements. 

North  Atlantic  right  whales,  and  mysticetes  in  general,  rarely  strand  alive 
in  New  England  waters;  they  are  often  discovered  floating  at  sea  in  advanced 
states  of  decomposition.  Therefore,  direct  comparison  of  NARW  biopsy 
samples  with  necropsy  samples  was  impractical  in  this  study;  however,  a 
sympatric  mysticete  species  {humpback  whale)  was  found  live-stranded  in 
Chatham,  MA,  and  died  on  the  beach  in  the  evening  following  attempts  at 
release.  Necropsy  commenced  the  morning  after  death,  and  samples  were 
obtained  from  integument,  mammary,  and  ovary  tissues.  Total  RNA  was 
extracted  and  mRNA  purified.  The  quality  of  integument-derived  total  RNA  was 
equivalent  to  that  from  other  necropsied  cetaceans  in  this  study,  but  the  ovarian 
RNA  was  badly  degraded.  Actin  was  amplified  by  PCR  from  all  humpback  whale 
tissues  using  methods  described  above,  and  an  AHR  cDNA  was  subsequently 
fully  cloned  and  sequenced  from  integument-derived  RNA  (Chapter  3).  This  may 
be  the  closest  approximation  to  RNA  degradation  expected  in  a  recently  stranded 
NARW,  and  is  consistent  with  our  measurements  from  stranded  odontocetes. 
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Even  less  invasive  techniques  than  biopsies  have  been  developed  to 
assess  biological  parameters  in  large  whales.  Sloughed  skin  and  fecal  samples 
have  been  used  with  success;  however,  the  utility  of  such  samples  is  limited  to 
DNA  and  hormone  assays  (Amos  et  al.  1992;  Clapham  et  al.  1993;  Hunt  et  al. 
2006;  Parsons  2001;  Parsons  et  al,  1999;  Rolland  etal.  2005;  Valsecchi  et  al, 
1998).  DNA  from  these  sources  often  cannot  be  reliably  assigned  to  individuals, 
is  of  suboptimal  quality  and  quantity,  and  can  easily  become  contaminated  by 
prey  DNA  (Taberlet  and  Waits  1998).  While  we  have  not  assessed  the  presence 
of  RNA  in  such  samples,  we  are  not  optimistic  in  applying  our  approach  to  fecal 
or  sloughed  skin  samples.  The  minimally  invasive  nature  of  remote  live  biopsy  is 
justified  for  gaining  access  to  the  quality  and  quantity  of  RNA  shown  here,  and 
promises  to  increase  our  knowledge  of  cetacean  physiology. 

4.4  Toxicological  and  Physiological  Applications 

While  some  countries  continue  to  support  so-called  scientific  whaling, 
those  efforts  have  not  advanced  the  understanding  of  cetacean  physiology,  in 
part  due  to  poor  experimental  design  and  failure  to  publish  results  (Gales  et  al. 

2005) .  Some  have  called  for  expansion  of  non-invasive  methods  for  acquiring 
information  on  large  whale  physiology,  specifically  endocrine  function  (Hunt  ef  al. 

2006) .  The  application  of  minimally  invasive  sampling  with  qualitative  gene 
expression  techniques  outlined  here  promises  to  illuminate  characteristics 
relevant  to  understanding  toxicological  response,  reproductive  state,  and 
nutritional  status  in  cetaceans.  Understanding  of  cetacean  reproductive 
endocrinology  could  be  improved  using  gene  constructs  developed  from  biopsies 
for  a  number  of  potentially  expressed  genes  such  as  steroid  hormone  receptors 
(currently  under  investigation)  and  their  downstream  gene  targets.  This  report 
identifies  one  marker  of  nutritional  status  (leptin/obese);  however,  other 
biomarkers  of  lipid  metabolism  and  immune  response  may  be  investigated  in  a 
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similar  manner  (e.g.  perioxisome  proliferator-activated  receptors,  interleukins  2, 

6,  &  8  (Chinetti  et  al.  2000)). 

Detection  of  aryl  hydrocarbon  receptor  expression  in  skin  of  living  NARW 
is  important  in  establishing  a  mechanism  of  action  for  some  dioxin-like 
environmental  contaminants  and  PAHs  -  chemicals  known  to  be  present  in  the 
habitat.  Previous  studies  have  shown  that  NARW  carry  body  burdens  of  AHR 
agonists  including  PAHs,  PCBs,  and  DDT  derivatives  (Moore  et  al.  in 
preparation;  Weisbrod  et  al.  2000b).  Human  keratinocytes  express  AHR, 
CYP1A1,  and  other  AHR-regulated  CYPs,  and  are  responsive  to  exposure  to 
AHR  ligands  (Reiners  et  al.  1998;  Swanson  2004),  and  skin  is  a  potential  target 
tissue  for  PAH  and  dioxin-like  compounds  that  activate  CYP  enzymes  in  an  AHR- 
dependent  manner  (Reiners  etal.  1998;  Saarikoski  et  al.  2005;  Tauchi  etal. 
2005).  The  presence  of  AHR  transcripts  in  cetacean  skin  indicates  the  potential 
for  activation  of  the  CYP1A  pathway;  detection  of  CYP1A1  transcripts  supports 
this  hypothesis. 

Previous  studies  have  documented  the  presence  of  immunoreactive 
CYP1A1  protein  in  cetacean  vascular  endothelial  cells  (Godard  et  al.  2004; 
Wilson  et  al.  2005a),  including  in  NARW  integument  biopsies  (Angell  et  al.  2004). 
Although  we  did  not  specifically  investigate  the  presence  of  endothelial  tissue  in 
our  biopsy  samples,  the  presence  of  vasculature  is  likely  to  be  minimal  given  the 
oblique  angle  of  sampling  and  from  gross  examination,  but  histological  analysis 
will  be  necessary  to  confirm  the  degree  of  vascularization.  The  predominant 
route  of  exposure  of  cetaceans  to  AHR  agonists  is  probably  systemic  via  dietary 
consumption  because  measured  PAH  and  PCB  levels  in  water  are  very  low 
(Moore  et  al.  in  preparation);  however,  dermal  exposure  cannot  be  excluded. 
Alternatively,  autocrine  or  paracrine  signaling  may  have  induced  AHR  and 
CYP1A1  expression  within  adipocytes  or  in  neighboring  cells.  Quantitative 
measurement  of  expressed  AHR  and  CYP1A1  transcripts,  paired  with 
histological  examination  by  immunohistochemistry,  could  help  clarify  which  cell 
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types  express  each  protein  and  whether  the  source  of  AHR  agonists  was  internal 
or  external. 

Several  studies  of  estrogen  receptors  and  aromatase  in  human  skin  have 
localized  expression  to  dermal  papillae  associated  with  hair  follicles  and 
sebaceous  glands,  and  shown  greater  expression  of  ERfi  than  of  ERa  (Kwon  et 
at.  2004;  Thornton  et  at.  2006;  Thornton  et  at.  2003).  Given  the  absence  of  these 
features  in  cetacean  skin,  the  probability  of  ER  expression  (and  notably  the  alpha 
form)  is  reduced.  Yet  both  ER  isoforms  are  expressed  in  human  adipocytes 
(Pedersen  et  at.  2001),  and  CYP19  is  abundantly  expressed  in  human  adipose 
stromal  cells  (Price  and  O'Brien  1993),  suggesting  the  real  possibility  of 
identifying  these  proteins  in  cetacean  integument  biopsies  if  sampling  includes 
hypodermis.  Identification  of  these  sequences  would  be  important  for  further 
characterization  of  cetacean  reproductive  endocrinology,  and  holds  potential  for 
more  directly  assessing  issues  of  reproductive  failure  in  NARW.  We  are 
continuing  efforts  to  amplify  both  ERa  and  ERp  cDNAs,  using  more  degenerate 
primers  targeted  to  a  shorter  segment  of  the  cDNA,  and  less  stringent  PCR 
parameters.  Combined  with  new  methods  for  sampling  hormones  from  whale 
feces  (Hunt  et  at.  2006;  Rolland  ef  at.  2005),  characterization  of  steroid  hormone 
receptors  for  in  vitro  studies  of  gene  expression  and  biochemistry  could  prove  to 
be  a  powerful  approach  for  understanding  NARW  reproductive  endocrinology, 
even  on  an  individual  basis. 

The  structure,  function,  and  expression  of  leptin  in  NARW  are  of  particular 
interest  because  of  observed  reproductive  failure  and  poor  body  condition  of 
females  (Angell  2005;  Pettis  et  at.  2004).  Because  of  the  regulatory 
interconnectedness  of  estrogen  receptors,  CYP19,  and  leptin  (Catalano  etal. 
2003;  Catalano  et  at.  2004),  and  their  roles  in  lipid  storage  and  reproduction, 
understanding  how  these  proteins  function  (or  are  disrupted)  in  NARW  is  a 
priority  for  physiological  investigations.  We  successfully  amplified  and  cloned  a 
segment  of  leptin  cDNA  from  white-beaked  dolphin  integument.  This  is  the  first 


53 


indication  of  leptin  expression  in  cetaceans,  although  a  leptin  sequence  has  been 
deduced  from  beluga  genomic  DNA  (Doyon  et  al.  2001).  Right  whale  leptin 
cDNA  might  be  recovered  from  superficial  biopsies  under  slightly  different  PCR 
conditions,  or  with  alternative  primers.  That  we  were  able  to  amplify  a  segment 
of  leptin  cDNA  from  a  fresh  necropsy  sample,  but  not  a  biopsy,  points  to  an 
important  difference  in  the  tissues  collected.  Because  epidermis  is  much  thicker 
in  mysticetes  than  in  odontocetes,  and  the  total  depth  of  sampling  was  equivalent 
(-2  cm  deep),  our  odontocete  samples  contained  a  large  portion  of  hypodermis, 
whereas  the  biopsy  samples  predominantly  contained  epidermis  and  dermis. 

Our  inability  to  isolate  a  leptin  cDNA  from  NARW  may  be  due  to  a  lack  of 
appropriate  tissue  expressing  leptin;  full  depth  biopsies  including  more  adipose 
tissue  may  be  more  effective.  This  may  also  be  true  for  ER«  and  CYP19,  which 
should  be  more  highly  expressed  in  adipose  tissue.  However,  our  ability  to  clone 
leptin  from  a  stranded  cetacean  provides  hope  of  acquiring  leptin  sequence  from 
suitable  NARW  tissue  in  the  future. 

4.5  Summary 

The  isolation  of  good  quality  RNA  from  cetacean  biopsies  is  best  applied 
in  combination  with  other  methods  such  as  cDNA  cloning,  functional 
biochemistry,  immunohistochemistry,  environmental  contaminant  data,  hormone 
measurements,  and  gene  expression  analyses.  The  initial  qualitative  process  of 
gene  discovery  illustrated  here  can  be  followed  by  quantitative  methods,  such  as 
real-time  RT-PCR  (qRT-PCR).  Once  sequences  have  been  identified,  qRT-PCR 
can  be  used  to  determine  expression  of  genes  of  interest.  In  the  time  since  we 
first  reported  using  RNA  from  cetacean  biopsies  (Lapseritis  and  Hahn  2004), 
another  group  has  shown  its  use  in  qPCR  {Spinsanti  et  al.  2006).  The  actin 
sequences  shown  here  can  be  important  for  applying  this  technique  to  NARW, 
HW,  and  AWSD,  and  are  available  to  other  researchers  via  GenBank. 
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Our  results  show  that  the  quality  of  RNA  recovered  from  biopsies  is 
greater  than  from  necropsied  tissues,  and  validates  the  idea  that  gene 
expression  can  be  detected  in  cetaceans  using  minimally  invasive  techniques. 
Expressed  genes  can  be  cloned  from  biopsy  sources,  and  used  to  examine  the 
physiological  status  of  living  cetaceans.  Extension  of  the  methods  described 
here  should  include  qRT-PCR,  paired  with  complementary  techniques  like 
immunohistochemistry  to  determine  the  extent  of  transcript  and  protein 
expression,  and  tissue-  or  cell-specific  expression.  Identification  of  cetacean 
AHR  and  CYP1 A1  sequences  permits  in  vitro  examination  of  toxicological 
processes  otherwise  unachievable  in  living  whales  because  of  ethical,  legal  and 
logistical  obstacles.  These  advances  offer  insight  at  the  molecular  level,  and 
promise  to  increase  our  understanding  of  cetacean  physiology  and  toxicological 
responses. 


New  sequences  from  this  study  are  available  in  GenBank,  Accession 
#EF091829,  EF091831  -  EF091838. 
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Table  1.  Sample  types  and  RNA  yield. 


Species 

Sample  Type 

Sample 

Total  RNA 

Absorbance 

mRNA 

CCSN  Id. 

Mass  (g) 

Yield  (ng/g) 

260/280 

Yield  (pg/g) 

Common  Dolphin 

post-mortem 

i.o 

62 

1.60 

na 

(Delphinus  delphis) 
CCSN-02-01 1-Dd 

(>24  h) 

Humpback  Whale 

post-mortem 

0.4 

469 

1.66 

14 

(Megaptera  novaeangliae) 
CCSN-02-259-Mn 

{16  -  20  h) 

Pilot  Whale 

post-mortem 

0.7 

477 

1,65 

na 

(Globicephalus  metas) 
CCSN-02-004-Gm 

{-18  h) 

Atlantic  White-Sided  Dolphin 

post-mortem 

0.7 

1080 

1.71 

na 

(Lagenorhyncus  acutus) 
CCSN-01*168-La 

(12-18  h) 

White-Beaked  Dolphin 

post-mortem 

1.0 

837 

1.68 

na 

(Lagenorhyncus  albirostris) 
CCSN-01-073-La 

(<  6  h) 

North  Atlantic  Right  Whale 
(Eubalaena  glacial  is) 

live  biopsy 

0.3  K 

1087 

1.81 

27 

0.4  R 

993 

1.77 

23 

K,  Field  sample  "K"  collected  6  April  1999  in  Bay  of  Fundy. 

R,  Field  sample  “R"  collected  from  a  female  6  April  1999  in  Bay  of  Fundy. 
na  -  data  not  available 
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Table  2.  Gene-specific,  degenerate  oligonucleotide  primers  used  for 
polymerase  chain  reactions.  AHR,  aryl  hydrocarbon  receptor;  CYP1A1, 
cytochrome  P450-1a1;  CYP19,  cytochrome  P450/aromatase;  OB,  leptin;  ER- 
a,  estrogen  receptor  alpha. 


Primer  name 

Seauence  (5’-*  3’) 

V 

Actin  Deg  5'  b 

ACA  ACG  GYT  CSG  GYA  TGT  GC 

68°C 

Actin  Deg  3’  5 

GAA  GCA  YTT  GCG  RTG  WAC  RAT 

65°C 

AHR  D33F  c 

GTC  IAA  CCC  ITC  IAA  GMG  ICA  YMG 

64  °C 

AHR  D1t8Fc 

CTS  CTS  CAR  GCI  CTS  AAY  GGI  TT 

65°C 

AHR  D143F  c 

CAR  GAY  TAY  CTI  GGI  TTY  CAI  CA 

62  °C 

AHR  D328R  c 

GCI  CGC  TGR  ATR  AAY  TGR  TAI  CC 

62°C 

AHR  D348R  c 

CAT  ICC  RCT  YTC  ICC  IGT  YTT 

62°C 

CYP1  A1-300Fdd'e 

GGA  YAA  GAG  RCT  GGA  CGA  GAA  TGC 

66°C 

CYP1A1-667Rdde 

GCC  ACT  GGT  TYA  CAA  AGA  CAC  ARC 

65  °C 

CYP19  2-23f d  f 

TGN  TTT  TGG  AAA  TGC  TGA  ACC  C 

67°C 

CYP19  672-692f df 

GCA  AGC  TCT  CCT  YMT  CAA  AC 

ere 

CYP19  810“830rd  l 

TCC  ATR  BDG  TCT  TCC  AGT  TTC 

60°C 

CYP19  1353-1373rd,f 

AAK  CGT  CTC  ARA  AGT  GTR  ACC 

55°C 

OB-F1 d 

GYT  KCN  AGG  YCC  NAG  AAG  CAC  ATC  C 

64°C 

OB-R1  d 

GGN  GKC  TTM  TGC  CTT  TGG  AAG  AGT  RG 

70°C 

OB  S2  d  9 

AGT  CCA  GGA  TGA  CAC  CAA  AAC  C 

66°C 

OB  AS2  d  9 

GYT  CAG  RGC  CAC  CAC  CTC  YGT 

67°C 

FI-ER-a  d 

CAC  ACC  AAA  GCV  TCB  GGV  ATG  G 

65°C 

RI-ER-a  d 

ATG  CTC  CAT  GCC  TTT  GTT  ACT  CAT  GTC  CC 

75°C 

a  Tm  =  melting  temperature;  calculated:  Tm=  81.5  +  (0.41  *  %GC)  *  (675/length); t:  Invitrogen 

Corporation,  Custom  DNA  Oligos.  Carlsbad,  CA,  USA;  c  Integrated  DNA  Technologies,  Inc,, 
Coralville,  IA,  USA; d  Sigma-Genosys,  The  Woodlands,  TX,  USA;  sequences  provided  by; e  E.W. 
Montie  (Woods  Hole  Oceanographic  Institution); 1  J.Y.  Wilson  (McMaster  University); 9  (Hammond 
etal.  2005). 
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NARW  WBD  AWSD  PW  CD 
Live  <6  h  12-18  h  18  h  >24  h 

Figure  1.  Evaluation  of  RNA  quality  from  post-mortem  and  live-biopsied 
cetacean  integument  samples. 

Total  RNA  was  assessed  by  denaturing  gel  electrophoresis.  RNA  recovered  by 
biopsying  live  north  Atlantic  right  whales  (NARW)  was  found  to  be  of  highest 
quality.  A  trend  in  RNA  degradation  appears  related  to  time  of  sampling  post¬ 
mortem,  but  may  also  reflect  individual  or  species-specific  differences,  or 
environmental  conditions  during  a  stranding  event.  Abbreviations:  common 
dolphin  (CD),  pilot  whale  (PW),  Atlantic  white-sided  dolphin  (AWSD),  and  white- 
beaked  dolphin  (WBD). 
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NARW-K  NARW-R  WBD  AWSD  PW  CD 
1  kb  actin  AHR  actin  AHR  actin  AHR  actin  AHR  actin  AHR  actin  AHR  1  kb 


Figure  2.  Actin  and  aryl  hydrocarbon  receptor  are  expressed  in  cetacean 
integument. 

Reverse-transcriptase  PCR  amplified  fragments  of  actin  (1.2  kb;  1.0  liM ,  5’  &  3' 
actin  degenerate  primers)  and  AHR  (-945  bp;  1.0  pM,  D33F/D348R  primers) 
cDNAs  were  visualized  by  agarose  gel  electrophoresis  (1%)  with  1  kb  markers. 
PCR  products  were  isolated  from  north  Atlantic  right  whale  total  RNA  (NARW-K) 
and  messenger  RNA  (NARW-R).  and  from  Delphinid  total  RNA  (WBD  and  CD 
AHR  bands  are  faint,  but  present).  An  initial  AHR  cDNA  fragment  was  isolated 
and  sequenced  from  Delphinid  RNA  samples  (690  bp;  1.0  uM  primers 
D118F/D348R:  WBD.  AWSD.  PW). 
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Figure  3.  Nucleotide  sequence  alignment  of  partial  aryl  hydrocarbon 
receptor  cDNAs  from  four  cetaceans. 

The  consensus  nucleotide  sequence  for  the  cetacean  AHR  is  shown  below  the 
alignment  of  NARW,  WBD,  AWSD,  and  PW  fragments.  Identical  bases  are 
indicated  with  and  shaded.  The  NARW  AHR  has  subsequently  been  fully 
sequenced.  GenBank  Accession  #EF091829. 
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Figure  4.  Partial  cytochrome  P450  1A1  cDNA  isolated  from  NARW  biopsy. 

A.  Two  products  were  amplified  from  NARW  mRNA  using  CYP1A1  primers 
(-370  and  -700  bp;  standard  PCR;  40  cycles,  0.5  uM  each  primer).  The  smaller 
band  is  indicated  by  the  arrow.  Human  placental  cDNA  was  used  as  a  control; 
however,  the  degeneracy  of  one  primer  precluded  binding  and  amplification. 

B.  A  BLAST  search  identified  the  sequence  of  the  smaller  NARW  band  as 
CYP1 A1 ;  the  larger  band  contained  an  intron  (Altschul  et  al.  1997).  NARW 
CYP1 A1  nucleotide  sequence  (GenBank  Accession  #EF091834)  is  shown 
aligned  with  partial  CYP1A1  sequences  from  minke  whale  (MW),  AWSD,  and 
striped  dolphin  (SD);  GenBank  Accession  #A231891 ,  AY641536,  and  AF235141, 
respectively). 
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Figure  5.  Estrogen  receptor  alpha  and  CYP19  PCR  reactions:  mysticete  and 
human  cDNA. 

A.  Humpback  whale  (HW)  mRNA  isolated  from  mammary  tissue  (HW-mam),  but 
not  that  from  integument,  ovary,  or  NARW  mRNA,  produced  a  faintly  visible  PCR 
product  using  primers  designed  to  amplify  a  large  portion  of  ERa  (1.5  kb 
expected  size;  1 .0  pM,  F1-ER-a/R1-ER-a),  The  identity  of  this  fragment  has  not 
yet  been  verified,  and  further  efforts  to  clone  ERa  and  ERB  from  NARW  and  HW 
are  in  progress. 

B.  Two  pairs  of  primers  were  used  in  separate  PCR  reactions  to  amplify  NARW 
and  human  CYP19  cDNA  using  Advantage  2  polymerase  (a:  CYP19  2- 
23f/CYP  19  81 0-830r,  830  bp  expected  size;  b:  CYP19  672-69 If/CYP  19  1353- 
1373r,  700  bp  expected  size;  1.0  pM  each).  Sequencing  of  two  NARW  bands 
using  the  first  pair  of  primers  (-800  and  -900  bp)  yielded  non-specific  results. 

A  and  B  each  represent  single  gels,  but  image  enhancement  was  used  to  better 
illustrate  bands  for  viewing  purposes  in  cases  where  bands  in  some  lanes  were 
much  brighter  than  in  others  (i.e.  CYP19:  human»NARW),  or  when  bands  were 
only  faintly  visible  (HW  ERa). 
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Figure  6.  Leptin  is  expressed  in  white-beaked  dolphin  integument. 

A.  Two  sets  of  primers  were  used  to  amplify  leptin  from  cetacean  cDNAs  (a:  OB- 
F1/OB-R1 ;  b:  OB  S2/OB  AS2;  1 .0  uM  each),  but  only  the  latter  pair  produced 
visible  products.  NARW  PCR  products  were  larger  than  expected  and  were 
found  to  be  non-specific  when  sequenced.  Delphinid  PCR  products  were  all  of 
expected  size  (-360  bp). 

B.  The  WBD  band  was  verified  by  sequencing  and  a  BLAST  search  as  a  partial 
leptin  cDNA  (Altschul  ef  a/.  1997).  The  WBD  nucleotide  sequence  (GenBank 
Accession  #EF091835)  is  shown  aligned  with  partial  leptin  sequences  from  pig, 
harbor  seal  (HS),  and  human  (HU)(GenBank  Accession  #NM_213840, 
AJ618981,  and  NM_000230,  respectively). 


72 


1  H  3  *5 


73 


Chapter  3. 

Molecular  and  Biochemical  Characterization  of  Aryl  Hydrocarbon 
Receptors  in  North  Atlantic  Right  Whale  (Eubalaena  glacialis)  and 
Humpback  Whale  ( Megaptera  novaeangliae) 
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ABSTRACT 

The  aryl  hydrocarbon  receptor  (AHR)  is  a  ligand-activated  transcription  factor 
that  binds  with  some  environmental  contaminants  and  regulates  gene  expression 
and  physiological  responses.  To  investigate  the  hypothesis  that  North  Atlantic 
right  whales  ( Eubaiaena  glacialis,  NARW)  are  highly  susceptible  to  chemical 
pollutants  such  as  polycyclic  aromatic  hydrocarbons  (PAH)  and  PCBs,  RNA 
prepared  from  NARW  integument  biopsies  was  used  to  clone  and  sequence  a 
cDNA  encoding  AHR.  A  predicted  850  amino  acid  protein  (95.8  kD)  was 
identified  as  NARW  AHR.  An  AHR  cDNA  was  also  cloned  from  a  sympatric 
species,  humpback  whale  ( Megaptera  novaeangliae ;  HW)  using  the  same 
techniques.  Pairwise  ClustalW  analysis  of  the  predicted  amino  acid  sequences 
revealed  94-97%  identity  among  sequenced  cetacean  AHRs.  Specific  binding  of 
in  vitro  expressed  NARW  and  HW  AHR  proteins  to  a  prototypical  AHR  ligand 
([3H]TCDD;  2,3,7,8-tetrachlorodibenzo-p-dioxin)  was  shown  by  velocity 
sedimentation  using  sucrose  density  gradients.  Although  the  mean  relative 
binding  affinities  of  mouse  and  human  AHRs  were  statistically  different  (p<0.05; 
ANOVA/Tukey's  post-hoc),  whale  AHRs  did  not  statistically  differ  from  that  of 
either  the  high-affinity  mouse  AHRb_1  (Kd  =  0.36  nM),  or  the  lower  affinity  human 
AHR  (Kd  =  0.97  nM)  in  saturation  binding  assays.  Whale  AHRs  appeared  to 
have  intermediate  affinities  for  [3H]TCDD  (Kd  =  0.64  -  0.85  nM),  but  larger 
binding  capacities  (Bmax)  than  mouse  and  human  AHR  proteins,  which  may  be 
related  to  increased  whale  AHR  protein  expression  in  this  assay.  Whale  AHRs 
could  activate  transcription  of  a  reporter  gene  in  a  mammalian  cell  culture  system 
transiently  transfected  with  AHR  constructs,  independent  of  TCDD  treatment. 
Benzo[a]pyrene  treatment  induced  ligand-dependent  reporter  gene  expression 
for  beluga  and  mouse  AHRs  in  a  transient  transfection  assay,  offering  preliminary 
evidence  for  combined  ligand-  and  species-specific  activity  for  mammalian 
AHRs. 
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1.  INTRODUCTION 

North  Atlantic  right  whales  are  critically  endangered,  with  a  population  size 
considerable  smaller  than  that  of  other  baleen  whales.  Despite  a  70-year  respite 
from  commercial  whaling,  the  species  has  failed  to  recover  and  may  be  declining. 
Recent  estimates  of  population  growth  rate  range  from  -0.02%  to  1%  (Caswell  et 
al.  1999);  Tim  Frasier,  personal  comm.)  yet  do  not  approach  that  of  the  southern 
right  whale  species  that  appears  healthy  (7%;  source).  Collisions  with  ships  and 
entanglement  in  fishing  gear  cause  the  direct  mortality  of  NARWs,  and  have 
proportionally  larger  impacts  on  pregnant  and  lactating  females  and  their 
neonates.  However,  other  causes  of  the  species’  failure  to  recover  may  include 
nutritional  stress  due  to  variations  in  prey  availability,  reproductive  effects  of 
exposure  to  marine  biotoxins  or  anthropogenic  chemical  contaminants,  and 
inbreeding  suppression  brought  on  by  the  extreme  population  reduction  caused 
by  early  20th  century  whaling. 

Concern  has  grown  among  conservationists  and  scientists  over  the  levels 
of  persistent  anthropogenic  contaminants  present  in  the  marine  environment, 
including  pristine  areas  far  from  point  sources  of  industrial  waste  products. 

Some  classes  of  contaminants,  dioxin-like  compounds  (DLC)  such  as  chlorinated 
dibenzo-dioxins,  polychlorinated  biphenyls  (PCBs),  and  polycyclic  aromatic 
hydrocarbons  (PAHs),  can  cause  a  variety  of  non-lethal  pathologies  at  relatively 
low  but  chronic  doses.  Associations  between  hormonal  and  morphological 
disorders  at  all  age  classes  and  many  stages  of  gestation  have  been  observed  in 
many  cetacean  and  pinniped  species  (Reijnders  2003).  However,  it  is  often 
difficult  to  clearly  determine  whether  these  pathologies  in  wildlife  are  the  results 
of  contaminant  exposure  because:  a)  many  species  are  not  well  studied  and  the 
background  incidence  of  these  effects  is  unknown,  b)  species-specific 
differences  in  doses  required  for  toxicity  are  unknown,  and  in  some  cases 
impossible  to  determine  using  standard  methods,  and  c)  interaction  of  the 
multitude  of  environmental  and  natural  chemicals  evades  simple  explanations. 
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Another  approach  is  required  to  assess  the  potential  impacts  of  DLC  on 
endangered  wildlife  species. 

Most  vertebrate  species  possess  a  mechanism  for  detection  and 
metabolism  of  DLC.  The  aryl  hydrocarbon  receptor  is  a  cellular  protein  that  binds 
DLC,  and  through  its  action  as  a  transcription  factor,  activates  the  transcription  of 
genes  encoding  a  suite  of  metabolic  enzymes  and  other  genes.  While  in  some 
cases  this  results  in  the  detoxication  of  the  chemical,  in  others  the  chemical's 
toxicity  is  increased  via  reactive  intermediates.  The  rate  of  clearance  and  the 
dose  and  duration  of  exposure  can  also  affect  the  degree  of  toxicity  experienced 
by  an  organism.  The  genetic  background,  and  in  fact  species-specific 
differences  in  AH  receptor  sequence,  can  dramatically  influence  the  sensitivity  of 
an  organism  to  DLC.  Mice,  rats,  and  birds  show  AHR  sequence-specific 
differences  that  relate  to  sensitivity  (Karchner  et  al.  2006;  Pohjanvirta  et  al.  1998; 
Poland  et  al.  1994).  The  AHR  in  one  cetacean,  the  beluga  whale,  has  been 
studied  using  in  vitro  bioassays  to  determine  its  affinity  for  DLC,  structure-activity 
relationships,  and  tissue-specific  expression  patterns  (Jensen  2000). 

Using  methods  similar  to  those  employed  for  the  investigation  of  beluga 
whale  (BW),  we  assess  the  sensitivity  of  the  NARW  by  identifying  and 
biochemically  characterizing  the  AHR  in  comparison  to  BW  and  a  sympatric 
mysticete,  the  humpback  whale  (HW;  Megaptera  novaeangliae).  We  establish 
the  relative  binding  affinity  of  AHRs  from  these  cetaceans  in  comparison  with 
representative  high-  and  low-affinity  AH  receptors  from  mouse  and  human, 
respectively.  The  transcriptional  activities  of  these  receptors  are  also  examined 
in  a  cell  culture  system  to  approximate  the  in  vivo  condition.  Given  the 
constraints  on  systemically  investigating  contaminant  effects  in  large, 
endangered  cetaceans,  we  show  the  efficacy  and  limitations  of  using  in  vitro 
molecular  biology  and  biochemistry  techniques  to  approximate  the  toxicological 
effects  of  DLC  on  cetaceans. 
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2.  METHODS 


2. 1  Cloning  and  Sequencing  Mysticete  AHRs 

Biopsy  samples  collected  from  two  live  NARW,  and  dissection  of  a  stranded 
HW  provided  integument  tissue  to  determine  AHR  sequence  in  these  species 
(described  in  detail  in  Chapter  2).  Messenger  RNA  (mRNA)  was  isolated  and 
purified  from  NARW  and  HW  samples  (RNA  STAT-60;  MRNA  kit),  and  used  for 
cDNA  synthesis  with  oligo  d{T)  primers  and  Omniscript  RT  polymerase 
(SOURCE).  Another  pool  of  NARW  cDNA  was  synthesized  using  Thermoscript 
RT  (Invitrogen)  and  a  gene  specific  primer  (MM3r). 

An  initial  fragment  of  NARW  AHR  was  amplified  from  Omniscript  cDNA  by 
polymerase  chain  reaction  (PCR)  using  degenerate  primers  (described  in  detail 
in  Chapter  2).  Gene  specific  primers  were  designed  from  this  fragment  for  Rapid 
Amplification  of  cDNA  Ends  (RACE)  (Table  1).  Degenerate  primers  were  also 
designed  based  on  alignments  of  AHR  fragments  from  NARW,  white-beaked 
dolphin  (GenBank  Accession  (GBA)  #  EF091831),  Atlantic  white-sided  dolphin 
(GBA  #  EF091832),  pilot  whale  (GBA  #  EF091833),  mink,  and  polar  bear 
(Chapters  2  and  4). 

SMART  RACE  (Clontech)  was  used  according  to  manufacturer's  instructions 
to  synthesize  two  separate  pools  of  cDNA  from  NARW  mRNA:  5’  RACE  Ready 
and  3'  RACE  Ready  cDNAs.  5’  AHR  RACE  product  was  amplified  from  5’  RACE 
Ready  cDNA  with  a  degenerate  primer  (MM3r,  100  liM)  and  universal  primer  mix 
(UPM;  supplied  with  the  SMART  RACE  kit).  PCR  conditions  were  as  follows:  5 
cycles:  94°C/5  sec,  72°C/3  min;  5  cycles:  94°C/5  sec,  70X/3  min;  25  cycles: 
94°C/5  sec,  68°C/10  sec,  72°C/3  min;  72°C/3  min;  rapid  cooling  to  4°C.  A 
nested  reaction  was  applied  to  this  product  using  a  specific  primer  (RW-1r,  10 
pM)  and  nested  universal  primer  mix  (NUPM;  SMART  RACE).  Nested  PCR 
conditions  were  as  follows:  20  cycles:  94°C/5  sec,  68°C/10  sec,  72°C/3  min. 

Five  additional  primers  were  designed  and  used  to  amplify  the  NARW  AHR  5’ 
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sequence  (not  shown).  Synthesis  of  5’  RACE  Ready  cDNA  was  repeated  using 
a  gene  specific  primer  (MM3r)  in  lieu  of  the  CDS  primer  supplied  by  the 
manufacturer  to  enrich  for  AHR  transcripts,  and  this  cDNA  was  subsequently 
used  with  RW-lr  and  UPM  in  a  RACE  PCR. 

3’  RACE  was  conducted  under  identical  PCR  conditions,  using  3’  RACE 
Ready  cDNA,  first  with  RW-lf  (10  nM)  and  UPM,  followed  by  a  nested  PCR  with 
MM-4F  (10  (.iM)  and  NUPM.  A  primer  specific  to  the  NARW  3’  UTR  was 
designed  for  full-length  amplification  (RW  AHR-1 R). 

NARW  cDNA  synthesized  with  Thermoscript  RT  was  used  with  the  Marathon 
RACE  kit  (Clontech),  following  the  manufacturer’s  instructions.  RACE  PCR  was 
conducted  with  this  cDNA  using  API  and  MM3r  (100  uM)  primers,  and 
Advantage  2  polymerase  with  the  following  protocol:  94°C/30  sec;  5  cycles: 
94°C/5  sec,  72°C/2  min;  5  cycles:  94°C/5  sec,  70°C/2  min;  25  cycles:  94°C/5 
sec,  68°C/2  min;  68  °C/3  min;  rapid  cooling  to  4°C.  A  nested  PCR  reaction  was 
applied  to  this  product  using  AP2  and  MM2r  primers:  94°C/30  sec;  20  cycles: 
94°C/5  sec,  69°C/2  min;  68°C/5  min;  rapid  cooling  to  4°C. 

Full-length  sequence  was  amplified  from  NARW  3'  RACE  Ready  cDNA  with 
a  primer  designed  for  beluga  AHR  (Bel-nat),  RW  AHR-1  R,  and  the  first  PCR 
protocol  above.  Humpback  whale  AHR  was  identified  in  a  single  full-length  PCR 
reaction  with  the  same  gene  specific  primers  using  cDNA  synthesized  with 
Omniscript  RT  polymerase.  Thermocycling  conditions  were  as  follows:  94°C/30 
sec;  5  cycles:  94°C/5  sec,  72°C/4  min;  5  cycles:  94°C/5  sec,  70°C/4  min;  30 
cycles:  94°C/5  sec,  68°C/4  min;  rapid  cooling  to  4°C. 

All  PCR  fragments  and  full-length  AHR  cDNAs  were  gel  purified  (Gene 
Clean),  ligated  into  the  pGEM-T  Easy  plasmid  vector  (Promega),  and  propagated 
in  E.  coli  JM109  cells.  Plasmid  DNA  was  isolated  from  6  colonies  per  fragment 
using  the  boiling-lysis  method  (Sambrook  et  al,  1989).  Inserts  were 
independently  verified  by  restriction  endonuclease  digestion  (EcoRI;  Promega). 
Two  to  six  positive  clones  were  sequenced  bi-directionally  to  obtain  consensus 
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sequences.  RACE  products  were  sequenced  using  Ml  3  primers  and  Big  Dye 
Terminator  Ready  Mix  with  standard  parameters,  purified  using  G-50  Sephadex 
with  Centri-Sep  columns  or  isopropanol  precipitation,  and  visualized  on  an  ABI 
377  DNA  sequencer  (Applied  Biosystems).  Full-length  products  were  sequenced 
by  the  University  of  Maine  DNA  Sequencing  Facility  (Orono,  ME).  Sequences 
were  manually  edited  using  Sequencher  4.2.2  (Gene  Codes  Corp.)  and 
MacVector  8.1 .1  (Accelrys,  Inc.)  software,  and  examined  by  NCBI/GenBank 
tBLASTx  for  gene  identification  (Altschul  et  al.  1997).  Predicted  amino  acid 
sequences  were  determined  for  each  sequence,  and  aligned  with  known  cDNA 
sequences  from  other  species  using  Clustalw  in  MacVector  8.1 .1  (beluga  whale 
AHR  GenBank  Accession  #AF332999;  human  GenBank  Accession  #L1 9872.1; 
mouse  AHR  GenBank  Accession  #M94623). 

2. 2  Expression  Vector  Construction 

Two  sets  of  primers  were  designed  to  correct  a  PCR  error  in  a  NARW 
clone  (RW5;  full-length  sequence)  using  the  Quik  Change  XL  Site-Directed 
Mutagenesis  (SDM)  kit  (Stratagene).  The  first  attempt  failed  to  amplify  plasmid 
DNA  using  10  ng  of  RW5  plasmid  DNA,  and  RW  C2261 A-F  and  RW  G2261T-R 
primers  following  manufacturer’s  instructions  and  the  following  thermocycling 
conditions:  95°C/1  min,  18  cycles:  95°C/50  sec,  60°C/50  sec,  68°C/5.5  min; 
68°C/7  min;  rapid  cooling  to  4°C.  Template  DNA  (RW5)  was  increased  to  50  ng 
in  a  second  SDM  reaction  and  modified  thermocycling  conditions:  95°C/2min,  18 
cycles:  95°C/1min,  60°C/50  sec,  68°C/11  min;  68°C/1 1  min;  rapid  cooling  to4°C. 
A  final  SDM  reaction  used  RW  C2261A2-F  and  RWG2261T2-R  primers  and 
thermocycling  conditions  as  follows:  95°C/2min,  18  cycles:  95°C/1min,  60°C/50 
sec,  68°C/15.5  min;  68°C/15.5  min;  rapid  cooling  to  4°C.  Amplified  plasmids 
were  cloned  and  propagated  in  XL1 0-Gold  cells  and  sequenced  by  the  University 
of  Maine  DNA  Sequencing  Facility. 
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Full-length  PCR  products  of  NARW  AHR  and  HW  AHR  were  initially 
cloned  and  sequenced  in  pGEM-T  Easy  vectors.  The  NARW  AHR  was 
subcloned  into  the  Hindlll  and  Notl  sites  of  the  pcDNA  3.1  expression  vector 
(Invitrogen).  A  chimeric  AHR  was  constructed  from  two  HW  clones  by 
endonuclease  digestion  ( Xmnl/Hindlll  or  Xmnl/Apal)  and  ligation  of  purified  error- 
free  fragments,  then  subcloned  in  pcDNA  at  the  Hindlll  and  Apal  sites.  Correct 
orientation  and  sequence  were  confirmed  by  full-length  sequencing  of  the  final 
constructs.  The  beluga  whale  (BW)  AHR  expression  vector  (pcDNA)  was 
previously  constructed  in  this  laboratory  using  similar  methods  (Jensen  2000; 
Jensen  and  Hahn  2001).  Expression  constructs  for  the  human  AHR  (phuAHR; 
(Dolwick  et  al.  1993a)  and  mouse  AHR  (pmuAHR;  (Dolwick  et  al.  1993b)  were 
cloned  in  pSPORT  vectors  and  graciously  provided  by  Dr.  C.  Bradfield  (McArdle 
Center  for  Cancer  Research,  Madison,  Wl).  Human  (HU)  and  mouse  (MO)  AHR 
were  subcloned  into  pcDNA  vectors  by  Diana  Franks  to  assess  differences  in 
protein  expression  related  to  expression  vector,  but  the  original  constructs 
(pSPORT  vector)  were  used  in  the  ligand  binding  assays  shown  here. 

2.3  In  vitro  Protein  Synthesis 

In  vitro  transcription  and  translation  (IVTT)  of  NARW,  HW,  BW,  HU,  and 
MO  AHR  proteins  was  performed  using  the  TnT  Quick  Coupled  Reticulocyte 
Lysate  Systems  following  the  manufacturer's  protocols  (T7:  NARW,  HW,  BW; 
SP6:  HU,  MO;  Promega).  (35S]-labeled  methionine  was  added  to  IVTT  reactions 
for  visualization  and  quantification  of  protein  synthesis  by  SDS-PAGE  and 
autoradiography.  [35S]-labeled  proteins  were  excised  from  dried  gels  and 
measured  by  liquid  scintillation  with  a  Beckman  LS5000TD  scintillation  counter  to 
quantify  protein  expression  (-5  IVTT  reaction  per  band).  Unlabeled 
methionine  was  used  to  synthesize  IVTT  proteins  for  ligand  binding  assays. 

2.4  Velocity  Sedimentation 
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IVTT  expressed  AHR  proteins  from  NARW,  HW,  BW,  HU,  and  MO  were 
analyzed  by  velocity  sedimentation  on  sucrose  gradients  in  a  vertical  tube  rotor 
(Tsui  and  Okey  1981).  IVTT  reactions  (50  pL)  were  diluted  1:1  with  MEEDMG 
buffer  (25  mM  MOPS  pH  7,5  containing  1  mM  dithiothreitol,  1  mM  EDTA,  5  mM 
EGTA,  0.02%  NaN3,  20  mM  Na2Mo04,  10%  (v:v)  glycerol,  plus  a  mixture  of 
protease  inhibitors  (Hahn  et  at.  1994)).  AHR  proteins  and  unprogrammed  lysate 
(UPL)  were  incubated  with  2  nM  [3H]TCDD  for  18  hours  at  4°C  in  glass  culture 
tubes,  Final  [3H]TCDD  concentration  was  verified  by  sampling  5  pL  from  each 
tube  for  total  counts.  After  incubation,  95  pL  from  each  tube  plus  3  pL 
sedimentation  markers  ([14C]ovalbumin  and  [14C]catalase)  were  applied  to  10- 
30%  sucrose  gradients  prepared  with  MEEDMG  buffer  (Coombs  and  Watts 
1985).  The  gradients  were  centrifuged  at  60,000  rpm  (VTi  65.2  rotor)  for  140 
minutes  at  4°C.  Fractions  (-150  pL  each)  were  counted  using  a  Beckman 
LS5000TD  scintillation  counter.  Specific  binding  was  calculated  as  the  difference 
between  total  binding  and  non-specific  binding  (UPL).  The  sedimentation 
coefficient  for  each  AHR  was  calculated  based  on  the  elution  patterns  of  the 
markers  ([14C]ovalbumin  in  fractions  3-4  (3.6  S)  and  [14C]catalase  in  fractions  14- 
16  (1 1 .3  S))  (Martin  and  Ames  1961 ). 

2.5  Saturation  Binding  Assays 

We  slightly  modified  the  protocol  developed  by  Jensen  (2000,  2001; 
based  on  (Poland  and  Glover  1976))  to  determine  the  specific  binding  of 
[3H]TCDD  for  IVTT-expressed  AHR  proteins.  IVTT  proteins  were  synthesized 
with  unlabeled  methionine  in  two  50  pL  reactions,  combined,  then  diluted  1:9.5  in 
MEEDMG  with  protease  inhibitors.  Diluted  proteins  (105  uL)  were  incubated  in 
glass  culture  tubes  with  one  of  nine  or  ten  concentrations  of  [3H]TCDD  (0  -  10 
nM,  nominal)  in  DMSO  at  4°C.  Each  tube  was  sampled  (5  uL)  for  total  counts  at 
the  beginning  of  the  incubation.  After  18  hours  of  incubation,  28  pL  aliquots  were 
transferred  to  triplicate  1 .5  mL  polypropylene  tubes  with  equal  volumes  of  2.3 
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mg/mL  dextran-coated  charcoal  (DCC).  Tubes  were  vortexed  three  times  for  5 
seconds  at  5  minute  intervals.  DCC  was  sedimented  by  centrifugation  at  7000*g 
for  3  minutes  to  separate  free  ligand  from  protein-bound  [3H]TCDD.  Supernatant 
containing  protein-bound  [3H]TCDD  was  removed  (45  pL)  and  measured  in  the 
Beckman  LS5000TD  scintillation  counter.  NARW,  HW,  BW,  HU,  and  MO  IVTT 
AHRs  plus  UPL  (empty  pcDNA  vector)  were  included  in  each  of  five  assays. 

Radioactivity  from  SBA/DCC  samples  was  detected  by  scintillation 
counting  for  5  minutes  in  each  replicate.  However,  we  questioned  the  accuracy 
of  this  measurement,  especially  for  small  concentrations  of  [3H]TCDD  for  which 
the  detection  error  approached  20%.  We  compared  measurements  of  the  same 
samples  read  for  5,  10,  and  20  minutes,  and  determined  that  longer  counts 
increased  precision  but  did  not  change  accuracy  (data  not  shown). 

Binding  curves  were  plotted  as  free  (unbound)  versus  bound  [3H]TCDD 
(nM).  Disintegrations  per  minute  were  all  converted  to  nM:  [(dpm  - 
background)/specific  activity  of  [3H]TCDD]/volume  per  sample.  Total  binding  was 
determined  by  taking  the  mean  of  each  triplicate  counted  for  each  concentration 
of  [3H]TCDD.  Free  [3H]TCDD  was  calculated  by  subtracting  the  total  bound 
[3H]TCDD  from  the  total  [3H]TCDD  (total  counts)  for  each  tube,  Non-specific 
binding  for  each  concentration  of  [;iH]TCDD  was  determined  from  the  slope  of  the 
line  fit  to  UPL  data.  Specific  binding  was  calculated  as  the  difference  between 
total  binding  and  non-specific  binding.  The  affinity  of  each  receptor  for  [3H]TCDD 
(equilibrium  dissociation  constant,  Kd)  and  maximal  binding  (Bmax)  were 
determined  by  non-linear  regression  assuming  one-site  binding  in  Prism  4.0 
(GraphPad). 

2. 6  Transfection  Assays 

COS-7  (monkey  kidney)  cells  have  been  shown  to  lack  endogenous  AHR 
expression,  and  only  minimally  express  ARNT  (Ah  Receptor  Nuclear 
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Translocator)  (Ema  et  al.  1994a;  Ema  et  a/,  1994b).  We  transiently  transfected 
COS-7  cells  with  an  AHR  construct  {NARW,  HW,  BW,  HU,  or  MO  AHR),  human 
ARNT,  a  firefly  luciferase  reporter  construct  (pGudLuc6.1),  and  a  Renitla 
luciferase  control  plasmid  (pRL-TK).  The  pGudLuc6.1  plasmid  contains  the 
firefly  luciferase  gene  with  an  MMTV  promoter  under  the  control  of  dioxin 
response  elements  derived  from  the  mouse  Cyplal  gene  (Long  et  al.  1998).  A 
total  of  300  ng  of  plasmid  DNA  was  cotransfected:  50  ng  ARNT,  20  ng 
pGudLuc6.1 , 3  ng  pRL-TK,  AHR  (NARW:  0.5,  1 , 5,  1 0,  or  25  ng;  HW:  0.5,  1 , 5, 

10,  or  25  ng;  BW:  0.5,  1, 5, 10,  or  25  ng;  HU:  50  ng;  MO:  25  or  10  ng),  and  empty 
vector  pcDNA3.1  (177  -  227  ng). 

Ceil  culture  and  transfection  were  performed  in  48-well  plates,  with  40,000 
cells/well  in  complete  medium  (Eagle’s  Medium  with  fetal  bovine  serum).  After 
24  hours  of  growth,  the  medium  was  replaced  with  serum-free  medium  and  the 
cells  transfected  with  Lipofectamine  2000  (Gibco)  according  to  manufacturer's 
directions.  Five  hours  later,  cells  were  dosed  with  DMSO,  TCDD  in  DMSO  (10 
nM  final  concentration),  or  benzo[a]pyrene  in  DMSO  (BaP;  10  uM  final 
concentration),  and  incubated  for  18  hours.  Cells  were  collected  and  lysed,  and 
luciferase  activity  was  measured  using  the  Dual  Luciferase  Reporter  Assay 
System  (Promega)  and  a  Turner  TD-20/20  luminometer  (Turner  Designs). 
Relative  luciferase  units  (RLU)  were  determined  by  the  ratio  of  emission  of  firefly 
versus  Renilla  luciferase  in  the  presence  of  substrate.  Data  were  analyzed  by 
two-way  ANOVA  in  Prism  4.0. 

3.  RESULTS 

3. 1  NARW  and  HW  Aryl  Hydrocarbon  Receptor  Cloning 

An  initial  fragment  of  AHR  cDNA  was  obtained  by  RT-PCR  using 
degenerate  primers  and  mRNA  derived  from  NARW  biopsies  (Chapter  2).  We 
acquired  additional  5’  coding  sequence  from  NARW  5’  RACE  Ready  cDNA 
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synthesized  from  integument  mRNA  (SMART  RACE;  MM3r/UPM,  then  nested 
RWIr/  NUPM);  however,  reactions  failed  to  progress  5’  of  nucleotide  55  (coding 
valine  19).  Other  researchers  have  had  similar  difficulties  with  5’-RACE 
extension  in  the  same  region  of  fish  AHR  cDNAs  that  were  overcome  using 
adaptor-ligated  cDNAs  made  from  heat-denatured  RNA,  and  high  temperature 
polymerases  like  Thermoscript  RT  (Karchner  et  al.  1999),  However,  these 
methods  were  unsuccessful  for  NARW  AHR.  The  region  3'  of  the  initial 
fragment,  including  a  short  UTR  and  polyadenylation  signal,  was  identified  using 
3'  RACE  Ready  cDNA  and  specific  primers  generated  from  the  initial  fragment 
(RW-2f/UPM  then  nested  MM-4-F/NUPM).  RACE  reactions  with  Thermoscript 
RT-derived  cDNAs  and  Marathon  RACE  did  not  yield  any  visible  products. 

We  designed  a  NARW  AHR-specific  reverse  primer  (RW  AHR-1 R)  and 
used  it  with  a  forward  primer  specific  to  the  beluga  whale  AHR  (Bel-nat)  to 
amplify  the  full-length  NARW  AHR  cDNA.  The  Bel-nat  primer  includes  the  first 
nine  coding  bases  from  BW  AHR,  all  of  which  are  highly  conserved  in 
mammalian  AHRs  (identical  nucleotides  in  beluga,  human,  and  guinea  pig;  a 
single  difference  in  rabbit  AHR).  This  primer  combination  resulted  in  full-length 
amplification  of  AHR  from  NARW  3’  RACE  Ready  cDNA  (SMART  RACE),  as  well 
as  from  HW  cDNA  (Omniscript).  Amplicons  were  gel  extracted,  purified,  cloned, 
and  multiple  clones  (6)  were  sequenced  bi-directionally.  Consensus  sequences 
for  NARW  and  HW  AHR  cDNAs  were  determined  by  alignment  of  clones,  An 
open  reading  frame  (ORF)  encoding  850  amino  acids  (95,8  kD  predicted  mass) 
was  identified  for  NARW  AHR  (Figure  1 ;  GenBank  Accession  #EF091829),  The 
HW  AHR  was  similar,  with  an  851  amino  acid  ORF  (95.9  kD;  Figure  2;  GenBank 
Accession  #EF091830). 

Amino  acid  sequences  for  NARW  and  HW  AHRs  were  aligned  with  those 
from  BW  (Jensen  and  Hahn  2001),  human  (Dolwick  et  al.  1993a),  and  mouse 
(Burbach  et  al.  1992)  using  ClustalW  in  MacVector  8.1.1  (Figure  3).  The  high 
degree  of  amino  acid  conservation  throughout  the  entire  coding  region  enabled 
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the  identification  of  functional  domains  established  from  studies  of  mouse  AHR 
(Burbach  et  al.  1992;  Dolwick  et  al.  1993b;  Fukunaga  et  al.  1995;  Poland  et  al. 
1994).  Pairwise  comparisons  of  nucleotide  and  amino  acid  sequences  for  each 
AHR  are  shown  in  Table  2.  The  basic  helix-loop-helix  (bHLH)  region  is  very 
highly  conserved  among  these  mammalian  AHRs  {95.6  -  100%  amino  acid 
identity),  and  contains  a  nuclear  localization  signal  important  for  binding  co¬ 
factors  (HSP90  and  ARNT)  as  well  as  DNA  response  elements  (Ikuta  et  al. 

1998).  The  PAS  region  is  a  defining  feature  of  proteins  in  this  family  (Period- 
ARNT-Single-minded)  and  encompasses  much  of  the  AHR  ligand  binding 
domain  (LBD).  Whale  AHRs  share  98.3%  amino  acid  identity  in  the  PAS,  and 
slightly  less  conservation  is  observed  between  whales  and  human  or  mouse 
(86.9  -  92.1  %).  Only  three  residues  within  the  LBD  differ  among  whales  (K/R 
250;  M/L  347;  &  M/L  399).  A  high  level  of  conservation  is  also  seen  among 
whale  AHRs  in  the  carboxy-terminus  encoding  the  transactivation  domain 
(defined  here  as  amino  acid  490  (rel.  to  mouse  AHR)  to  the  stop  codon):  -89% 
mino  acid  identity  in  the  TAD,  compared  with  -74%  whales  to  human,  and  -56% 
whales  to  mouse.  An  extensive  analysis  of  amino  acid  sequence  differences  and 
their  putative  functional  effects  in  marine  mammal  versus  other  mammalian 
AHRs  are  discussed  in  Chapter  IV. 

3.2  Expression  Vectors 

Sequence  alignments  from  multiple  clones  also  identified  potential  PCR 
errors  in  each  clone.  A  HW  AHR  construct  of  the  consensus  sequence  was 
made  from  two  clones  by  endonuclease  digestion,  purification,  and  ligation  of 
error-free  sections,  then  subcloning  in  an  expression  vector  (pcDNA). 

The  same  method  was  not  possible  for  NARW  AHR  because  of  the 
locations  of  errors  relative  to  restriction  sites.  A  NARW  AHR  clone  with  a  single 
PCR  error  in  the  third  position  of  a  codon  (C2261  rather  than  consensus  A) 
resulted  in  a  non-synonymous  but  similar  amino  acid  (glutamine  (Q)  at  position 
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749  rather  than  lysine  (K)).  Site-directed  mutagenesis  was  attempted  to  alter  the 
sequence  to  match  the  NARW  AHR  consensus.  The  first  SDM  reaction  failed  to 
amplify.  Amplification  proceeded  in  our  second  SDM,  but  a  short  inverted  repeat 
within  the  primer  sequence  resulted  in  an  insert  of  three  concatenated  primer 
sequences  in  the  plasmid.  We  redesigned  the  primers  to  remove  this  motif,  but 
amplification  was  unsuccessful. 

We  decided  to  proceed  with  biochemical  assays  using  the  NARW  AHR 
construct  with  the  single  error;  first,  because  of  its  location  in  the  more  variable 
C-terminus.  Second,  the  side  chains  in  lysine  and  glutamine  are  both  polar  and 
tend  to  be  oriented  on  protein  surfaces.  Finally,  lysine  and  glutamine  have 
similar  molecular  weights  (128.18  and  128.14,  respectively),  which  in  total  limits 
the  effect  of  this  substitution.  Although  large-scale  deletions  or  insertions  have 
been  shown  to  affect  the  transcriptional  activity  in  other  AHRs  (Korkalainen  et  al. 
2000),  we  do  not  believe  the  single  point  mutation  present  in  our  NARW  AHR 
construct  has  an  impact  on  transcriptional  activation.  While  the  possibility 
remains  that  amino  acid  749  is  important  to  recapitulate  NARW  AHR  function, 
the  data  from  transient  transfections  in  COS-7  cells  shows  similar  activity  for 
NARW  and  HW  AHRs  (below,  section  3.6).  Further  work  to  correct  this  error  in 
the  plasmid  and  replication  of  biochemical  function  is  necessary  to  confirm  this 
assertion. 

3.3  Characterization  of  in  vitro  Expressed  Whale  AHR 

AHR  constructs  were  expressed  using  a  cell-free  in  vitro  transcription  and 
translation  system  (IVTT).  Whale  AHR  constructs  were  under  the  control  of  the 
T7  promoter,  whereas  human  and  mouse  AHRs  used  the  SP6  promoter.  Protein 
expression  was  confirmed  by  autoradiography  of  SDS-PAGE  separated  IVTT 
reactions  including  [35S]-labeled  methionine  (Figure  4A).  Gel  separation  was 
consistent  with  predicted  molecular  masses  of  95.8,  95.9,  95.4,  96.1,  and  90.3 
kDa  for  NARW,  HW,  BW,  HU,  and  MO  AHR,  respectively.  Bands  were  excised 
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and  counted  (disintegrations  per  minute;  dpm)  to  quantify  protein  expression. 
Variation  in  protein  levels  was  compared  in  four  independent  IVTT  reactions. 
Whale  AHRs  were  consistently  expressed  at  higher  levels  than  human  or  mouse 
AHRs;  this  may  have  been  a  consequence  of  vector  or  promoter  differences 
(whale  AHRs  in  pcDNA  (T7);  human  and  mouse  in  pSPORT  (SP6);  Figure  4B). 
However,  when  subcloned  into  pcDNA,  expression  of  human  and  mouse  AHRs 
remained  lower  than  that  of  whale  AHRs  (Figure  4C).  This  may  be  explained  by 
differences  in  post-translational  modifications,  or  sequence-specific  differences 
that  affect  the  processivity  of  the  IVTT  system. 

3.4  Specific  Binding  of[3H]TCDD  to  Whale  AHRs 

The  ability  of  IVTT  expressed  AHRs  to  specifically  bind  [3H]TCDD  (2nM) 
was  determined  by  velocity  sedimentation  on  sucrose  gradients  (Figure  5). 
Specific  AHR  binding  was  observed  with  peaks  between  fractions  10-19,  and 
was  differentiated  from  non-specific  protein  binding  by  comparison  with 
unprogrammed  lysate  (UPL;  IVTT  with  empty  pcDNA  vector).  UPL  has 
previously  been  shown  to  lack  specific  binding  to  [3H]TCDD  in  competition 
binding  experiments  (Jensen  and  Hahn  2001;  Karchner  et  al.  1999).  NARW  and 
HW  AHRs  sedimented  with  peaks  at  10. 7S  and  10. OS,  respectively.  Beluga 
whale,  human,  and  mouse  AHRs  sedimented  with  peaks  at  11. 3S,  10. 7S,  10. 7S, 
respectively.  Differences  in  specific  binding  maxima  may  be  related  to 
differences  in  the  amount  of  protein  per  reaction  as  observed  above,  but  also 
may  reflect  species-specific  differences  in  ligand  binding. 

3.5  Relative  Binding  Affinities  of  Right  Whale,  Humpback  Whale,  Beluga  Whale, 
Human,  and  Mouse  AHRs  by  Saturation  Binding  Analysis 

We  sought  to  compare  the  relative  binding  affinities  of  NARW  and  HW 
AHRs  with  high  and  low  affinity  AHRs  from  other  species.  Previous  studies  have 
described  the  binding  characteristics  of  the  beluga,  human,  and  mouse  AHRs 
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(Burbach  etal.  1992;  Dolwick  etal.  1993a;  Jensen  and  Hahn  2001;  Poland  et  al. 
1994),  Jensen  and  Hahn  (2001)  reported  that  the  mouse  (AH^1  allele)  and 
beluga  whale  have  high  affinity  AHRs  in  comparison  with  the  human.  We  used 
the  modified  saturation  binding  assay  with  dextran-coated  charcoal  (SBA/DCC) 
established  by  Jensen  (Jensen  and  Hahn  2000;  Jensen  and  Hahn  2001)  to 
determine  binding  affinities  of  NARW,  HW,  BW,  HU,  and  MO  AHRs  synthesized 
by  1VTT  (Figure  6). 

Specific  binding  was  calculated  as  the  difference  between  total  binding 
and  non-specific  binding  (UPL).  Nonlinear  regression  analysis,  equilibrium 
dissociation  constants  (Kd),  and  maximum  binding  (Bmax)  were  derived  assuming 
one-site  specific  binding  (Prism  4.0).  To  determine  the  optimal  time  of  incubation 
we  compared  binding  constants  calculated  from  2  hours  and  18  hours  incubation 
of  NARW  AHR  with  [3H]TCDD.  Longer  incubation  time  resulted  in  slightly  higher 
maximal  binding,  but  was  similar  to  that  from  2  hours  incubation  (Figure  7).  We 
decided  that  18  hours  incubation  was  both  logistically  easier  and  most  likely  to 
capture  the  receptor  behavior  in  the  assay  (whether  kinetic  or  dynamic). 
Ultimately,  binding  constants  and  maximum  binding  derived  from  2  and  18  hour 
incubations  were  within  one  standard  deviation  of  the  means  from  five  replicate 
assays  conducted  at  18  hours  (Table  3).  These  results  indicated  that  if  logistical 
constraints  limited  the  length  of  the  experiment,  similar  results  could  be  achieved 
with  shorter  incubations;  however,  we  recommend  a  time  course  of  assays  be 
conducted  before  using  a  shorter  incubation  time  on  AHRs  of  unknown  binding 
function. 

The  means  of  five  independent  SBA/DCC,  each  with  all  five  AHRs,  are 
shown  in  Table  3.  Differences  in  Bmax  are  correlated  with  the  amount  of  AHR 
protein  in  each  assay  (Figure  8;  ^=0.86).  However,  the  differences  in  binding 
observed  among  the  cetacean  AHRs  are  not  related  to  differences  in  AHR 
expression,  and  may  be  a  function  of  post-translational  processing,  or  species- 
specific  differences  in  amino  acid  sequence  and  function.  When  normalized  for 
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protein  content,  maximum  binding  is  similar  across  all  AHRs  investigated  (Table 
3,  Bmax  norm). 

3.6  Transcriptional  Activation  of  Reporter  Genes  by  Whale  AHR 

We  performed  transient  transfection  assays  to  investigate  the 
transcriptional  activities  of  NARW  and  HW  AHRs.  Mammalian  cells  deficient  in 
endogenous  AHR  (COS-7)  were  cotransfected  with  an  AHR  construct,  human 
ARNT,  a  firefly  luciferase  reporter  construct  (pGudLuc6.1),  and  a  Renilla 
luciferase  control  plasmid  (pRL-TK).  TCDD  induced  reporter-gene  activity  with 
the  MO  AHR  in  all  assays  (>2. 5-fold;  p<0.001 ),  but  HU  AHR  had  no  effect  on 
reporter  gene  transcription,  and  luciferase  activity  was  indiscernible  from  empty 
vector  controls  both  with  and  without  TCDD.  Our  laboratory  has  not  observed 
measurable  transcriptional  activity  with  this  HU  AHR  construct  in  transient 
transfections,  but  the  source  of  this  inactivity  has  yet  to  be  determined.  Other 
human  AHR  constructs  are  now  being  tested. 

Whale  AHRs  were  transcriptionally  active  in  the  presence  and  absence  of 
TCDD.  NARW,  HW,  and  BW  AHRs  produced  greater  reporter-gene  expression 
than  MO  AHR  at  the  same  plasmid  concentration  when  treated  with  either  DMSO 
or  TCDD  (25  ng  transfected;  Figure  9A).  No  significant  increase  in  reporter  gene 
expression  was  seen  with  whale  AHRs. 

Luciferase  activity  was  related  to  the  concentration  of  transfected  AHR. 
AHR  was  titrated  from  0.5  -  50  ng  to  establish  the  optimum  conditions  of 
detection  to  distinguish  background  from  induced  luciferase  reporter  expression. 
Luciferase  activity  for  the  highest  concentration  of  AHR  (50  ng)  for  whale,  but  not 
human  AHR,  exceeded  the  upper  limits  of  detection  for  the  luminometer.  Even 
lower  concentrations  of  whale  AHR  produced  luciferase  units  above  the  upper 
range  of  measurement,  necessitating  a  ten-fold  reduction  in  the  standard  volume 
read  on  the  luminometer.  As  little  as  0.5  ng  whale  AHR  was  sufficient  to  induce 
reporter-gene  transcription,  independent  of  ligand.  Still,  there  was  no  difference 
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in  the  transcriptional  activation  of  whale  AHRs  with  or  without  TCDD.  Although 
high  concentrations  of  transfected  plasmid  can  cause  constitutive  reporter  gene 
activation,  reduction  of  whale  AHRs  did  not  eliminate  constitutive  expression. 

While  whale  AHRs  were  not  sensitive  to  TCDD  under  the  conditions 
tested,  BaP  induced  luciferase  transcription  2.4-fold  with  BW  AHR  (1.0  ng 
transfected  only;  pO.OOl)  and  1.8-fold  with  MO  AHR  (p<0.05)  in  a  single 
experiment  (Figure  9B).  The  high  variation  observed  in  these  data  warrants 
additional  experiments  to  confirm  reporter  gene  induction  by  BaP;  however, 
these  preliminary  data  suggest  that  there  might  be  ligand-specific  differences  in 
transcriptional  activity  of  whale  AHRs. 

4.  DISCUSSION 

Direct  mortality  from  violent  interactions  with  ships  and  fishing  gear 
remains  the  primary  conservation  focus  for  NARW.  A  population  growth  model 
suggested  that  saving  the  lives  of  just  2  adult  females  annually  could  reduce  the 
risk  of  species  extinction  currently  predicted  in  about  200  years  (Fujiwara  and 
Caswell  2001).  However,  there  is  a  body  of  evidence  suggesting  that  NARW 
may  be  chemically  stressed,  potentially  contributing  to  the  slow  population 
growth  rate  by  disrupting  reproduction.  Our  principle  objective  was  to  determine 
if  NARW  possess  species-specific  characteristics  that  make  them  particularly 
sensitive  to  the  effects  of  chemical  contaminants  in  the  marine  environment. 
Specifically,  we  questioned  whether  environmental  contaminants  might  influence 
the  reproductive  function  of  the  species,  in  light  of  the  variable  fecundity 
observed  in  recent  decades.  To  approach  this  question,  we  first  sought  to 
determine  if  NARW  have  mechanisms  to  detect  and  respond  to  environmental 
pollutants  that  exert  toxicological  effects  via  the  AHR-mediated  pathway. 

In  a  previous  study  (Chapter  2),  we  showed  that  NARW  express  AHR  and 
CYP1A1  in  superficial  integument,  and  are  thus  capable  of  producing  a  response 
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to  AHR  ligands  such  as  PHAH  and  PAH.  Laboratory  studies  inform  our 
understanding  of  species-specific  sensitivity  to  AHR  agonists,  showing  a  causal 
relationship  between  in  vitro  ligand-binding  affinity  and  in  vivo  sensitivity  (Ema  et 
al.  1994b;  Karchner  et  al.  2006;  Poland  et  al.  1994;  Ramadoss  and  Perdew 
2004).  Two  early  experiments  found  that  inbred  strains  of  mice  with  low  affinity 
receptors  (AHd  allele)  were  less  sensitive  to  TCDD,  but  those  with  a  high-affinity 
receptor  (AH^1  allele)  were  highly  susceptible  to  the  physiological  effects  of  the 
chemical  (Birnbaum  et  al.  1990;  Poland  et  al.  1994).  Similar  types  of 
experiments  are  rarely  possible  with  marine  mammals,  and  are  impossible  for 
large  baleen  whales.  Jensen  and  Hahn  (2001)  established  that  beluga  whale 
AHR  function  could  be  investigated  using  an  in  vitro  biochemical  assay  thereby 
permitting  controlled  experimentation  for  a  class  of  animals  otherwise  difficult  to 
study.  Given  the  range  of  species-specific  responses,  we  sought  to  characterize 
the  in  vitro  biochemical  function  of  NARW  AHR  in  comparison  with  representative 
high-  and  low-affinity  AHRs  to  determine  where  it  may  be  on  the  continuum  of 
ligand  binding  affinity,  and  presumable  organismal  sensitivity. 


4. 1  Whale  aryl  hydrocarbon  receptor  sequences  are  highly  conserved 

NARW,  HW,  and  BW  AHRs  share  a  high  degree  of  amino  acid  sequence 
identity.  The  conservation  seen  among  whales,  human,  and  mouse  AHRs  also 
indicates  that  functional  domains  identified  in  the  latter  two  species  are  likely  to 
have  similar  functions  in  the  whale  AHRs.  The  amino-terminal  end  of  the  AHR 
protein  contains  the  bHLH,  PAS  A,  and  PAS  B  sequences  that  are  highly 
conserved  in  AHR  proteins.  Conversely,  large  sequence  divergence  has  been 
noted  for  the  carboxy-terminal  AHR  region  (Dolwick  et  al.  1993a;  Jensen  and 
Hahn  2001 ;  Ramadoss  and  Perdew  2004),  but  for  the  three  whale  species  this 
region  is  exceptionally  similar  (—89%).  While  striking  at  first,  these  data  conform 
to  expectations  based  on  evolutionary  divergence  of  the  Cetacea.  The  first 
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cetaceans  evolved  approximately  50  million  years  ago  (mya)  (Fordyce  and 
Barnes  1994;  Gingerich  etal.  2001),  and  combined  morphological  and  molecular 
data  sets  place  the  divergence  of  mysticetes  and  odontocetes  between  19  and 
34  mya  (Arnason  and  Gullberg  1996;  Messenger  and  McGuire  1998;  Milinkovitch 
1995);  balaenopterids,  including  HW,  only  diverged  in  the  last  7  my  (Demere  ef 
al.  2005).  For  comparison,  primates  and  rodents  diverged  approximately  96  mya 
(Nei  ef  al.  2001),  and  human  and  mouse  AHRs  share  only  68%  amino  acid 
identity,  whereas  chicken  and  tern  diverged  82-90  mya  and  maintain  92% 
identical  residues  throughout  the  AHR  (Karchner  ef  al.  2006).  The  few 
differences  that  exist  between  cetacean  AHRs  at  the  amino  acid  sequence  level 
are  analyzed  in  detail  in  Chapter  IV. 

4.2  Mysticete  AHRs  bind  dioxin  with  high  affinity 

In  this  study  we  systematically  compared  the  biochemical  functions  of  two 
newly  identified  mysticete  AHRs  and  a  previously  characterized  odontocete  AHR, 
with  those  from  human  and  mouse.  The  mouse  AHb1  allele  used  here  has  been 
shown  to  have  a  high  ligand  binding  affinity  (Poland  ef  al.  1994).  The  human 
AHR  represents  a  low-affinity  form  with  an  approximately  4-10  fold  lower  affinity 
than  the  mouse  AHR  (Ema  et  al.  1 994b;  Harper  et  al.  1 988;  Poland  ef  al.  1 994; 
Ramadoss  and  Perdew  2004).  A  statistical  difference  in  Kd  was  only  detected 
between  human  and  mouse  AHRs.  NARW  and  HW  AHRs  both  specifically  bind 
to  TCDD  with  a  relative  affinity  intermediate  to  those  of  mouse  and  human  AHRs. 
The  beluga  whale  AHR  was  previously  shown  to  have  a  TCDD-binding  affinity 
that  was  significantly  higher  than  that  of  the  human  AHR,  but  not  different  from 
that  of  the  mouse  (Ahb_1)  AHR  (Jensen  and  Hahn  2001).  In  our  hands,  BW  AHR 
was  also  statistically  indistinguishable  from  MO  AHR,  but  was  also  not  different 
from  HU  AHR  in  saturation  binding  assays  conducted  with  nearly  identical 
methods.  The  high  variability  in  Kd  values  among  assays  may  be  responsible  for 
the  lack  of  statistical  differences;  in  particular,  the  variability  of  the  human  and 
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humpback  whale  AHRs  in  these  experiments  reduced  our  ability  to  identify 
differences  among  receptors.  We  were  able  to  detect  statistical  significance  in 
the  ~3-fold  difference  between  mouse  and  human  AHRs,  therefore  the  <2-fold 
difference  between  whale  and  mouse  or  human  is  outside  the  power  of  the 
analysis. 

Whale  AHRs  can  also  activate  gene  transcription  of  a  reporter  gene  in  a 
cell  system.  Transcriptional  activity  was  observed  for  transiently  transfected 
whale  AHRs  independent  of  TCDD  exposure,  as  previously  shown  for  BW  AHR 
{Jensen  and  Hahn  2001).  Reduction  in  the  amounts  of  whale  AHR  constructs 
transfected  did  reduce  reporter  gene  expression,  but  did  not  yield  differentiation 
of  DMSO  and  TCDD  treatments.  A  possible  explanation  is  that  whale  AHR 
constructs  express  more  protein  than  the  mouse  in  the  COS-7  system,  resulting 
in  higher  protein  concentrations  with  subsequent  constitutive  transcriptional 
activation.  Quantification  of  protein  expression  in  COS-7  lysates  is  required  to 
address  this  hypothesis.  Mouse  AHR  consistently  activated  reporter  gene 
transcription  2.5  -  6.1  fold  in  a  TCDD-dependent  manner,  but  induction  was  only 
statistically  significant  for  assays  using  25  ng  plasmid  (2-way  ANOVA,  p  <  0.05  - 
0.001;  three  of  four  assays).  We  observed  benzo[a]pyrene-dependent  activation 
for  BW  AHR  but  not  the  other  whales,  suggesting  both  species-  and  ligand- 
specific  differences  in  transcriptional  activity.  BaP  also  induced  reporter  gene 
transcription  for  mouse  AHR.  Such  differences  may  be  related  to  differential 
adaptations  in  cetacean  suborders,  or  indicate  that  environmental  contaminants 
may  have  variable  effects. 

4.3  Whale  AHR  proteins  in  vitro 

The  in  vitro  expression  and  binding  capacities  of  whale  AHRs  differ 
dramatically  from  those  of  human  and  mouse  AHR.  Quantitative  measurement 
of  IVTT  indicated  that  whale  AHR  proteins  are  expressed  at  a  level  ~ 2-5  fold 
higher  than  that  of  human  or  mouse  AHR  constructs  used  in  saturation  binding 
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assays  (Table  3).  A  strong  correlation  exists  between  the  amount  of  AHR  protein 
and  maximal  binding  detected  in  these  assays  (Figure  8).  One  explanation  might 
be  different  efficiencies  of  transcription  and  translation  related  to  different 
promoters  in  the  two  expression  vectors  (whale  AHRs  in  pcDNA  (T7),  human  and 
mouse  AHR  in  pSPORT  (SP6)).  However,  this  difference  was  only  minimally 
alleviated  when  mouse  and  human  AHRs  were  subcioned  into  the  same 
expression  vectors  (Figure  4C),  suggesting  that  there  is  a  sequence-specific 
influence  over  in  vitro  expression.  The  length  of  AHR  sequence  was  probably 
not  a  factor  limiting  protein  expression  given  that  whale  AHRs  are  long  proteins 
(845  -  851  a. a.)  and  were  highly  expressed,  while  mouse  AHR  was  the  shortest 
protein  we  used  (805  a. a.)  and  expressed  the  least  protein  in  the  pcDNA  vector. 
Furthermore,  the  similar  activities  of  the  whale  AHRs  leads  us  to  believe  there  is 
something  intrinsic  to  whale  AHR  proteins  such  that  they  express  better  in  vitro. 

It  is  unknown  whether  this  factor  increases  transcription,  translation,  or  the 
stability  of  the  whale  AHR  proteins.  The  question  remains:  why  are  some  AHR 
proteins  more  highly  expressed  than  others  in  vitro ?  The  answer  to  this  question 
may  be  important  in  accurately  applying  in  vitro  biochemical  methods  to 
understand  mechanisms  of  action  of  AHR  agonists. 

The  difference  in  protein  content  may  impact  the  results  of  our  binding 
assays.  In  vitro  ligand  binding  is  sensitive  to  protein  concentrations,  and  high 
protein  concentrations  may  decrease  the  apparent  binding  affinities  of  AHRs 
(Bradfield  etal.  1988).  Although  we  show  large  differences  in  AHR  expression 
among  receptors,  the  total  protein  content  of  each  IVTT  reaction  was 
approximately  equal  therefore  mostly  avoiding  the  issue  developed  by  Bradfield 
et  al.  (1998).  Another  study  illustrated  problems  associated  with  in  vitro 
expression  and  high  protein  concentrations  (Ramadoss  and  Perdew  2004). 

Given  the  high  expression  of  our  whale  AHRs,  we  may  be  underestimating  the 
ligand  binding  affinities  of  these  receptors  in  these  experiments.  This  also  may 
explain  the  difference  in  beluga  AHR  affinity  in  this  study  compared  with  Jensen’s 
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initial  characterization  that  showed  a  significant  difference  between  human  and 
both  BW  and  mouse  AHRs.  Reducing  the  amount  of  whale  AHR  used  in  the 
binding  assays  by  dilution  with  UPL  could  reduce  the  apparent  binding  constants, 
and  would  perhaps  allow  statistical  differentiation  of  whale  AHRs  from  human. 

Another  explanation  for  differences  in  protein  expression  may  lie  in  the 
sequences  themselves.  Recent  work  on  the  MDR1  in  HeLa  cells  has  suggested 
that  codon  usage  may  impact  protein  function  by  altering  the  rate  of  co- 
translational  folding  (Kimchi-Sarfaty  et  al.  21  December  2006).  The  prevalence 
of  rare  codons  in  a  cDNA,  even  with  conserved  amino  acid  sequence,  may 
further  influence  the  rate  of  protein  construction  at  a  more  simplistic  level.  A 
comparison  of  codon  usage  for  within  the  constructs  used  here  remains  to  be 
completed. 

4.4  Summary 

The  framework  developed  by  the  methods  presented  in  this  study  offers  a 
means  for  comparing  the  primary  mechanism  of  sensitivity  to  DLC  in  cetaceans 
(i.e.  the  AHR).  We  have  shown  that  mysticete  AHRs  bind  to  TCDD  as  well  as 
that  from  BW,  and  with  binding  constants  intermediate  to  those  of  the  high-affinity 
mouse  and  lower-affinity  human  AHRs.  Cetacean  AHRs  are  capable  of 
transcriptional  activity  in  a  reporter-gene  cell  system,  but  further  work  is  required 
to  determine  whether  whale  AHRs  activate  gene  expression  in  a  ligand- 
dependent  manner.  Competition  binding  studies  are  needed  to  determine 
structure-activity  relationships  of  various  AHR  ligands  and  test  whether 
environmentally  relevant  contaminants  may  have  more  important  effects  than 
studies  using  TCDD  might  indicate,  especially  given  our  results  with  BaP  in 
transfection  assays.  Finally,  we  emphasize  that  understanding  AHR  activity  in 
whales  is  only  one  small  part  of  the  complex  mechanism  of  response  to 
environmental  contaminants,  which  may  be  influenced  by  many  other  proteins, 
hormones,  and  synergistic  activities  in  this  complex  pathway. 
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Table  1.  Oligonucleotide  primers  for  aryl  hydrocarbon  receptors. 

Primer_ Sequence  (5'  ->  3’) Tm 

MM-3-r  GAR  GAG  TTY  TCT  GGA  GGA  AGC  TG  65.1  °C 

MM-4-F  GTT  TTA  GGC  TAT  ACT  GAA  GCA  GAG  59.2°C 

RW-lr  CAA  AGC  CAT  TCA  GTG  CCT  GCA  GTA  A  71.5°C 

RW_PB-5f  GCC  ACC  ATC  CAT  ACT  TGA  AAT  CCG  70.7°C 

Bel-nat  CCC  AAG  CTT  GGG  CAC  CAT  GAA  CAG  CAG  C  81 .6°C 

RW  AHR- 1 R  AAC  CAA  GAT  GAA  AAG  TGG  GC  62.8°C 

RW  C2261 A-F  CCT  CAA  AAC  CAA  [a]AG  CAC  GGA  CTC  AAT  CCA  CAG 

TCA  GCC  C  84.8°C 

RW  G2261T-R  GGG  CTG  ACT  GTG  GAT  TGA  GTC  CGT  GCT  [t]TT  GGT 

TTT  GAG  G  84.8°C 

RW  C2261A2-F  CCA  A[a]A  GCA  CGG  ACT  CAA  TCC  ACA  GTC  AGC  CC  81.6°C 

RW  G2261T2-R  GGG  CTG  ACT  GTG  GAT  TGA  GTC  CGT  GCT  [t]TT  GG  81 ,6°C 
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Table  2.  Comparison  of  mammalian  aryl  hydrocarbon  receptors:  percent 
sequence  identity. 

Nucleotide  and  amino  acid  sequences  were  aligned  separately  using  ClustalW 
1 .8  in  MacVector  8.1.1 .  Pairwise  comparisons  are  shown  for  full-length  RW,  HW, 
BW,  HU,  and  MO  AHRs.  Percent  nucleotide  identity  is  shown  above  the 
diagonal,  and  amino  acid  identity,  below. 


RW  AHR 

HW  AHR 

BW  AHR 

HU  AHR 

MO  AHR 

RW  AHR 

97,7 

95.7 

85.7 

72.1 

HW  AHR 

97.1 

95.5 

85.8 

72.1 

BW  AHR 

94.0 

93,8 

86.3 

72.7 

HU  AHR 

83.5 

83.0 

83.2 

73.1 

MO  AHR 

68.5 

68.1 

67.9 

68.3 
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Table  3.  Dissociation  binding  constants  (IQ)  and  maximal  binding  (Bmax) 
determined  for  IVTT  expressed  aryl  hydrocarbon  receptors. 

The  protein  content  for  each  100  pL  reaction  (mean  +  SD)  was  determined  by 
[35S]-iabeled  methionine  IVTT/SDS-PAGE  (n=4;  Fig.  4B).  Kd  and  Bmax  (means  + 
SD)  were  calculated  from  saturation  binding  assays  of  IVTT  AHRs  with  [3H]TCDD 
(n=5).  In  a  one-way  ANOVA  of  mean  KDs,  only  mouse  and  human  AHRs  were 
statistically  different  (£,  Tukey’s  post-hoc  test,  p<0.05).  Maximum  binding  (Bmax) 
was  normalized  taking  into  account  differences  in  AHR  protein  content  for  each 
receptor. 


fmol  AHR 

in  100  nL  IVTT 

Kd  +  SD 
(nM) 

Bmax  i  SD 
(nM) 

Bmax 

(norm) 

NARW AHR 

828  +  160 

0.64  +  0.19 

0.41+0.14 

0.42 

HW  AHR 

852  +  376 

0.85  +  0.44 

0.65  +  0.24 

0.65 

BW  AHR 

804  + 362 

0.78  +  0.22 

0.49  +  0.09 

0.52 

HU  AHR 

172  +  35 

0.97  +  0.46 

0.09  +  0.04 

0.45 

MO  AHR 

425  +  83 

0.36  +  0.13  * 

0.27  +  0.10 

0.54 
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ATGAACAGCAGC  GGCGCCAACATC  ACCTACGCCAGT  CGCAAGCGGCGC  AAGCGGGTGCAG  AAAACTGTCAAG 
M  N  S  S  G  ANI  T  Y  A  S  R  K  R  R  KPVQ  K  T  V  K  24 

73  CCAGTCCCAGCT  GAAGGAATCAAG  TCGAATCCTTCA  AAGCG  GC  AT  AG  A  GACCGACTTAAT  ACCGAATTGGAC 

P  V  P  A  E  G  IK  SNPS  K  R  H  R  DRLN  TELD43 
14  5  CGTTTGGCCAGC  CTGCTGCCTTTT  CCACAAGATGTT  ATTAATAAGCTG  GACAAACTTTCA  GTTCTTAGGCTC 

R  L  A  S  LLPF  PQDV  INKL  DKLS  V  L  R  L  72 

217  AGTGTCAGTTAT  CTAAGAGCCAAG  AGCTTCTTTGAT  GTTGC AT TACAG  TCTACTCCAGCT  G ACCG AAATGG A 

SVSY  LRAK  SFFD  VALG  S  T  P  A  DRNG96 

289  GTCCAGGACAAC  TGTAGAACAAAA  TTCAGAGAAGGC  C TG AACT TGCAG  GAAGGAGAATTC  TTACTGCAGGCA 

V  g  D  N  CRTK  F  R  E  G  L  N  L  0  EGEF  L  L  0  A  120 

361  C TG AATGGCTTT  GTATTGGTTGTC  ACTACAGATGCT  TTGGTCTTTTAT  GCTTCT  TCTACT  ATACAAGATTAC 

L  N  G  F  VLVV  TTDA  LVFY  ASST  X  Q  P  Y  144 

4  33  CT GGGGTT TC AG  CAGTCTGATGTC  ATCCATCAGAGT  GTGTATGAACTG  ATCCATACTGAA  GACCGAGCTGAA 

LGFQ  OSDV  IHQS  V  Y  E  L  I  H  T  E  D  R  A  E  168 

505  tttcagcgccag  ctacactgggca  ttaaacccttca  cagtgtccagac  TCTGGACAAAAA  ATTAATGAAGCT 

f  g  r  o  lhwa  lnps  q  c  p  d  sgqk  I  N  E  A  192 

57  7  AATGGCCTCTCA  C AGCCAGCAGTC  T ATT AT AACC C A  GACCAGCTTCCT  CCAGAGAACTCC  TCGTC TATGGAA 

NGLS  QPAV  YY  NP  DQLP  PENS  S  S  M  E  21.6 

64  9  AGGTGCTTCGTT  TGCCG AC T AAGG  TG TTTG CTGGAT  AAT TC ATCTGGT  TTTCTGGCAATG  AATTTCC AAGG G 

RCFV  CRLR  C  L  L  D  NSSG  FLAN  N  F  Q  G  240 

72 1  AGGTTGAAGTAT  C TTC ATGGAC AG  AACAAGAAAGGG  AAAGATGGATCA  ATACTTCCACCT  CAGTTGGCTTTG 

R  L  K  ¥  LH  GO  NKKG  KDG  S  ILPP  QLAL264 

793  TTTGCAATAGCT  AC TCC AC TGCAG  CCACCATCCATA  CTTGAAATCCGA  ACCAAAAATTTC  ATCTTTAGAACC 

F  A  I  A  TPLQ  PPSI  LEIR  TKNF  IFRT  288 

865  AAACACAAGCTA  G AC  TTTAC ACCT  AC  TGGTTG  TG AT  GCC  AAAGGAAG A  ATTGTTTTAGGC  TAT AC  TG AAGC  A 

KH  KL  D  F  T  P  TGCD  AKGR  IVLG  YTEA312 

937  G AGC T ATGCATG  AGAGGATCAGGA  TATC AATT T ATT  CAT GC TGCT GAT  ATGCTTTA TTGT  GCTGAGTACCAT 

ELCW  RGSG  YQFI  HA  AD  HLYC  A  E  Y  H  336 

1009  ATCCGGATGATT  AAGACTGGAGAG  AGTGGCATGATA  GTGTTC AGC  CTT  CTTACCAAAGAC  AATCG AT  GG ACT 

IRHI  KTG  E  SGML  VFRL  LTKD  NRWT360 

1081  TGGG  TGCAGTC A  AATGC AC  GC  TT A  GTTTATAAGAAT  GGAAGACCAGAT  TATATCATTGCA  ACTCAGAGACCT 

WVQS  NARL  V  ¥  K  N  GRPD  YIIA  TQRP384 

1153  CTAACAGATGAA  GAAGGAACAGAG  CATCTACGAAAA  CGAAATATGAAG  TTGCCTTTTATG  TTTACCACTGGC 

ltde  EGTE  HLRK  RNNK  LPFM  FTTG  408 

1225  GAAGCTGTTTTA  TATGAGGTAACC  AACCCTTTTCCT  CCCATAATGGAT  CCCTTACCAATA  AGG ACT AAAAC T 

EAVL  Y  6  V  T  NPFP  PI  MD  PLP1  RTKT432 

1297  GGTGCTGGTGGA  AAAGGTTCTGCT  ACCAAATCAACT  CTAAATAAGGAT  TCCCTCAATCCC  AGTTCCCTCCTG 

GAGG  KG  S  A  T  K  S  T  L  N  K  D  £  L  N  F  S  S  L  L  456 

1369  AATGCCATGATG  CAACAAGATGAA  TCTATTT ATC TC  T ATCC TGCTTC A  AGTAGTACACCT  TTTGAAAGAAAC 

NANM  Q  Q  D  E  S  IYL  YPAS  S5TP  FERN  4B0 

1441  TTTTTCACTGAC  TCTCAGAATGAG  TGCAGTAGTTGG  CAAAGCAATGTT  GCACCAATGGGA  AGTGACAATATC 

F  FTD  S  Q  N  E  CSSW  0  $  N  V  APMG  SDNI504 

1513  CTGAAACATGAG  C AG AT  TGGC  CAG  TCTCAGGAAATG  AACCCCACAGTC  TCCG  GAGATC  AC  GCG  GGGCTC  TT  T 

LKHE  QIG  Q  S  0  E  M  N  P  T  V  S  G  D  H  A  G  L  F  52  8 

1585  CC  AG AT AAT AGA  AAT AG  TG AC  TTG  TACAGCATTATG  AAACA CCTAGGC  ATTGATTTTGAA  GAT  ATCAAACAC 

PDNR  NSDL  YSIM  KHL  G  IDEE  DIKH552 

1657  ATGCAACAGAAT  G AGGAAT TTTT C  AGAACTGACTTT  T CTGGTG AAGAT  GACTTCAGAGAT  ATTGACTTAACA 

MQQN  EEFF  RTDF  S  G  E  □  D  FRD  1DLT  576 

1729  GACGAAATTCTG  ACCTATGTCGAA  GACTCTTTAAAT  AAGTCTGCCTTT  GGCTGTTCAGGT  TGCCATCAGCAA 

DEIL  TYVE  D  S  L  N  K  £  A  F  GCSG  C  H  g  g  600 

1801  CAGTCCATGGCT  CTGAACCCCAGC  TGC ATG  GT AC  AG  GAG  C ACC  TCCGG  TTAGAACAGCAA  G AGCAGCT ACAG 

Q  5  M  A  LNPS  CNVQ  EHLR  L  E  Q  Q  B  0  L  Q  624 

1873  CACCATCAGCAG  CACACAGCAGCA  GACCAGCAACAG  CTGTGTCAGAAA  ATGAAGCACATG  CAAGTTAACGGC 

HHQC  H  T  A  A  EQQQ  L  C  0  K  MKHN  0  V  N  G  648 

194  5  ATGTTTGCAAAC  TGGAACTCCGAC  CAATCCGTGCCT  TTTAACTGTCCT  CAGCAAGACTTA  C AGCAGT ACAGC 

MEAN  WNS  D  O  S  VP  FNCP  Q  Q  D  L  0  0  Y  S  672 

2017  CTCTTTTCAGAC  TTACCTG  GG AC  C  AG  TC AAGAGTTT  CCCTACAAACCC  GAGTTTGATCCT  ATGCCTT  ACGCA 

VFSD  LP  GT  SOEF  PYKA  EFDP  MPYA696 

2089  C AG AAC TTTATC  CCC TGTAGTC AG  TC TGTG TTGCCG  CCACATTCTAAG  GGTACACAGTTA  GACTTTCCCATT 

QNF1  PCSO  SVLP  PHSK  GTQL  DFFI720 

2161  GGGGATTTCGAA  CCACCCCCATAC  CCTACAACTTCT  TCTCATTTAGAA  GATTTTGTCACA  TGTTTACAAGTT 

GDFE  PAPY  PTTS  SHLE  DFVT  CLQV744 

2233  CCTCAAAACCAA  AAGCACGGACTC  AATCCACACTCA  GCCCTACTAACT  CCTCAGACATGT  TATGCGGGGGCT 

PQNQ  KHGL  NP  OS  ALLT  POTC  YAG  A  76B 

2305  GTGTCAATGTAC  CAGTGCCAGCCG  GAAGCTCAGCAC  AGCCATGTG GC T  CAGATGCAGTAC  AACCCGACGATG 

VSHY  QCQP  EACH  SHVA  Q  H  Q  Y  NPTN  792 

2377  CCAGGCCCGACG  ATGCCAGGCCCG  C AGGC ATT  TTT A  AACAAGTTTCAG  AAC  G  G AGG  AG  TC  TTAAATGAAACC 

PGPT  MPGP  QAFL  NKFQ  NG  GV  LNET816 

2449  TATCCACCTGAA  TTAAATAGCATA  AATAACACTCAG  CCTACCGCCCAT  CTTCACCCGTCA  GAAGCCAGACCT 

YPAE  LNSI  NNTQ  PTAH  LHPS  EARP84  0 

2521  TTCCCTGACTTG  ACATCCAGTGGA  TTCC  TGT AAT  TC  CAAGCCCACTTT  TCATCTTGGTTT  TTTGATTAAATT 

FPDL  TSSG  FL*  850 

2593  CGTTTG TG AAAC  ATCGAAAAAAAA  AAAAAAAAAAAA  AAAAAAAAA 

Figure  1 .  Full  length  north  Atlantic  right  whale  aryl  hydrocarbon  receptor  cDNA. 

The  nucleotide  and  predicted  amino  acid  sequences  (numbered  left  and  right,  respectively)  for 
the  coding  sequence  of  NARW  AHR.  The  stop  codon  is  indicated  with  an  asterisk  (*). 
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1 

AT  GAAC  AGC AGC 

M  H  S  S 

73 

CCAGTCCCAGCT 

P  V  P  A 

14  5 

CG TTTG GCC AGC 

R  L  A  S 

217 

AG  TGTC  AGTTAT 

S  V  S  Y 

289 

GTCCAGGACAAC 
V  Q  D  H 

361 

CTGAATGGCTTT 
L  N  G  F 

433 

CTGGGGTTTCAG 
L  G  F  Q 

505 

TTTCAGCGCCAG 
F  0  R  Q 

577 

AATGGCCTCTCA 

N  G  L  S 

649 

AGGTGCTTCG  TT 

R  C  F  V 

721 

AGGTTGAAGTAC 

R  L  K  Y 

793 

TTTGCAATAGCT 

F  A  I  A 

865 

AAACACAAGCTA 
K  H  K  L 

937 

GAGCTATGCATG 
E  L  C  M 

1009 

ATCCGGATGATT 

1  R  M  I 

1081 

TGGGTGCAGTCA 
W  V  Q  S 

1153 

CTAACAGATGAA 
L  T  D  E 

1225 

GAAGCTGT  TTT A 

E  A  V  L 

1297 

GG  TGCTGG  TGGA 

G  A  G  G 

1369 

AATGCCATGATG 
N  A  M  M 

1441 

TTTTTC  AG  TG AC 

F  F  S  D 

1513 

CTGAAAC  ATG AG 

L  K  H  E 

1585 

CCAGATAATAGA 

POUR 

1657 

ATGCAACAGAAT 

M  Q  Q  N 

1729 

GATGAAATTCTG 
D  E  I  L 

IB  01 

CAGTCCATGGCT 
0  S  M  A 

1873 

CACCGTCAGCAG 

H  R  0  Q 

1945 

GGCATGTTTGCA 
G  M  F  A 

2017 

AATGTC  TTTTCA 
N  V  F  S 

2089 

GC  ACAG  AACTTT 
A  Q  N  F 

2161 

AT  TGGG  C ATTTC 

I  G  H  F 

2233 

GTTCCTCAAAAC 
V  P  0  N 

2305 

GCTGTGTCAATG 
A  V  S  H 

2377 

ATGCCAGGCCCA 
M  P  G  P 

2449 

ACCT AT  CC AGCT 
T  Y  P  A 

2521 

CCTTTCCCTGAC 

P 

F  P  0  L 

A GCGCCAACATC 
SAMI 
GAAGGAATCAAG 
E  G  I  K 
CTGCTGCCTTTT 
L  L  P  F 
CTAAG AG  CC AAG 
L  R  A  K 
TGTAGAACAAAA 
C  R  T  K 
GTATTGGTTGTC 

Y  L  V  V 
CAGTCTGATGTC 

Q  S  D  V 
CTACACTGGGCA 
L  H  W  A 
C AGCC  AGCAGTC 
Q  P  A  V 
TGCCGACTAAGG 
C  R  L  R 
C  TTCATGG AC AG 
L  H  G  0 
ACTCCACTGCAG 
T  P  L  Q 
GACTTTACACCT 
D  F  T  P 
AGAGGATCAGGA 
R  G  S  G 
AAG AC  TGGAG AG 
K  T  G  E 
AATGC ACGC  TT A 
N  A  R  L 
GAAGGAACAGAG 
E  G  T  E 
TATGAGGTAACC 

Y  E  V  T 
AAAGGTTCTGCT 

K  G  S  A 
CAACAAGATGAG 
CODE 
TCTCAGAATGAG 
S  Q  N  E 
CAGATTGGCCAG 
0  I  G  Q 
AATAGTGACTTG 
H  S  0  L 
GAGGAATTTTTC 
E  E  F  F 
ACTTACG  TCGAA 
T  Y  V  E 
CTGAACCCCAGC 
L  H  F  6 
C AG  AT AG  CAGCG 
H  I  A  A 
AACTGGAACTCC 
N  W  N  S 
GACTTACCTGGG 
D  L  P  G 
ATC  CCCTGT AGT 
IPGS 
GAACCAGCCCCA 
E  P  A  P 
C AAAAGC  ACGG A 
Q  K  H  G 
TACCAGTGCCAG 

Y  0  C  Q 
ACGGTGC  CAGGC 

T  V  P  G 
GAATTAAATAGC 
ELMS 
TTGAGATCCAGT 
T  S  S  C 


ACCT AC  GCCAGT  C  GC AAGC  GGCG  G 
T  Y  A  S  R  K  R  R 

TCGAATCCTTCA  AAGCGGC AT AG  A 
SNPS  KRBR 
CC  AC  AAG ATG  TT  ATTAATAAGCTG 
P  Q  D  V  I  N  K  L 

AGCTTC  TTTG AT  G  TTGC  ATT ACAG 
SFFD  VALQ 
TGCAGAGAAGGC  CTGAACTTGCAG 
C  R  E  G  L  N  L  Q 

AC  T AC AG ATGCT  T  TGGTCTTTTAT 
TTDA  L  V  F  Y 

ATCCATCAGAGT  G TGTATG AACTG 
IHQS  VYEL 
GTAAACCCTTCA  CAGTGTCCAGAC 

V  N  P  5  Q  C  P  D 

TATTATAACCCA  GAG  CAGGTT  CCT 

V  Y  N  P  D  Q  V  P 

TGTCTGCTGGAT  AATTCATCTGGT 

C  L  L  D  N  5  5  G 

AACAGGAAAGGG  AAGGA  TGGATC A 
NRKG  KPGS 
CCACCGTCCATA  CTTGAAATCCGA 
PPSI  LE1R 
ACCGGCTGTGAT  GCCAAAGGAAGA 
TGCD  AKGR 
tatcaatttatt  catgctgctgat 

YQFI  H  A  A  D 

AGTGGCATGATA  GTG  TTC AGGCTT 
SGM1  VFRL 
GTTTATAAGAAT  GGAAG AC  C AG AT 

V  Y  K  N  GRPD 

C ATT  T ACG AAAA  CGAAATATGAAG 

H  L  R  K  RNHK 

AACCCTTTTCCT  TCCATAATGGAT 
K  PFP  SIHD 
ACCAAATCAACT  CTAAATAAGGAT 
TEST  L  N  K  D 

TCTATTTA TC  TC  T AT CCTGCT TC  A 
S  I  Y  L  YPAS 

TGCAGTAGTTGG  CAAAGCAATGTT 
C  S  S  W  QSNV 

TCTCAGGAAATG  AACCCCACAGTC 
SQEM  UP  T  V 

TACAGCATTATG  AAACACCTAGGC 
YSIM  KHLG 
AGAACTGACTTT  TCTGGTGAAGAT 
RTDF  SGE£> 
GACTCTTTAAAT  AAGTCTGCCTTT 
D  S  L  N  K  S  A  F 

TGCATGGTACAG  G AGCACC  TCCAG 
CMVQ  E  R  L  0 

GAGCAGCAGCAA  CAGCTGTGTCAG 
E  Q  Q  Q  Q  L  C  Q 

GACCAATCCGTG  C  CT  TTTAACTGT 
0  0  S  V  P  F  N  C 

ACC AGT  CAAG  AG  TTTCCCTACAAA 
TSQE  FPYK 
C AGTCT  GTGTTG  CCGCCACATTCT 
Q  S  V  L  P  P  H  5 

TACCCTACACCT  TCT  TCTAATTT A 
YPTP  SSNL 
CTCAATCCACAG  TCAGCCC  TAGT A 
LNPQ  S  A  L  V 

ccggaagctcag  cacagccatgtg 

PEA  Q  HSRV 
CCGCAGGCGTTT  TTAAACAAGTTT 
P  0  A  F  LNKF 

ATAAATAACACT  CAGCCTACCGCC 
INNT  QPTA 
GGATTCCTGTAA 
F  L  *  851 


AAGCCGGTGCAG 
K  P  V  Q 
GACCGACT  T AAT 
D  R  L  n 
GACAAACTTTCA 
D  K  L  S 
TCCACTCCAGCT 
S  T  P  A 
GAAGGAGAATTC 
E  G  E  F 
GC  TTCTTCT ACT 
ASST 
ATCCACACTGAA 
I  H  T  E 
TCTGGACAAAAA 
S  G  OR 
CC AG AG  AACTCT 
PENS 
TTTCTGGCAATG 
FLAM 
AT ACTTCCGC  CT 
I  L  P  P 
ACC AAG AACT  TC 
T  K  N  F 
ATTGTTTTAGGC 
I  V  L  G 
ATGCTTTATTGC 
ML  Y  C 
CTTACCAAAGAC 
L  T  K  D 
TAT ATC  ATTGC A 
Y  I  I  A 
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Figure  2.  Full  length  humpback  whale  aryl  hydrocarbon  receptor  cDNA. 

The  nucleotide  and  predicted  amino  acid  sequences  (numbered  left  and  right,  respectively)  for 
the  coding  sequence  of  HW  AHR.  The  stop  codon  is  indicated  with  an  asterisk  (*). 
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Figure  3,  Alignment  of  aryl  hydrocarbon  receptor  amino  acid  sequences. 

Predicted  amino  acid  sequences  of  whale,  human,  and  mouse  AHR  coding  regions  used  in  this 
study  were  aligned  using  ClustalW  (1.8)  in  MacVector  8.1.1.  Identical  residues  in  three  or  more 
sequences  are  outlined  and  shaded;  similarities  are  in  bold. 
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Figure  4.  Whale  AHR  constructs  express  protein  by  in  vitro  transcription 
and  translation  (IVTT). 

A,  Right  whale  and  humpback  whale  AHR  coding  sequences  were  ligated 
into  pcDNA  expression  vectors  (T7  promoter)  and  expressed  by  IVTT  in  the 
presence  of  pS]-labeled  methionine;  unprogrammed  lysate  (UPL;  empty 
pcDNA  vector),  and  beluga  (pcDNA),  human,  and  mouse  AHR  constructs 
(pSPORT  vector;  SP6  promoter)  were  also  expressed.  Protein  was  detected 
and  quantified  using  SDS-PAGE  and  autoradiography.  8,  Variability  in  the 
amount  of  protein  expressed  in  individual  IVTT  reactions  was  assessed  in 
replicate  experiments  (means  +  SD,  n=4).  C,  Higher  expression  of  whale  AHRs 
was  not  completely  explained  by  different  expression  vectors;  when  subcloned 
into  pcDNA,  the  pattern  of  lower  expression  of  human  and  mouse  AHR 
expression  remained  (means  +  SD,  n=3). 
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Figure  5.  In  vitro  expressed  AHR  proteins  show  specific  ligand  binding  by 
velocity  sedimentation. 

Specific  ligand  binding  was  assessed  by  velocity  sedimentation  analysis  on 
sucrose  gradients  after  incubation  of  [JH]TCDD  (2  nM)  with  IVTT  expressed 
AHR  proteins.  Specific  binding  was  calculated  for  each  receptor  as  the 
difference  in  total  binding  and  non-specific  binding  (as  measured  by 
unprogrammed  lysate,  UPL).  The  elution  positions  of  sedimentation  markers 
[14C]ovalbumin  (3.6  S)  and  [14C]catalase  (1 1.3  S)  are  shown  with  left  and  right 
arrows,  respectively. 


108 


Figure  6.  Representative  binding  curves  for  aryl  hydrocarbon  receptors. 

Binding  curves  from  one  of  five  replicate  saturation  binding  assays  are  shown  for 
North  Atlantic  right,  humpback  and  beluga  whale,  human,  and  mouse  AHR. 

IVTT  expressed  AHR  proteins  were  diluted  in  buffer  and  incubated  with  one  of 
nine  or  ten  concentrations  of  [3H]TCDD.  Specific  binding  (open  squares)  was 
calculated  as  the  difference  between  total  binding  (closed  squares)  and  non¬ 
specific  binding  (open  circles;  UPL,  unprogrammed  lysate).  Nonlinear  regression 
analysis  and  binding  constants  were  derived  assuming  one-site  binding  in  Prism 
4.0.  (Y=Bmax*X/(Kd+X);  GraphPad). 
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FREE  [3H]TCDD  (nM) 


Figure  7.  Saturation  binding  assay  incubation  time  does  not  affect  binding 
constants. 

NARW  IVTT  AHR  protein  was  incubated  with  ten  concentrations  of  [3H]TCDD  for 
2  or  18  hours  to  determine  optimum  assay  conditions.  Greater  maximal  binding 
was  achieved  from  the  longer  incubation,  but  measures  of  Kd  and  Bmax  are  within 
one  standard  deviation  of  the  means  of  5  other  assays. 
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Figure  8,  Maximal  binding  is  correlated  with  the  amount  of  AHR  protein. 

Differences  in  Bma*  are  correlated  with  the  amount  of  AHR  protein  in  each  binding 
assay  (r2-  0.86).  The  same  volume  of  IVTT  AHR  (100  uL)  was  used  for  each 
saturation  binding  assay.  However,  each  construct  does  not  express  the  same 
amount  of  AHR  protein  (Table  3);  whale  AHR  expression  constructs  consistently 
express  greater  amounts  of  specific  protein. 
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Figure  9.  Whale  aryl  hydrocarbon  receptors  are  transcriptionally  active. 

Mammalian  cells  deficient  in  endogenous  AHR  expression  (COS-7)  were 
cotransfected  with  0.5  -  50  ng  NARW,  HW,  BW,  HU,  or  MO  AHR  expression 
plasmids,  plus  human  ARNT,  a  reporter-gene  (pGudLuc6.1 ),  and  a  transfection 
control  (pRL-TK)  in  48-well  culture  plates.  We  compared  transcriptional 
activation  of  each  AHR  with  DMSO  or  10  nM  TCDD  in  four  experiments 
(representative  results  shown  in  A).  MO  AHR  consistently  induced  reporter-gene 
expression  (2.5-fold)  in  the  presence  of  TCDD,  whereas  whale  AHRs  were 
transcriptionally  active  with  and  without  ligand.  A  single  assay  (B)  consisting  of 
two  plates,  one  dosed  with  TCDD  (10  nM),  the  other  with  benzofa] pyrene  (10 
uM),  and  each  with  separate  DMSO  controls,  produced  similar  results  for  TCDD 
(MO  AHR:  2.5-fold  induction);  BaP  induced  reporter-gene  expression  with  BW 
AHR  (2.4-fold)  and  MO  (1.8-fold)  AHRs.  Two-way  ANOVA,  with  Bonferroni  post- 
hoc  tests;  p<0.05  (*);  p<0.001  (***), 
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Chapter  4:  Species-Specific  Differences  in  Mammalian  Aryl  Hydrocarbon 
Receptor  Ligand  Binding  Function 
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ABSTRACT 

The  aryl  hydrocarbon  receptor  (AHR)  is  a  ligand-activated  transcription  factor 
best  known  for  regulating  expression  of  genes  encoding  biotransformation 
enzymes  (e.g.  cytochrome  P450s)  and  other  genes  important  in  development, 
cell  cycle  control,  and  reproduction.  Subtle  species-specific  differences  in  AHR 
proteins  can  affect  the  cellular  response  to  chemical  exposure,  and 
subsequently,  organismal  sensitivity  to  environmental  contaminants.  Species- 
specific  differences  in  AHR  structure  and  function  were  compared  to  investigate 
variations  in  sensitivity  of  aquatic  and  marine  mammals  to  persistent  organic 
pollutants,  relative  to  other  model  mammalian  species.  Partial  AHR  cDNAs  were 
cloned  and  sequenced  from  mink  and  polar  bear,  and  were  compared  with  the 
ligand  binding  domains  of  other  marine  and  terrestrial  mammalian  AHRs  using 
amino  acid  alignment  and  phylogenetic  analysis.  The  relative  binding  affinities  of 
several  marine  and  terrestrial  mammal  AHRs  were  compared  in  saturation 
binding  assays  using  in  vitro  expressed  proteins.  Among  these  were  guinea  pig 
and  hamster  AHRs,  species  commonly  believed  to  have  high  and  low  affinity 
receptors  (respectively)  responsible  for  observed  differences  in  adult  dioxin 
lethality,  but  with  diverse  pathologies  resulting  from  sub-lethal  exposures. 
Homology  modeling  was  used  to  infer  the  tertiary  structures  of  AHRs  for 
representative  cetartiodactyl  (whale  +  cow),  carnivoran  (seal  +  dog),  and 
euarchontogliren  (human  +  guinea  pig  +  rodent)  AHR  ligand  binding  domains. 
While  the  models  illustrated  potential  structural  differences  that  might  influence 
the  AHR  ligand  binding  pocket,  the  major  differences  among  clades  did  not 
explain  differences  in  binding  observed  by  saturation  binding  assays.  Other 
regions  of  the  AHR  protein  likely  play  a  role  both  in  the  tertiary  conformation  of 
the  protein  and  its  overall  function. 
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1.  INTRODUCTION 


While  it  is  important  to  consider  the  ligand  binding  and  gene  activation 
functions  of  aryl  hydrocarbon  receptors  in  North  Atlantic  right  whales  for  reasons 
specific  to  the  species’  recovery,  placing  it  in  context  with  other  mammalian 
species  helps  in  the  interpretation  of  those  results.  The  study  of  marine 
mammal  toxicology  is  relatively  young,  whereas  model  organisms  such  as 
rodents  have  been  extensively  studied.  Integrating  these  lines  of  research  using 
a  comparative  approach  increases  the  power  of  individual  analyses,  and 
provides  a  more  in-depth  portrait  of  aryl  hydrocarbon  receptor  function  in 
mammals. 

Laboratory  studies  inform  our  understanding  of  species-specific  responses 
to  toxicant  exposure.  Rodent  studies  dominate  the  toxicology  literature, 
providing  basic  information  on  lethal  doses  and  phenotypes  associated  with 
controlled  toxicant  exposures.  Acute  toxicity,  wasting  syndrome,  endocrine 
imbalances,  immune  suppression,  lesions,  developmental  abnormalities,  and 
behavioral  changes  among  other  pathologies  have  been  observed  in  rats,  mice, 
guinea  pigs,  and  hamsters  following  exposure  to  2,3,7,8-tetrachlorodibenzo-p- 
dioxin  (TCDD,  dioxin)  and  dioxin-like  compounds  (DLC)  (Pohjanvirta  and 
Tuomisto  1994;  Poland  and  Knutson  1982).  More  limited  data  exist  for 
unconventional  and  wildlife  species  due  to  ethical  and  logistical  concerns 
governing  toxicological  experimentation,  and  anecdotal  evidence  or  accidental 
releases  often  prevail  in  this  arena.  Toxicological  experimentation  in  marine 
mammals  has  primarily  been  conducted  on  pinnipeds.  Several  experiments 
have  shown  that  harbor  seals  exposed  to  dietary  polychlorinated  biphenyls 
(PCBs)  suffer  reproductive  and  endocrine  effects  (Brouwer  et  al.  1989;  Reijnders 
2003;  Reijnders  1986).  Reproductive  dysfunction  and  morphological  disorders 
have  also  been  possibly  or  probably  associated  with  organochlorine  exposure  in 
other  pinnipeds  and  several  cetacean  species  (Reijnders  2003).  More  recently 
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mink  has  emerged  as  a  popular  subject  of  toxicological  studies  owing  to  its 
susceptibility  to  dioxin  toxicity,  closer  physiological  similarities  to  humans  than 
those  of  rodents,  and  ease  of  handling  (Calabrese  et  al.  1992;  Heaton  et  al. 

1995;  Hochstein  et  al.  1988);  such  studies  add  breadth  to  our  understanding  of 
phylogenetic  relationships  salient  to  toxicological  sensitivity. 

These  studies  show  the  importance  of  investigating  individual  species' 
sensitivities,  but  classical  toxicological  methods,  such  as  dosing  a  set  of 
individuals  and  measuring  specific  endpoints,  are  not  easily  transferable  to 
endangered  wildlife  or  large  animals.  There  is  a  need  for  indirect  measurements 
to  infer  susceptibility  to  environmental  toxicants.  Characterization  of  the 
structures  and  functions  of  different  species’  AHRs  can  help  to  guide  the  choice 
of  model  species  and  facilitates  the  extrapolation  of  results  from  one  species  to 
another. 

Halogenated  aromatic  hydrocarbons  (HAH),  polycyclic  aromatic 
hydrocarbons  (PAH),  and  DLC  can  exert  toxicological  effects  in  vertebrates  via 
the  aryl  hydrocarbon  receptor  (AHR).  The  AHR  is  a  cellular  transcription  factor 
and  a  member  of  the  PAS  family  of  proteins.  Sources  of  species-  and  tissue- 
specific  differences  in  organismal  sensitivity  to  environmental  toxicants  can 
include  ligand  pharmokinetics,  pharmocodynamics,  and  metabolism;  AHR 
functionality;  differences  in  AHR  cofactor  function,  expression  and  concentration; 
and  other  biochemical  and  physiological  factors.  However,  in  some  cases,  AHR 
binding  characteristics  can  be  directly  related  to  organismal  sensitivity. 

Two  issues  underlie  the  need  for  more  diverse  representation  of 
toxicological  sensitivity.  First,  model  organisms  (typically  rodents)  are  often 
studied  as  proxies  for  human  biomedical  applications.  While  this  approach  is 
sensible,  it  overlooks  species-specific  differences  that  can  be  critically  important 
in  predicting  the  effects  of  drugs  or  toxicants  in  humans.  It  is  vital  to  understand 
such  differences  in  order  to  appropriately  extrapolate  findings  from  one  organism 
to  another.  Second,  and  related,  is  the  application  of  the  same  principles  to 
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wildlife  health.  Although  often  discounted  as  unimportant  by  governing  bodies 
and  funding  agencies,  wildlife  health  is  important  for  maintaining  an  ecosystem 
able  to  support  human  populations,  and  is  sometimes  a  sentinel  for  human 
health  issues.  Understanding  the  relationships  of  model  or  surrogate  species  to 
d ifftcu It-to-study  wildlife  species  (especially  endangered  ones)  is  necessary  if  we 
are  to  apply  appropriate  and  timely  tools  for  conservation  and  resource 
management. 

AHR  acts  in  concert  with  many  other  cellular  proteins,  hormones,  and 
pathways.  Because  of  this,  it  is  not  always  possible  to  make  definitive 
connections  between  AHR  functional  characteristics  and  organismal  sensitivity  to 
DLC.  However,  studies  of  AHR  knock-out  mice  have  shown  that  binding  of  DLC 
to  AHR  is  the  necessary  proximal  step  for  most  chemical  effects  (Fernandez- 
Salguero  etal.  1996;  Mimura  et  al.  1997).  Yet  any  number  of  downstream 
factors  may  be  responsible  for  conferring  species-specific  susceptibility  to  DLC. 

It  remains  difficult  to  determine  whether  species-specific  variations  in  response  to 
toxicant  exposure  result  from  these  potential  downstream  differences,  or  from 
differences  in  AHR  function. 

In  this  chapter  we  use  a  comparative  approach  to  examine  AHR  ligand 
binding  domain  sequences  to  understand  species  that  exhibit  a  range  of 
sensitivities  to  dioxin-like  compounds  (DLC).  We  expand  the  range  of  species 
from  the  traditionally  studied  (mouse,  human,  guinea  pig,  hamster)  to  the  wild, 
marine  realm  (whales  and  seals).  Through  the  use  of  a  biochemical  assay, 
primary  protein  structural  analysis,  and  prediction  of  tertiary  protein  folding,  we 
attempt  to  infer  the  range  of  sensitivity  to  dioxin  along  a  continuum  of  low-  and 
high-affinity  AH  receptors.  These  analyses  are  admittedly  ambitious,  and 
inherently  fraught  with  errors;  however  we  hope  to  put  forth  a  model  that  may  be 
tested  with  more  rigorous  in  vitro  methods  or  field  data  in  the  future. 
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METHODS 

2.1  Aryl  Hydrocarbon  Receptor  Sequences 

In  order  to  examine  marine  mammal  AHRs  in  the  context  of  those  in  other 
mammals,  we  compiled  full  and  partial  AHR  amino  acid  sequences  from  this 
(Chapter  3)  and  other  studies  (Table  1 ).  In  addition,  several  of  these  AHR 
cDNAs  were  acquired  as  expression  plasmids,  kindly  provided  by  the 
researchers  who  developed  them:  mouse  (Burbach  et  al.  1992),  hamster 
(Korkalainen  et  al.  2000),  guinea  pig  (Korkalainen  et  al.  2001),  human  (Dolwick 
ef  al.  1993),  beluga  (Jensen  and  Hahn  2001),  and  harbor  seal  (Kim  and  Hahn 
2002).  Each  of  these  constructs  was  subcloned  into  the  pcDNA  vector  in  order 
to  standardize  the  conditions  under  which  each  protein  was  expressed  in  the  in 
vitro  system. 

2.2  RT-PCR  &  Sequencing  of  AHR  in  Mink  and  Polar  Bear 

Polyadenylated  RNA  was  isolated  and  purified  from  liver  tissues  from  mink 
and  polar  bear  samples.  Tissues  had  been  collected  and  stored  15  months 
(mink)  to  6  years  (polar  bear)  at  -80°C  prior  to  processing.  Complementary  DNA 
(cDNA)  was  synthesized  using  random/oligo  d(T)  primers  and  Omniscript  reverse 
transcriptase.  PCR  was  performed  with  degenerate  primers  D33F  and  D328R 
(Table  1)  with  the  following  thermocycling  protocol  with  AmpliTaq  Gold 
polymerase:  94°C  10  min;  35  cycles:  94°C  15  sec,  50°C  30  sec,  72°C  30  sec; 
72°C  7  min;  4°C  hold.  Bands  of  expected  size  were  excised  following  separation 
by  agaraose  gel  electrophoresis,  then  ligated  in  pGEM-T  Easy  cloning  vector. 
Positive  clones  were  again  size  selected  by  endonuclease  digestion  and  gel 
electrophoresis,  then  sequenced  using  Big  Dye  Terminator  cycle  sequencing 
with  isopropanol  purification.  Sequences  were  edited  using  Sequencher  4.2.2 
and  MacVector  8.1.1. 
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2.3  Sequence  Alignment  and  Analysis 

Aryl  hydrocarbon  receptor  amino  acid  sequences  from  Table  1 ,  plus  new 
partial  sequences  from  mink  and  polar  bear  AHRs,  were  aligned  using  the 
ClustalW  algorithm  in  MacVector  8.1.1 .  Phylogenetic  trees  were  constructed 
using  PAUP*  4.0. 

2.4  Biochemical  Characterization:  Ligand  Binding 

In  vitro  transcription  and  translation  (IVTT)  produced  AHR  proteins  from 
cDNAs  indicated  with  an  asterisk  (*)  in  Table  2.  The  amounts  and  sizes  of  each 
AHR  were  assessed  by  conducting  additional  IVTT  in  the  presence  of  [35S]- 
labeled  methionine,  separating  5  pi  of  each  reaction  by  SDS-PAGE,  performing 
fluorography  of  the  dried  gels,  then  cutting  and  counting  each  band  in  a 
scintillation  counter. 

Unlabeled  IVTT  AHRs  were  subjected  to  saturation  binding  analysis  with 
[3H]TCDD  as  described  in  Chapter  3.  Briefly,  proteins  were  diluted  in  buffer  and 
aliquoted  into  glass  vials  containing  one  of  8  concentrations  of  [3H]TCDD,  and 
incubated  overnight  at4°C.  Incubations  were  separated  into  duplicate  1.5  mL 
polypropylene  tubes  with  equal  volumes  of  buffer  plus  dextran  coated  charcoal 
(DCC).  Each  tube  was  then  vortexed  three  times  for  5  seconds  at  5  minute 
intervals.  Following  sedimentation  of  DCC  by  centrifugation,  supernatant  was 
removed  and  counted  for  presence  of  [3H]TCDD.  Specific  binding  was 
determined  as  detailed  in  Chapter  3,  using  a  series  of  tubes  containing  UPL 
(unprogrammed  lysate)  to  account  for  non-specific  binding.  Binding  curves  and 
constants  were  calculated  using  Prism  4.0. 

2.5  Homology  Modeling 

AHR  LBD  amino  acid  sequences  for  the  species  listed  in  Table  2  were  submitted 
directly  to  the  SWISS-MODEL  database  employing  the  First  Approach  method 
(http://www.exDasv.ch/swissmod/SWISS-MODEL.html).  HIF2«  was  manually 
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selected  as  the  template,  but  was  also  chosen  by  the  automated  SWISS-MODEL 
server  as  the  best  template.  Models  were  viewed  using  the  free  software  Swiss- 
Pdbviewer  for  Mac  OSX  (http://www.expasv.ch/SDdbv/mainpaae.html)  (Guex  and 
Peitsch  1997;  Peitsch  1995;  Schwede  etal.  2003).  RMS  deviations  were 
calculated  before  and  after  each  three-dimensional  model  was  adjusted  with  Best 
Structural  Fit  relative  to  the  mouse  AHR  model  to  assess  improvements  in 
alignment.  Previous  studies  have  shown  the  PAS  structures  are  useful  for 
predicting  functional  site  locations  and  general  features  of  unknown  mechanisms 
(Denison  etal.  2002;  Pandini  and  Bonati  2005;  Procopio  etal.  2002). 

3.  RESULTS 

3. 1  Mink  and  Polar  Bear  AHR 

Unique  nucleotide  sequences  were  obtained  from  mink  and  polar  bear 
RNAs  (Figure  1),  sharing  high  identity  with  other  AHR  sequences  published  in 
GenBank  in  BLASTn  searches  (mink:  93%  identity  to  Carnivore  and  88%  identity 
to  Cetacean  AHRs;  polar  bear:  97%  identity  to  Primate,  and  91-95%  identity  to 
Cetacean  and  Carnivore  AHRs).  BLASTp  searches  of  GenBank  revealed  94- 
98%  amino  acid  identity  to  known  AHRs  included  in  this  study.  The  top  three 
database  matches  for  mink  were  dog,  harbor  seal,  and  Baikal  seal  AHRs;  for 
polar  bear  they  were  white-beaked  dolphin,  Atlantic  white-sided  dolphin,  and  pilot 
whale  AHRs,  These  sequences  correspond  to  a  highly  conserved  portion  of  the 
AHR  protein,  the  PAS  region,  and  degenerate  primers  have  been  successful  in 
amplifying  this  region  from  many  vertebrate  species  including  cetaceans 
(Chapter  2),  fish,  a  turtle,  and  amphibians  (Karchner  et  al.  2000;  Karchner  et  at. 
1999;  Lavine  et  al.  2005).  Repeated  attempts  at  extending  the  sequence 
information  by  rapid  amplification  of  cDNA  ends  were  unsuccessful  for  both  mink 
and  polar  bear  (data  not  shown).  The  presence  of  introns  impeded  our  ability  to 
obtain  additional  coding  sequence,  and  may  have  been  the  result  of  incompletely 
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spliced  RNA.  Genomic  DNA  contamination  in  the  RNA  solutions  may  have 
reduced  amplification  efficiency.  Contamination  with  genomic  DNA  was 
combated  with  DNase  treatment  and  further  mRNA  purification,  but  these 
procedures  did  not  yield  better  results. 

3.2  Mammalian  AHR  Amino  Acid  Sequence  Alignment 

Functional  domains  have  been  determined  experimentally  for  mouse  and 
human  AHRs.  The  ligand  binding  domain  (LBD)  roughly  corresponds  to  the 
PAS-B  domain,  and  includes  amino  acids  230-402  in  mouse  and  236-408  in 
human.  Alignment  of  full-length  AHR  sequences  listed  in  Table  2  illustrates  the 
high  conservation  of  amino  acid  sequences  among  mammals  (Figure  2),  and 
therefore  we  used  the  mouse  and  human  AHRs  as  templates  for  identifying 
functional  domains  in  new  sequences.  The  variability  in  the  carboxy-terminus 
noted  by  other  researchers  is  less  pronounced  in  this  alignment,  probably  due  to 
the  increased  taxa  sampling  within  the  eutherian  mammals.  Ligand  binding 
domains  from  full-length  and  partial  AHR  sequences  were  aligned  using  the 
functionally  determined  mouse  and  human  LBDs  as  guides  (Figure  3).  The  high 
degree  of  amino  acid  identity  among  LBDs  (84  -  99%  within  mammals,  76-82% 
to  killifish)  suggests  conservation  of  function  among  vertebrate  AHRs. 

Maximum  parsimony  analysis  of  the  LBD  alignment  identified  22 
parsimony  informative  characters  among  the  173  amino  acid,  ungapped  data  set 
(indicated  in  Figure  2B  with  arrows  and  green  shading).  The  most  parsimonious 
phylogenetic  tree  was  determined  with  branch  and  bound  analysis.  The  resulting 
tree,  rooted  with  killifish  AHR1,  illustrates  the  number  of  similarities  among  these 
species  (Figure  4A).  In  particular,  the  cetacean  branches  are  completely 
unresolved  and  essentially  identical.  Similar  analysis  of  full-length  AHR  amino 
acid  sequences  from  these  same  species  introduces  more  areas  of  variability 
(229  parsimony  informative  characters  of  769  residues,  excluding  gaps), 
permitting  differentiation  of  each  taxon.  The  tree  derived  from  full-length 
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sequence  more  closely  approximates  the  accepted  mammalian  phylogeny 
(Figure  4B;  (Murphy  et  al.  2001 )).  These  results  indicate  that  the  LBD  contains 
limited  information  for  determining  phylogenetic  relationships  of  each  AHR,  and 
may  be  more  suitable  for  identifying  similarities  among  species. 

3.3  Biochemical  Characterization 

At  this  point,  assessments  of  ligand  binding  activity  cannot  rely  solely  on 
sequence  similarities  to  predict  the  binding  characteristics  of  AHRs  from 
unstudied  species.  Experimental  data  are  required  to  define  ligand  binding 
affinities,  and  ideally  the  experimental  design  should  include  all  species  of 
interest  within  the  same  assay  for  comparison.  In  the  previous  chapter, 
extensive  comparisons  were  made  between  North  Atlantic  right  whale, 
humpback  whale,  beluga  whale,  mouse,  and  human  AHRs  to  determine  specific 
binding  and  binding  constants.  The  approach  here  is  similar,  but  includes  more 
species,  on  a  smaller  scale. 

Two  independent  saturation  binding  assays  were  conducted,  the  first 
using  all  8  concentrations  of  [3H]TCDD,  the  other  with  a  single  high  concentration 
(5,5  nM,  nominal,  2-3  nM  final)  plus  a  DMSO  control.  The  same  volume  of  diluted 
IVTT  AHR  (100  pL)  was  used  for  each  receptor  in  the  saturation  binding  assay. 
The  humpback  whale  and  hamster  AHR  constructs  using  pcDNA  vector  were 
only  available  for  the  latter  experiment.  Specific  binding  was  calculated  as  the 
difference  between  total  binding  and  non-specific  binding  (UPL).  Nonlinear 
regression  analysis,  equilibrium  dissociation  constants  (Kd),  and  maximum 
binding  (Bmax)  were  derived  for  the  first  experiment  assuming  one-site  specific 
binding  (Figure  5). 

A  single  concentration  of  [3H]TCDD  was  incubated  with  IVTT  expressed 
AHR  proteins  to  determine  the  maximum  specific  binding  capacity  of  each 
species’  AHR.  Whale  and  hamster  AHRs  exhibited  the  greatest  binding 
capacities  (>0.9  nM;  Figure  6),  but  also  produced  higher  amount  of  protein  than 
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the  other  constructs.  In  general,  greater  concentrations  of  IVTT  AHR  resulted  in 
higher  specific  binding  when  incubated  with  a  single  high  concentration  of 
[3H]TCDD  (r2  =  0.59;  Figure  7).  However,  even  in  the  same  expression  vector, 
each  construct  does  not  express  the  same  amount  of  AHR  protein.  We  suggest 
that  future  experiments  account  for  the  difference  in  AHR  expression  by 
modifying  the  proportion  of  specific  protein. 

Previous  work  (e.g.  Chapter  3)  has  shown  variability  among  replicate 
saturation  binding  assays,  complicating  the  statistical  differentiation  of  individual 
AHRs  within  a  narrow  range  of  relative  binding  affinity  (Kd)  values.  In  order  to 
place  the  single  KD  value  obtained  in  this  study  in  context  with  such  variation,  we 
included  the  relative  binding  affinities  calculated  from  multiple  saturation  binding 
assays,  including  the  one  represented  in  Figure  4;  these  data  are  shown  in 
Figure  8.  The  Kd  for  the  harbor  seal  AHR  is  similar  to  that  obtained  for  the  initial 
characterization  of  that  protein  (Kim  and  Hahn  2002).  While  some  receptors 
exhibit  a  greater  range  of  values  than  others  in  these  assays,  the  distributions  of 
Kds  are  balanced  above  and  below  the  means.  The  single  values  for  harbor  seal 
and  guinea  pig  AHRs  are  within  the  expected  range  of  values,  but  may  not  fully 
represent  the  range  of  binding  affinities  specific  to  each  species. 

3. 4  Homology  Modeling 

Tertiary  structures  of  the  ligand  binding  domain  were  predicted  using  the 
SwissMode!  PDB  First  Approach  method  for  right  whale,  humpback  whale, 
beluga  whale,  cow,  dog,  harbor  seat,  Baikal  seal,  mouse,  hamster,  and  guinea 
pig  aryl  hydrocarbon  receptors.  All  structures  were  based  on  the  HIF2«  (PDB# 
1p97a)  peptide  characterized  by  NMR  (Erbel  etal.  2003).  Fragments  from  mink, 
polar  bear,  and  other  cetaceans  (sequences  from  Chapter  2)  were  too  small  for 
modeling.  Computations  for  energy  minimization  of  each  model  were 
accomplished  with  GROMOS96  implementation,  and  each  structure  was  aligned 
relative  to  the  mouse  AHR  model  by  Best  Structural  Fit  in  the  Swiss-PDBViewer. 
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Similarity  of  tertiary  structure  was  the  overwhelming  result  of  these  first 
approximations.  All  models  possess  a  large  six-stranded,  anti-parallel  beta-sheet 
formed  by  a  combination  of  amino-  and  carboxy-ends  of  the  LBD.  One  major 
alpha-helix  occurs  at  the  approximate  mid-point  in  the  LBD,  and  two  other 
smaller  helices  (or  a  single  “broken"  or  twisted  helix)  flank  the  beta-sheet.  An 
irregular  loop  connects  these  helices  and  sometimes  includes  a  beta-turn  or 
partial  helix,  depending  on  the  species.  Together,  these  structures  define  a 
hollow  cylindrical  space  that  is  the  likely  ligand  binding  “pocket”  (Figure  9,  Figure 
10A). 

The  morphology  of  the  predicted  tertiary  structure  of  AHR  LBDs  differs  at 
two  irregularly  structured  regions.  A  Minor  Loop  within  the  beta-sheet  shows 
slight  variation  among  some  species  (e.g.  it  is  the  site  of  the  only  apparent 
differences  in  human  versus  rodent  models),  and  may  be  dependent  on  the 
overall  beta-sheet  construction  specific  to  each  sequence.  A  second  irregular 
feature  (Major  Loop)  varies  from  an  undifferentiated  loop  to  a  loop  with  a  beta- 
turn.  This  feature  spans  the  region  between  the  major  alpha-helix  and  the  small 
helices,  occupies  a  larger  relative  space  than  the  other  variable  loop  in  the  beta- 
sheet,  and  as  such  may  be  capable  of  greater  flexibility.  The  latter  site  is  found  to 
be  the  site  of  the  greatest  variation  among  models. 

The  predicted  tertiary  structures  of  all  three  whale  AHR  LBDs  were 
identical,  as  expected  from  the  paucity  of  primary  sequence  variation  in  the  LBD 
seen  in  the  alignment  and  accentuated  in  the  phylogenetic  tree.  Cow,  dog,  and 
harbor  seal  also  share  identical  tertiary  structures  with  the  whale  LBD,  but  the 
Baikal  seal  differs  noticeably  in  one  location.  Among  the  Cetartiodactyls  and 
Carnivora,  the  Major  Loop  contains  a  more  rigid  predicted  structure  that  be  may 
considered  a  beta-turn,  with  the  exception  of  the  Baikal  seal,  which  lacks  this 
feature  and  remains  an  irregular  and  highly  flexible  loop  (Figure  10B).  The 
functional  relevance  of  this  potential  mobility  cannot  be  inferred  in  the  absence  of 
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in  vitro  AHR  binding  data  or  in  vivo  studies  for  the  Baikal  seal.  This  Major  Loop 
structure  is  otherwise  unique  to  Cetartiodactyla  and  Carnivora. 

Interestingly,  the  mouse  and  hamster  AHR  LBDs  are  identical  to  one 
another,  and  both  differ  from  the  human  model  only  slightly  at  the  Minor  Loop 
(Figure  IOC).  Although  large  differences  in  KD  were  shown  experimentally,  the 
Bmax  values  for  these  species  were  very  similar  and  this  may  be  related  to  the 
similar  structure  of  the  ligand  binding  pockets. 

The  model  for  guinea  pig  AHR  LBD  is  decidedly  different  (Figure  9D).  A 
much  larger  portion  of  the  protein  lacks  tertiary  structure,  including  the  loss  of 
one  beta  strand  and  a  truncated  second  beta-strand  at  the  amino-end  of  the 
sequence,  an  abbreviated  minor  alpha-helix,  and  the  absence  of  the  Major  Loop 
beta-turn  (the  latter  as  in  the  other  Euarchontoglires).  The  major  alpha-helix  is 
more  loosely  coiled.  The  end  result  is  a  very  unstructured  LBD  that  may  be  more 
susceptible  to  the  influences  of  other  parts  of  the  full-length  protein,  and  may 
perhaps  have  consequences  for  ligand  binding. 


4.  DISCUSSION 

While  others  have  shown  that  certain  differences  in  the  amino  acid 
sequence  of  AHR  result  in  changes  in  ligand  binding  affinity,  none  have 
questioned  whether  these  differences  can  be  used  to  infer  a  species-level 
predilection  to  POP  sensitivity  and  increased  risk  of  biological  effects.  This 

approach  could  become  a  first-approximation  tool  for  regulating  human  activities 
and  improving  environmental  conditions  for  endangered  species. 

4. 1  New  AHR  Sequence  Acquisition 

As  seen  in  the  alignment  of  full-length  AHR  sequences,  an  expanded 
series  of  closely  related  taxa  improves  the  overall  alignment.  Specifically,  the 
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carboxy-terminal  region  of  the  protein  appears  much  less  variable  than  other 
authors  have  observed  when  comparing  more  distantly  related  species.  While 
our  knowledge  of  full-length  AHR  sequences  is  still  limited,  the  partial  sequences 
recently  acquired  for  odontocete  cetaceans  (Chapter  2),  mink  and  polar  bear  do 
provide  more  information  on  the  diversity  of  AHRs  within  Mammalia.  Subtle 
differences  in  binding  characteristics  and  inferred  tertiary  structure  may  not  be 
fully  explained  by  these  partial  sequences;  therefore,  continued  efforts  to  clone 
and  sequence  full-length  AHR  cDNAs  remains  worthwhile  to  better  understand 
this  pathway. 

4.2  Relative  Binding  Affinities  of  Mammalian  AHRs 

Although  individual  species  are  often  studied  in  depth,  few  studies  have 
directly  compared  the  binding  characteristics  of  multiple  species  in  order  to 
understand  the  relative  differences  in  organismal  sensitivity  to  DLC.  Here  we 
show  the  relative  binding  affinities  of  AHRs  from  both  endangered  wildlife 
species  and  model  organisms,  but  replicate  assays  are  necessary  to  increase 
confidence  in  the  mean  binding  constants  and  account  for  assay  variability.  A 
single  experiment  was  conducted  using  the  hamster  AHR  cDNA  in  an  alternative 
expression  vector;  however,  the  results  were  ambiguous  (data  not  shown).  We 
have  found  that  consistent  use  a  single  expression  vector  reduces  intra-assay 
variability  and  eliminates  a  factor  of  uncertainty;  future  experiments  should 
include  the  use  of  hamster  AHR  in  pcDNA. 

4.3  Predicting  Sensitivity  by  Amino  Acid  Analysis:  functionally  significant 
polymorphisms  and  critical  residues 

The  significance  of  any  one  amino  acid  in  AHR  function  has  been  explored 
primarily  in  rodent  and  human  studies  (Okey  et  al.  2005).  Eight  point  mutations 
are  known  among  four  strains  of  mice,  two  of  which  occur  in  the  PAS  regions 
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causing  single  amino  acid  changes.  The  switch  from  alanine  to  valine  at  position 
375  (A375V)  results  in  a  reduction  in  affinity  for  ligand;  accompanied  by  a  L348F 
change,  binding  is  reduced  to  a  “dioxin-resistant"  phenotype  (Poland  etal.  1994). 
A  similar  situation  was  found  in  the  comparison  of  the  high-affinity  mouse  AHR 
with  low-affinity  human  AHR.  The  mouse  amino  acid  375  (A)  corresponds  to 
residue  381  (V)  in  human,  and  a  V381D  mutation  in  the  human  AHR  abolishes 
ligand  binding  (Ema  et  at.  1994;  Ramadoss  and  Perdew  2004). 

Han/Wistar  rats  are  among  the  most  dioxin-resistant  mammals  tested,  and 
possess  a  truncated  AHR  protein  resulting  from  a  single  point  mutation  at  the 
start  of  intron  10,  which  disrupts  proper  splicing  and  eliminates  approximately  40 
amino  acids  from  the  carboxy-terminus  (Pohjanvirta  et  at.  1998).  While  the  rat 
AHR  mutation  is  outside  the  range  of  the  LBD,  it  illustrates  the  power  of  a  single 
site  to  influence  the  structure  of  the  protein.  Hamsters  are  similar  to  Han/Wistar 
rats  in  exhibiting  a  more  dioxin-resistant  phenotype,  whereas  another  popular 
rodent  model,  the  guinea  pig,  is  among  the  most  dioxin-sensitive.  In  hamsters, 
the  glutamine  (Q)  rich  carboxy-terminus  is  substantially  expanded  and  enriched 
with  Q  compared  with  other  AHRs  (Korkalainen  et  at.  2000).  In  contrast,  the 
same  region  in  the  guinea  pig  AHR  is  only  half  the  size  of  that  from  hamster.  A 
distinct  correlation  can  be  seen  across  published  mammalian  AHR  sequences 
between  the  number  of  glutamine  residues  in  the  carboxy-terminal  region  and 
sensitivity  to  acute  dioxin  lethality  (Korkalainen  et  at.  2001).  These  studies 
illustrate  the  potential  importance  of  regions  besides  the  LBD  in  the  function  of 
AHR  in  mammals. 

A  recent  study  in  birds,  however,  addresses  the  specific  importance  of 
single  amino  acid  differences  the  AHR  LBD  (Karchner  et  at.  2006).  Chicken  and 
common  tern  are  closely  related  species  that  have  highly  similar  AHR  amino  acid 
sequences,  yet  terns  are  80-  to  250-fold  less  sensitive  than  chickens  to  the 
effects  of  HAH  (Brunstrom  and  Halldin  1998;  Hoffman  et  at.  1998;  Lorenzen  et  at. 
1 997).  Only  three  amino  acids  differ  in  the  LBD  of  these  AHRs.  Site-directed 
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mutagenesis  experiments  show  that  two  of  these  sites  in  the  chicken  AHR  (1324 
and  S380)  are  responsible  for  its  high-affinity  dioxin  binding,  as  well  as  its  ability 
to  activate  gene  expression  in  a  reporter  gene  system  (Karchner  et  al.  2006).  In 
contrast,  the  common  tern  AHR  LBD  possesses  valine  and  alanine  residues  at 
the  analogous  sites  and  has  reduced  binding  and  transactivation  functions.  This 
study  succinctly  illustrates  the  relationship  between  a  few  amino  acids  in  the 
AHR  LBD  and  dramatic  differences  in  organismal  sensitivity  to  dioxin. 
Additionally,  the  chimeric  chicken-tern  AHR  constructs  used  in  this  study  show 
that  not  only  differences  in  the  carboxy-terminal  "transactivation  domain"  impact 
transcriptional  activity  but  that  the  LBD  may  influence  this  function  as  well. 

4.4  Homology  Modeling  AHR 

Homology  models  are  only  as  good  as  the  original  structural  data  used  as  a 
template.  Because  neither  an  NMR  structure  nor  an  x-ray  crystal  structure  exists 
for  any  species'  AHR,  we  based  our  models  on  a  related  PAS-family  protein, 
hypoxia-inducible  factor  2-alpha  (HIF2a;  PDB  #  1 P97).  Although  PAS-family 
proteins  do  not  share  a  large  amount  of  direct  amino  acid  sequence  identity, 
certain  regions  -  namely  the  PAS  regions  -  contain  regular  motifs  suitable  for 
inferring  structural  elements.  Appropriate  templates  for  homology  models  need 
only  share  only  20-25%  amino  acid  sequence  identity  with  the  target  sequence 
(Gale  Rhodes,  website,  unreferenced).  The  structure  for  the  PAS-B  region  of 
HIF2a  is  based  on  NMR  spectroscopy,  and  is  the  synthesis  of  26  models  (Erbe! 
et  al.  2003).  The  HIF2a  PAS-B  model  may  be  the  best  template  available  until 
crystal  structures  for  an  AHR  become  available.  In  the  meantime,  attempts  at 
homology  modeling  AHRs  are  only  approximations  of  reality,  and  serve  best  as 
illustrations  of  hypothetical  conditions,  or  for  guiding  experimental  design. 

Our  goal  was  to  compare  the  relative  structures  of  AHRs  from  different 
mammalian  species.  In  pursuit  of  this  goal,  whether  or  not  the  model  accurately 
represents  the  true  structure  of  AHR  is  less  relevant.  As  long  as  all  the  AHR 
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models  here  were  constructed  in  like  fashion  (as  they  indeed  were),  comparisons 
among  and  between  species  should  remain  valid  approximations  of  species- 
specific  variations.  Combined  with  the  experimental  work  of  many  others,  these 
comparisons  contribute  to  the  hypothesis  that  discrete  residues  in  the  LBD  are 
critical  for  AHR  function.  Refinement  of  homology  modeling  techniques  or  the 
addition  of  AHR  crystal  structures  to  existing  databases  could  advance  the 
accuracy  of  our  comparative  models. 

Others  have  recently  attempted  to  use  homology  modeling  techniques  to 
infer  the  tertiary  structure  of  the  aryl  hydrocarbon  receptor  based  on  the  crystal 
structures  of  HERG  and  FixL  (Denison  ef  a/.  2002;  Procopio  et  al.  2002).  The 
AHR  models  presented  here  are  consistent  with,  and  only  subtly  differ  from, 
these  models.  Ours  is  the  only  AHR  structure  modeled  after  HIF2a,  and  may 
represent  a  more  accurate  representation  as  both  the  template  and  AHR 
dimerize  with  ARNT  via  the  beta-sheet  in  the  LBD.  We  expect  that  the  remainder 
of  the  AHR  protein  influences  the  final  folding  conformation  and  consequently  the 
function;  thus,  there  remains  a  need  for  additional  work  including  crystallization  of 
the  AHR.  While  our  results  are  approximations,  this  data  set  increases  the 
availability  of  sequence  information  and  increases  confidence  in  previous 
homology  models  for  AHR  by  approaching  the  model  via  a  new  template. 

The  combination  of  sequence  analysis,  biochemical  characterization,  and 
homology  modeling  presents  another  tool  for  explaining  differences  in  AHRs. 

We  suggest  that  this  approach  provides  a  good  framework  for  presenting 
testable  hypotheses  using  increasingly  available  sequence  data  in  public 
databases,  and  could  be  extended  for  use  in  other  systems.  Several  questions 
have  emerged  from  this  work  thus  far  for  which  we  could  imagine  developing 
experiments.  First,  for  species  that  show  high  conservation  of  amino  acid 
sequence  and  whose  tertiary  structural  models  indicate  identical  LBDs, 
differences  in  ligand  binding  may  not  be  related  to  the  PAS-B  region.  What  is  the 
impact  of  the  few,  specifically  different  residues  on  ligand  binding?  Could  these 
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variable  sites  confer  differences  in  ligand  specificity  that  are  not  seen  in 
experiments  with  TCDD?  Alternatively,  other  portions  of  the  protein  may 
influence  the  accessibility  of  the  binding  pocket.  Furthermore,  does  the  size  and 
accessibility  of  the  binding  pocket  alter  the  probability  of  ligand  insertion?  Could 
a  wider  potential  binding  pocket  relax  constraints  on  ligand  docking,  and  if  so, 
would  there  be  consequences  for  binding  efficacy?  These  questions  could  be 
addressed  first  with  virtual  mutagenesis  and  ligand  docking  studies  using  the 
existing  models,  followed  by  site-directed  mutagenesis  and  chimeric  or  truncated 
AHR  constructs  in  functional  bioassays. 

Acknowledgements:  R.J.  Aulerich  kindly  provided  the  mink  tissue  used  for 
partially  cloning  and  sequencing  the  AHR  cDNA. 
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Table  1.  Degenerate  aryl  hydrocarbon  receptor  primers. 

Name  Sequence  (5’  ->  3’)  Tm  (°C) 

D33F  GTC  IAA  CCC  ITC  IAA  GMG  ICA  YMG  64 

D143F  CAR  GAY  TAY  CTI  GGI  TTY  CAI  CA  62 

D328R  GCI  GCG  TGR  ATR  AAY  TGR  TAI  CC  63 

D348R  CAT  ICC  RCT  YTC  ICC  IGT  YTT  62 

Designed  by  Karchner/Hahn;  synthesized  by  Integrated  DNA  Technologies. 
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Table  2.  Aryl  hydrocarbon  receptor  sequences  used  for  multiple  alignment 
and  homology  modeling.  Asterisks  indicate  those  sequences  used  for  in  vitro 
assays. 


Species 


GenBank  Accession  # 


North  Atlantic  right  whale  (NARW;  Eubalaena  glacialis) 
Humpback  whale  { Megaptera  novaeangliae) 

Beluga  whale  ( Delphinapterus  leucas) 

Pilot  whale  ( Globicephalus  melas) 

Atlantic  white-sided  dolphin  (AWSD;  Lagenorhynchus  acutus) 
White-beaked  dolphin  (WBD;  L  albirostris) 

Cow  ( Bos  taurus ) 

Harbor  seal  ( Phoca  vitulina) 

Baikal  seal  ( Phoca  siberica) 

Mink  (Mustela  vison)  this 

Polar  Bear  ( Ursus  maritimus)  this 

Dog  { Canis  familiaris) 

Human  ( Homo  sapiens) 

Mouse  ( Mus  musculus) 

Guinea  Pig  ( Cavia  porcellus) 

Hamster  (Mesocricetus  auratus) 

Killifish  AHR1  ( Fundulus  heteroclitus) 


EF091829  * 
EF091830  * 
AF332999 * 
EF091833 
EF091832 
EF091831 

XP612996;  predicted 
AB056700 * 
AB072432 

study;  not  yet  submitted 

study;  not  yet  submitted 
XP532485;  predicted 
LI  9872  * 

M 94623  * 

AY028947 * 
AF275721  * 
AF024591  * 


HIF2ri  (PDB  #1 P97) 
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Figure  1.  Partial  eDNA  sequences  of  mink  and  polar  bear  aryl  hydrocarbon 
receptor. 

Sequences  derived  from  nested  RT-PCR  yielded  partial  aryl  hydrocarbon 
receptor  cDNAs  (from  polyadenylated  RNA,  primers  D33F/D348R  followed  by 
D143F/D328R). 
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Figure  2.  Alignment  of  full-length  mammalian  AHR  amino  acid  sequences. 

Aryl  hydrocarbon  receptor  amino  acid  sequences  for  which  the  entire  cDNA  is 
known  from  cloning  or  genome  analysis  were  aligned  using  ClustalW  in 
MacVector  8.1.1.  Killifish  ( Fundulus  heteroclitus)  AHR1  is  shown  in  the  bottom 
row.  Identical  residues  are  outlined  and  shaded  grey.  The  full  alignment  spans 
two  pages. 
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Figure  2  continued. 
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Figure  3.  Aryl  hydrocarbon  receptor  ligand  binding  domains  are  well 
conserved  in  mammals. 

(A.)  A  general  schematic  of  a  prototypical  AHR  is  shown;  lower  scale  bar 
indicates  cDNA  length  in  kilobases,  and  boxes  represent  protein  sequence  with 
the  basic-helix-loop-helix  (bHLH),  PAS,  and  glutamine-rich  (Q-rich)  regions 
shaded.  Based  on  functional  studies  of  mouse  and  human  AHRs,  the  ligand 
binding  domain  (LBD;  boxed  in  A)  encompasses  the  PAS  B  domain;  homologous 
regions  of  other  mammalian  AHRs  were  selected  for  alignment.  Amino  acid 
sequences  of  16  mammalian  and  one  fish  AHR  LBD  were  aligned  using 
ClustalW  in  MacVector  8.1.1  (B).  Identical  residues  are  outlined  and  shaded 
grey.  The  high  degree  of  amino  acid  similarity  suggests  conservation  of  function 
among  vertebrate  AHRs.  Parsimony  informative  characters  identified  using 
PAUP*4.0  are  indicated  with  arrows  and  green  shading.  A  single  amino  acid 
responsible  for  much  of  the  binding  difference  in  human  and  mouse  AHRs  is 
indicated  with  the  upward-pointing,  dashed  arrow.  Irregular  structural  regions  or 
loops  showing  the  greatest  variation  in  homology  models  are  identified  in  aqua 
and  orange  boxes,  corresponding  to  the  sequential  secondary  structure  coloring 
in  Figure  9. 
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Figure  4.  Phylogenetic  analysis  of  mammalian  AHR. 

Maximum  parsimony  analysis  of  either  the  LBD  (A)  or  the  full  length  AHRs  (B) 
generally  recapitulated  the  accepted  view  of  mammalian  phylogeny  (Murphy  et 
al.  2001).  Analysis  of  just  the  LBD  identified  22  informative  residues  of  the  173 
ungapped  amino  acid  characters  (as  highlighted  in  Figure  3),  but  was  unable  to 
differentiate  between  the  whale  AHRs,  as  shown  by  the  collapsed  branches.  The 
LBD  alone  is  not  sufficient  to  decipher  phylogenetic  relationships.  Outgroup  was 
killifish  ( Fundulus  heteroclitus)  AHR1. 


144 


A 


harbor  seal 

Baikal  seal 

dog 

cow 

right  whale 
humpback  whale 
beluga 

human 
guinea  pig 

mouse 

hamster 


CARNIVORA 


CETARTIODACTYLA 


PRIMATES 


RODENTIA 


kllltflsh 


right  whale 
humpback  whale 
beluga  whale 

cow 

harbor  seal 

Baikal  seal 

dog 

human 

guinea  pig 

hamster 

mouse 

killtflati 


CETARTtODACT 


CARNIVORA 


PRIMATES 


RODENTIA 


145 


EUARCHONTOGURES  >  EUARCHONTOGLIRES 


0.55n 


■  RW 
*  BW 


— I - 

1 

2 

3 

4 

5 

Free 

3H-TCDD  (nM) 

RW 

BW 

GP 

HS 

MO 

HU 

BMAX 

0.5709 

0,3199 

0  3748 

0.2291 

0.1414 

0  1527 

KD 

08028 

0.6292 

08222 

0,6337 

0,2998 

0  8692 

Figure  5.  Saturation  binding  analysis  of  IVTT  AHRs  from  marine  mammals 
and  rodents. 

Specific  binding  curves  were  derived  from  incubation  of  IVTT  AHRs  with  nine 
concentrations  of  [3H]TCDD  in  a  saturation  binding  assay.  Maximal  binding 
(Bmax)  and  relative  binding  affinities  (Kd)  specific  to  this  experiment  are  shown 
below. 
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Figure  6.  Maximum  binding  capacity  of  mammalian  AHRs  in  vitro. 

A  single  concentration  of  [3H]TCDD  was  incubated  with  IVTT  expressed  AHR 
proteins  to  determine  the  maximum  specific  binding  capacity  of  each  species’ 
AHR.  Whale  and  hamster  AHRs  exhibited  the  greatest  binding  capacities  (>0.9 
nM),  but  also  produced  higher  amount  of  protein  than  the  other  constructs. 
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Figure  7.  Specific  binding  is  related  to  AHR  protein  concentration. 

In  general,  greater  concentrations  of  IVTT  AHR  resulted  in  higher  specific  binding 
when  incubated  with  a  single  high  concentration  of  [3H]TCDD  (5.5  nM  nominal;  r 
=  0.59).  The  same  volume  of  IVTT  AHR  (100  pL)  was  used  for  each  receptor  in 
the  maximum  binding  capacity  assay.  However,  even  in  the  same  expression 
vector,  each  construct  does  not  express  the  same  amount  of  AHR  protein. 
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Figure  8.  Relative  binding  affinities  for  IVTT  expressed  AHRs  in  saturation 
binding  assays. 

Relative  binding  affinities  of  multiple  saturation  binding  assays  with  [3H]TCDD  are 
combined  in  a  boxplot  to  illustrate  the  range  of  values  acquired.  Sample  sizes 
vary:  MO  n  =  6,  RW  n  =  9,  HS  n  =  1 ,  BW  n  =  6,  HW  n  =  5,  GP  n  =  1 ,  HU  n  =  6. 
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Minor  Loop 


Figure  9.  Predicted  tertiary  structure  of  AHR  LBD  as  represented  by  the 
NARW  homology  model. 

The  basic  morphology  of  the  AHR  LBD  is  illustrated  by  with  the  homology  model 
of  the  north  Atlantic  right  whale  sequence.  Sequential  secondary  structure 
coloring  begins  at  the  amino  terminal  end  (N)  in  blue,  and  concludes  at  the 
carboxy  terminal  end  of  the  LBD  (C)  in  red.  All  AHR  LBDs  modeled  after  the 
HIF2a  PAS  domain  share  major  structural  features  including  a  six-stranded  beta 
sheet,  a  highly  conserved  structure  among  PAS  family  proteins  important  for 
cofactor  binding.  The  beta  sheet  is  flanked  by  a  large  alpha  helix  on  one  side 
(green),  and  several  smaller  helices  on  the  other  (aqua  shades).  Pairwise 
comparison  of  homology  models  revealed  two  regions  of  variability:  a  Minor  Loop 
between  beta  strands,  and  a  Major  Loop  between  the  helices. 
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Figure  10.  Comparisons  of  AHR  LBD  predicted  tertiary  structure  among 
mammals  by  homology  modeling. 

An  alternative  perspective  (A,  top  view)  of  the  NARW  AHR  model  suggests  the 
space  in  which  a  ligand  might  fit.  With  the  exception  of  the  Baikal  seal,  ail 
Cetartiodactyl  and  Carnivore  models  were  identical,  and  possess  a  small  beta- 
turn  or  partial  helix  in  the  midst  of  the  Major  Loop  (B).  In  contrast,  human  and 
rodent  models  had  no  such  structural  element  in  the  Major  Loop,  and  were 
virtually  identical  except  for  subtle  deviations  in  the  Minor  Loop  <C).  Additional 
differences  are  seen  in  the  guinea  pig  model  (D),  such  as  truncated  beta  strands, 
and  more  undifferentiated  loops,  These  models  cannot  resolve  the  differences 
observed  in  binding  properties,  such  as  those  seen  in  mouse  and  human,  but 
offer  an  interesting  perspective  on  the  relative  importance  of  individual  amino 
acid  differences  on  receptor  activity,  and  help  to  generate  many  more  questions 
about  the  relationship  between  AHR  structure  and  function. 
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Chapter  5:  Summary  and  Future  Research 
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Extensive  environmental  contamination  combined  with  the  inherent 
complexity  of  the  natural  world  together  present  a  difficult  problem  to  people 
concerned  with  regulating  and  reducing  harmful  inputs,  and  for  mitigating  effects 
on  both  wildlife  and  humans.  Multiple,  complementary  approaches  are 
necessary  to  connect  environmental  contaminant  exposures  with  biological 
effects.  Progress  has  been  slow,  in  particular  in  terms  of  wildlife  health,  and  the 
"weight  of  evidence"  approach  {Ross  2000)  has  been  the  most  powerful,  but 
weakly  effective,  means  of  producing  a  change  in  industrial  practices  and  human 
behavior.  Some  attempts  to  link  environmental  exposure  with  systematic 
responses  have  begun  to  emerge  recently,  specifically  addressing  the  health  of 
marine  mammals  (Wilson  et  al.  2005).  However,  much  remains  to  be  done. 
Projects  involving  pinnipeds  or  small  cetaceans  that  can  be  repeatedly  captured 
alive  and  sampled  extensively  (for  example,  the  HERA  Project)  may  provide  the 
best  information  on  individual,  population,  and  temporal  factors  important  for 
assessing  pollutant  exposures  and  risks  to  marine  mammals.  In  the  absence  of 
such  data  sets  for  endangered  species  and  large  whales,  the  use  of  in  vitro 
systems  and  model  organisms  is  a  more  accessible  means  of  predicting  or 
inferring  the  biological  effects  of  environmental  contaminants.  The  work 
presented  in  this  thesis  takes  advantage  of  available  molecular  and  biochemical 
techniques,  as  well  as  tools  from  model  species,  to  learn  more  about  marine 
mammal  toxicology  through  the  examination  of  the  aryl  hydrocarbon  receptor 
function  and  diversity. 

1.  Major  Conclusions 

1.1  Integument  biopsies  from  North  Atlantic  right  whale  may  be  used  to  clone 
genes  of  toxicological  and  physiological  significance. 

We  first  showed  that  the  application  of  two  established  techniques  - 
minimally  invasive  integument  sampling  from  living  whales,  and  cDNA  cloning  - 
could  be  combined  to  successfully  amplify  fragments  of  several  genes  of 
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toxicological  and  physiological  significance  from  the  North  Atlantic  right  whale. 
While  these  methods  have  existed  for  some  time,  the  use  of  this  tissue  source  for 
RNA  extraction  had  not  been  performed  (Lapseritis  and  Hahn  2004).  From  this 
RNA  we  acquired  initial  cDNA  sequence  information  for  AHR  and  CYP1 A1 .  The 
quality  of  the  biopsy-derived  RNA  was  compared  with  that  from  stranded 
cetaceans  at  various  time  points  post-mortem;  the  former  was  of  higher  quality 
and  yield,  but  the  latter  was  also  useful  for  amplifying  these  sequences,  including 
a  sequence  for  leptin.  The  use  of  biopsied  tissues  to  identify  expressed  genes 
opens  a  new  avenue  of  research  to  explore  the  physiological  status  of  living 
whales. 

1.2  First  characterization  of  aryl  hydrocarbon  receptor  in  mysticetes. 

The  initial  fragments  of  AHR  reported  in  Chapter  2  ultimately  yielded  full- 
length  cDNA  sequences  for  North  Atlantic  right  whale  and  humpback  whale 
(Chapter  3).  These  sequences  were  ligated  in  expression  vectors  for  a  variety 
of  in  vitro  assays  to  determine  biochemical  functions  relative  to  previously 
characterized  beluga  whale,  human,  and  mouse  AHRs.  Through  radiolabeled 
ligand  binding  assays,  we  showed  that  NARW  and  HW  AHRs  specifically  bind 
dioxin.  The  relative  binding  affinities  of  cetacean  AHRs  are  intermediate  to  those 
of  the  high-affinity  mouse  AHR  and  the  low-affinity  human  AHR.  The  biological 
activity  of  whale  AHRs  was  demonstrated  by  reporter  gene  activation  in 
transfected  mammalian  cells.  While  TCDD-dependent  activation  was  not 
observed,  some  ligand  specificity  was  suggested  by  differential  reporter-gene 
activation  by  benzo[a]pyrene. 

1.3  Ligand  binding  domain  sequences  may  be  used  to  infer  binding 
characteristics  in  mammalian  AHRs. 

While  many  researchers  have  investigated  the  effects  of  different  ligands 
on  the  binding  and  function  of  individual  AHRs,  few  have  compared  the 
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differential  behavior  of  multiple  AHRs.  Here,  we  compared  the  sequences  of 
AHRs  from  many  mammalian  species  to  identify  areas  of  similarities  and  unique 
sites  that  might  suggest  their  function  (Chapter  4).  We  contributed  new  partial 
sequences  for  mink  and  polar  bear  AHR  cDNAs.  Preliminary  results  from  a 
ligand  binding  assay  provided  relative  binding  affinities  for  marine  mamma!  and 
rodent  species  in  comparison  with  human  AHR,  and  provide  a  context  for 
predicting  the  sensitivity  of  these  AHRs  to  dioxin  in  vitro.  From  these  analyses 
we  suggested  a  set  of  amino  acids  within  the  LBD  that  may  be  critical  for 
differences  in  dioxin  binding  and  subsequent  AHR  function.  Homology  models 
were  produced  for  AHR  LBDs  that  illustrate  the  similarities  and  differences  in 
tertiary  structure  that  may  help  to  infer  binding  characteristics  attributable  to 
specific  amino  acids.  Some  of  these  residues  appear  more  important  than  others 
in  determining  the  shape  and  tertiary  structure  of  the  LBD  in  homology  models  of 
these  AHRs. 

2.  Future  Research 

At  this  stage  it  is  almost  ironic  that  there  remains  so  much  more  left  to 
learn  about  the  topics  addressed  in  the  six  years  of  work  that  were  put  into  this 
thesis.  However,  in  reality  it  seems  only  natural  that  a  thesis  should  inspire  at 
least  as  much  future  work  as  it  entailed.  As  has  been  said  many  times,  “the 
more  we  know,  the  more  we  know  we  don’t  know."  !  know  so  little!  And  the  list 
of  projects  to  continue  is  extensive. 

In  sum,  the  work  contained  in  this  thesis  contributes  many  details  of  the 
mechanism  of  AHR  action  in  marine  mammals.  However,  a  link  between 
environmental  contaminants,  AHR  function,  and  reproductive  effects  remains 
elusive.  This  is  necessary  to  connect  north  Atlantic  right  whale  (NARW) 
reproductive  failure  with  an  AHR-dependent  (or  -independent)  mechanism  of 
action.  Several  current  lines  of  research  may  shed  light  on  this  relationship. 
Among  these  are  efforts  to  better  measure  the  environmental  contaminants 
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present  in  baleen  whale  diets  (Moore  et  al.  in  preparation),  to  sample  NARW 
steroid  hormone  levels  (Hunt  et  al.  2006;  Rolland  et  al.  2005),  to  culture  cetacean 
oocytes  (Iwayama  et  al.  2005),  and  to  characterize  whale  steroid  hormone 
receptors  and  their  affinities  for  hormones  versus  the  contaminants  measured  in 
their  diets  (see  below). 

First,  I  plan  to  aggressively  pursue  amplifying  sequences  for  NARW 
estrogen  receptors.  This  is  an  important  goal  as  other  techniques  for  addressing 
endocrine  function  in  cetaceans  are  emerging  (Hunt  etal.  2006;  Rolland  et  al. 
2005),  and  the  development  of  in  vitro  tools  for  characterizing  reproductive 
function  will  complement  hormone  profiles  collected  in  the  field.  New,  more 
degenerate  primers  have  been  designed  to  amplify  either  the  alpha  or  beta  forms 
of  the  ER,  and  will  increase  our  chances  of  acquiring  sequence  information  if 
transcripts  are  present.  We  will  apply  these  primers  both  to  existing  RNA 
preparations  as  well  as  to  fresh  tissues.  We  hope  to  obtain  samples  of  deep 
blubber  and  liver  archived  at  the  NIST  from  a  NARW  calf  that  was  necropsied 
rapidly  post-mortem  in  January  2006.  We  will  also  collect  any  new  tissue  that 
becomes  available  during  the  2007  calving  season  when  most  vessel  strikes 
claim  the  lives  of  calves  and  females  near  shore.  These  same  deep  blubber 
tissues  will  also  be  used  to  continue  efforts  at  leptin  characterization  in  NARW. 

In  order  to  complete  the  work  of  Chapter  3,  we  intend  to  perform  two 
additional  experiments.  The  initial  transactivation  experiments  (Chapter  3,  Figure 
9)  must  be  confirmed  with  additional  replicates,  and  have  shown  that  10  ng  of 
whale  AHR  construct  are  likely  to  produce  the  most  unambiguous  results.  The 
surprising  differences  in  responses  to  B[a]P  compared  with  TCDD  in  a  single 
transactivation  experiment  allude  to  an  interesting  tale  of  two  compounds  with 
different  activities  in  related  cetacean  species,  We  will  continue  to  investigate 
this  observation. 

The  second  experiment  contributing  to  the  results  in  Chapter  3  will  relate 
in  vitro  biochemical  activity  of  NARW  AHR  to  specific  PCBs  and  PAHs  measured 
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in  blubber  and  NARW  prey.  A  filter-based  hydroxyapatite  assay  has  been 
designed  to  measure  competitive  binding  of  5  PAHs  (phenanthrene, 
naphthalene,  dibenzofuran,  benzo[a]pyrene,  and  anthracene)  and  5  PCBs  {PCB: 
1 01 , 1 1 8,  1 26, 1 53,  and  1 70)  shown  to  be  common  in  copepods  collected  from 
Georges  Bank  to  Cape  Cod  Bay  and  the  Gulf  of  Maine  {Moore  et  al.  in 
preparation;  Weisbrod  et  al.  2000).  We  will  also  include  astaxanthin,  a  potential 
AHR  agonist  naturally  present  in  copepods.  This  experiment  will  more  directly 
address  the  hypothesis  that  NARW  AHR  is  more  sensitive  to  ligands  more 
prevalent  in  its  ecological  niche. 

An  additional  set  of  ligand  binding  assays  is  required  to  improve  the 
results  given  in  Chapter  4.  Following  the  subcloning  of  hamster  AHR  in  the 
pcDNA  vector,  only  a  single  binding  assay  was  conducted  with  a  single  high 
dose  of  [3H]TCDD,  In  order  to  more  fully  develop  the  relative  binding  affinities  for 
marine  mammal  and  model  rodent  AHRs,  several  replicate  saturation  binding 
assays  will  be  conducted.  These  assays  will  help  account  for  differences  in  AHR 
expression  observed  in  the  quantitative  IVTT  that  influence  maximal  binding 
capacity  (Chapter  4  Figure  6).  In  the  absence  of  x-ray  crystallography-based 
determination  of  AHR  LBD  tertiary  structure,  the  homology  models  must  suffice 
to  illustrate  potential  structural  consequences  of  subtle  differences  in  this  region. 
Improvements  can  be  made,  however  to  refine  these  models,  perhaps 
incorporating  or  comparing  models  generated  from  other  templates  such  as 
ARNT  {Erbel  et  al.  2003)  or  FixL  (Pandini  and  Bonati  2005).  Virtual  mutagenesis 
experiments  with  these  models  will  also  benefit  the  overall  analysis  of  the 
importance  of  individual  amino  acids  on  AHR  ligand  binding  function. 

Those  last  experiments  would  be  greatly  improved  with  the  addition  of 
mink  and  polar  bear  AHR  sequences  and  constructs.  A  tremendous  amount  of 
effort  was  invested  towards  this  end,  with  many  pages  in  notebooks  as  the  only 
tangible  outcome.  Completion  of  the  LBD  sequence  would  at  least  permit 
homology  modeling  and  sequence  analysis  to  infer  binding  characteristics. 
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Efforts  at  increasing  the  sequence  data  for  these  two  species  may  only  yield 
results  if  new  RNA  is  prepared  from  fresh  tissues. 

These  continuing  projects  promise  to  provide  a  wealth  of  research 
opportunities  as  well  as  to  maintain  my  collaboration  with  the  Hahn  laboratory  at 
Woods  Hole  Oceanographic  Institution  for  several  years  to  come.  More 
importantly,  the  work  submitted  in  this  thesis  has  been,  and  continues  to  be,  the 
passionate  pursuit  of  something  greater.  Amidst  the  terrible  environmental 
circumstances  that  humans  have  created,  I  remain  hopeful  that  this  and  any 
future  work,  will  contribute  to  the  improvement  of  life  for  all  creatures  on  this 
planet. 
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